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SYNOPSIS 


In the present v^ork, we have studied low-lying 

energy levels of light end inteimediate nuclei (6 _< A _< 92) 

with realistic as well as phenomenological interactions. The 

thesis consists mainly of four chapters, first chapter 

gives a review of recent developments in nuclear theory, 

particularly the structure aspect. In Chapter II we first 

present a review of conclusions as axrived at by different 

researchers who ha^e made microscopic calculations of 

renormaliaiation corrections to free N-dl interaction. Purely 

microscopic approach turning out a failure, they suggest 

the introduction of phenomenological and/or experimental 

information into microscopic calculations. Taking this 

clue, we start with Sussex matrix elepients (SI'®) as effective 

interaction in first order. Pnder the assumption that 

the potential which the S^® represent is momentum independent 

and is a combination of central, spin-orbit and tensor parts, 

we first calculate the spectra of '^li, ^^Sc, '"^^Ca, 

58 . 92 

ITi and Zr nuclei with each of these components and with 
the combinations central + spin - orbit, central + tensor, 
and also with the total SME. We conclude that the spin- 
orbit part is almost ineffective for two-body spectra. 

However, the tensor part is more important, particularly 

for T = 0 states, and. finer details of the snectra are dependent 



XVI 


on presence or otherwise o” the tensor part in the inter- 
ection although the central p=>rt is the most important and 
dominant contributor to level energies. The calculated 
level energies are not in agreement with ex-cerimental 
values. The off-diagonal tensor force only effec- 

tively renorraali'^es the triplet central interaction, the 

3 

T = 0 low-lying states are predominantly rel'^tive states 

and the experimental values of matrix elements therein also 
bear the maximum uncertainties, 'v'e find that fevv of* the 
low-lying T = 0 levels of ^li, and ^^Sc nuclei can be 

reproduced if the relative matrix elements are multi- 

plied by a suitable factor ( 1 + a), i.e., a]l the renormali- 
zation corrections to the T = 0 SAIE can be effectively 

3 

simulated by a times the bare Interaction matrix elements in 
relative states. The calculated ground state T = 1 level 
in ^^Ni and ^^Zr nuclei turns out to be above 

the experimental position. A.ddition o-p a constant strength 
pairint-: force is shown to reproduce the correct ground state 
energy with an appropriate choice of the strength. 4n 
improvement over t'lis is the .addition o-^ a 6-function force 
(instead of the pairing force) which besides reproducing 
correct ground state energy, also improves some of the ’ 

higher levels. Further addition of a quadrupole force of | 

suitable strength is found to simulate the renormalization j 

r 

tf' e long range part as needed for hiaher levels. In : 

I 

[' 

! 
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general, J-dependent effective interactions can be established, 
for both 1=0 and T = 1 levels, through introduction of 
phenomenology. In Chapt-^r III,Hhe problems of and 

■TO 

Cl nuclei are investigated with phenomenological effective 

interactions of Oaussian radial shape and containing central, 

spin-orbit and tensor dependences .y The forms of effective 

potential tried are ceatral, central + spin - orbit, central + 

tensor and central + spin-orbit + tensor. The parameters 

2 

are deduced by making a fit to experimental level energies 
AO 

of ' K. Validity of the para'ieters is checked by applying 
them to energy level calculation for the 01 nucleus. The 
lowest 3 state is not correctlv described in both and 

'Z Q 

Cl nuclei. Except for this, the predicted energy levels 
38 

in Cl are in good agreement with experimental energies. 

In Chapter IV, ^we raaksy/a general study of the two-nucleon 
interaction. We assume a Yukav/a shaped central force of 
arbitrary strength and range parameters and calculate the 
effect of configuration mixing on T=0 and T=1 levels of 

configuration in Sc as the space dimensionality 
is increased by adding confisruration-after-configuration. 

Within fp-shell we calculate the splitting between centroids 
of T = 1 and T = 0 levels of the (Ofyy' 2 ) configuration as a 
function of the various strength parameters and the 

space dimensionality. /For a randomly chosen set of parameters. 
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we find th^t th .9 levels under consideration do not 
approach a saturated energy value as the space is expanded 
and sometimes the inclusion of configuration is found 

to have larger effect than that of few of f'le other low- 
Iving configurations. This puts suspicion over effective 
interaction derived within fp-shell since it appears that 
the single particle energy of an orbit is not the only 
decisive factor in determining its effect on lom’'-lving 
energy levels and it might not be an honest attempt to 
par-ametrize an effective interaction wherein the effects 
of the raore-important-than-the-truncated-model-sToace- 
conf igurations are absorbed by a forcible fitting of 
parameters. We justify it by approximately calculating 
the configuration mixing co ntributions to 0 and 1 (0fr^y2) 
configuration leyels, of all configurations lying up to 
0i^^^2 single particle level'. 



CHAPTER I 


IITTRODUCTION 


1.1 DEVEIOPM51\TT 0'^ FORM !IL ISM FOR RTJCLEAR STRUCTURE 
CALCULATIOUS : 

Underst finding the structure of finite nuclei 
requires the understanding of the basic- nucleon-nucleon 
(N-N) force. The nuclear forces between two nucleons 
are believed to be mediated by exchange of mesons which, 
hence, involve finite tim^e for interaction to take place. 
The first assumption that is made is that ever^hough the 
interaction might be meson-mediated but we assume it to 
be instantaneous i.e. we assume potential interactions. 
Then,- the many-body Hamiltonian for a system of A. par'bicles 
is obtained as the sum of individual kinetic energies and 
sum of all interactions 


H = T + Y (1.1) 

A 

where T = S T. (1.2) 

i=1 ^ 


One, then, solves the many-body Schrodinger 


equation 
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H y = E y (1.3) 

for the system's total energy E correspon'^ing to the 
A-particle v/ave function 


In principle, the problem is solved if we Icnow 
the mathematical form of V and the method of exactly 
solving the maiiy-body Schrodinger ea nation. Since none 
of these is available, we have to resort to approximate 
methods. We postulate models which might in some v/ay reduce 
troubles. 

In view of strong interactions, liquid drop 
model was the first to be proposed which satisfactorily 
explained the bulk properties of nuclei e.g. binding 
energies. One important conclusion that was derived on 
this model V'/as that for most part of it the nuclear inter- 
action is two-body in character i.e. the total potential 
energy can be written as 

A 

V = Z (1.5) 

K.j 

However, the existence of magic nmbers demanded 
a model wherein the nucleons might be considered as moving 
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almoBi; independently of each other so as to retain 
their individuality in order to be able to fill quantised 
orbits. Thus came the independeat parf'^le model (IR^) 
with infinite harmonic oscillator well and a one-body spin- 
orbit potential constituting the average field. However, 
this model could not successfully account for magnebic 
moments and failed completely in describing excited states, 
particularly in case of nuclei ’which have more than one 
particle in last partially filled orbit. It w^s, therefore, 
assumed that the two-nucleon interaction is not completely 
exhausted in giviny rise to the average potential and a 
part of it, the residual interaction, still operates 
be’bwe.en valence nucleons. The IllCP with residual inter- 
actions was called 'the shell model'. Thus, the Hamiltonian’ 
(1.1) can be rewritten as 

H = T + U + (7 - TJ) 

= + H' (1.6) 


where H is the average potential , 

Hq = T + U (1.7) 

is the sum of single particle Hamiltonians and 


H' = y - u 


( 1 . 8 ) 
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is the residufil interaction. 

Under the assumption of. weak; interactions, the 
many- nucleon wave function (1.4) can be written as 
product of w.'^^'ve functions ;Lor indiTidual n-^rticles 

y = (r^) 4^2 ^^2^ ^^*9) 

Dividing the problem into « core part and the 
valence nucleons part, the wave function (1.9) can be 
regrouped as 


y = H^c ^v 

where 4^^ is the product of single particle wave functions 

of V valence nucleons and 49 is the wave function for 

o 

(A-v) core nucleons. Then, the total energv g can be 
written .as a sum of core and valence nucleons' energies, 
ihe interaction between valence nucleons and core being 
included in the energy of valence nucleons. Tf the core 
is not disturbed while considering different levels and one 
is interested in splittings of energy levels rather than 
the absolute level energies, tlie core energy simply drops 
out. Thus, finally, the equation to be solved is 


(Hq + Y)^ = 


( 1 . 11 ) 



5 


where the vnriables pertein to T-^leiice nucleons only 
end the subscript v will not be put henceforth. We h=ve 
else sliffhtlj^ ch-nged the not^^tion.in thet now V denotes 
the residu^’l interaction between the valence nucleons 
(insteed of PI' as earlier). 

Starting with the lowest partially filled orbit, the 
unperturbed (i.e. no residual interactions) wave functions 
span an ini'inite dimensional Hilbert space such that 

^0 ^i = ®i ‘t^i i = 1, 2, . . . . oo (1.12) 

where 6^ is the sum of single particle energies of 
valence nucleons in the configuration The problem, 

thus, finally reduces to that of understanding the potential 
V such that eq. (1.1 1) when solved with V gi\res eigen- 
values ID which are in agreement with experimental level 
energies and wave functions V^'s which also properly 
describe other nuclear properties o.g. electromagnetic 
transition rates etc. 

While in a stationary orbit, the nucleons, as a 
result of residual interaction, can undergo virtual 
excitations for ■ time limits as prescribed by uncertainty 
principle. These virtual excitations , although experimentally 
undetectable, do affect xhe energies. If the virtual 
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excj.tntions '^re considered, "biie wn^re functions for 
Valence nucleons can no longer be t^icen ns purs (|)^'s, 
A.lso , beforehand tliei'e is no way of knowing how many 
orbits ere mixed. Hence, the wave function, in presence 
of residual interactions, h-'^s to be written ns 

CXD 

f = s n h 0.13) 

i=1 ^ 

The eo, (1.11') is then equivalentlv written as 


I n.(€^ 6^. + < (j) |7l (f). > ) = Ea. 


3 = 1 , 2 ,.... 


(1.14) 


Solving th? set of infinite number of enuntions 
(1.14) is Gnuive.lent to di^gonalisinv the matrix of 
(Hj^ + V) in the infinite diraensionnl Hilbert space spanned 
bv th'-i basis vooto-”s If Y is vsr"’’ weak, one could 

try to solve tho- problam hj retainin'? onl.y a finite number 
of terms in the expansion ( 1 . 15 ). Fow'ever, such hope 
fails to niaterislize if Y induces strong short range 
correlations. In presence of strong repulsions, the 
nucleons can scattvor into any orbit and any number of 
times. It would be of no use finding some way of diagonal- 
izing infinite dimensional matrices since the perturbation 
theory v/ould not be applicable for calculating matrix 
elements of Y. 
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In presence of singularities in V, one uses 
reaction metrir derived from V, instead of using V 
directlv, 

GC'S) =7 + 1 G- (E) (1.15) 

E - 

where G-(E) is the reaction matrix operator, Q is the 
Pauli operator allowing interacting particles to scatter 
into unoccupied states only and is the iinperturbed 
Hamiltonian acting in the space of intermedia’ te states. 

Once the two-nucleon potential is chosen - it may 

be a smooth phenomenological potential or a reaction 

matrix derived from free interaction containing hard 

core, the next step is to solve the self-consistent 

H-^rtree-Pock (HP) problem for single particle energies 

and v/ave functions which are then used in solving eq . (1.11) 

with the potential Y between valence nucleons same as the 

i 

original two-nucleon potential. Such attempts have sho?7n 
that the results calculated with HP single particle wave- 
functions are sim.ilar to those obtained v/ith harmonic 
oscillator vfave ' functions. Thus, instead of solving the 
self-consistency problem, one generally assumes the 
single particle orbitals as pure harmonic oscillator 
wave functions and takes single particle energies from 
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experimental spectra of neighbouring odd-mass, onc-valenco~ 
nucleon nuclei. 

If Y contains singjulerities , ue c'^n avoid them 
in diagon-ali.aation by introducing reaction matrix in eq. 
but, then, G- has to be calculated first using eg. (1.15) where 
one again faces Y. 

,4 method, called the separation method, to handle 

the singularities in Y, was developed bv Mosskowshi and 
2 

Scott . In this method, the total potential Y is divided 
into a short range part Vg, which contains all singularity 
•and a long range part Yj^, which is smooth and weak, such 
that when tw^o nucleons interact,. Y^ j.tself does not 
produce anv phase shift i.e. if the two interacting 
nucleons do not come close enough so that be felt, the 

presence of is -as good as its absence. Thus, the 

residual interaction is ef fectival 3 '' just the long range 
part Y^. Yet another m.athod of handling singularities 

3 

is the unitaraz-modGl-ouerator approach of Shakin et.al. 
where one does not calculate reaction matrix. Instead, 
they introduce correlated wave functions which vanish 
v\^henever two nucleons come closer than hard core radius. 

The calculation of matrix elements of a hard core 
interaction becomes a two-step process: one has to first 
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construct the potential end then derive the reaction 

matrix. Experimentally, all information about the 

interaction is contained in the pha-^'f. sbifts and, in 

principle, it should bo possible to relate the matrix 

elements of the interaction to the measured phase shifts, 

thus shippin-^ the intermediate step of constructing 

a potential. Such an approach has successfully been 

attempted by Elliot et.al.^' emploving auxiliary potential 

method. The matrix elements obtained are generally 

called Sussex matrix elements (SJ/IE), The analogy of this 

method to Moszkowski-Scott separation method for deriving 

reaction matrix elem^ents from phase shifts has also been 
S 6 

established . Sanderson et.al. have also suggested a 
way of incorporating short range repulsion in the original 
SBIE. 

Once the single particle basis, single particle 
energies and the two-nucleon interaction is specified, 
ono should construct the Hamiltonian m.atrix in the infinite 
dimensional Hilbert space and diagonalize it to get the 
eigen-energies and eigen-vectors for the valence nucleons 
wherefrom other nuclear properties can be calculated. 
Hovifever, infinite dimensionai. space is urmianageable . 
Therefore, one looks for a formalism wherein one may 
work in a manageably large subspace of the full Hilbert 
space and still get the same result as one v/ould get in the 
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infini'tft (? imensiorLal sp=ice. This process of selectin.Q; 
out only few configurations to worh with is called 
truncation the truncated space ys nodel space. Let 
us denote the full space br T) and the model space hj d. 

An operator 9 that describes a nuclear phenomenon in I) 
will not do so in the same way if calculations are performed 
in d only. Hence the need of an 'effective’ operator 6 

ext 

such that calculations with ©g^-P iii yield same results 
as one would get with 0 in space D. 

Let us rewrite the basic enjdations in I) ’ 

00 

= I a. (|) 

1 

Ho 't'i = ®i i = 1 oo (1.16; 

H = + V 

H = E '-W 

lie 'define' the model space weve function 




(1.17) 


Let us define two pro.-jection operators: 

P operating on^p projects on to the space d: Q operating 
on projects onto rest of the Hilbert space. Mathe- 


matically, 
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P = S li > < i ! (1.18) 

led 

Q = I !i> < 1 I 

Q 

However, P + Q = I, tbe ir1enti^7 operetor. 

Thus , 

Piy> = 1 7a > 

We 'define ’ the effective Hamiltonian through 

-W 7d = E fd 

We write from the set of equations (1,16) 

■ S a. (H^ + 7) |). = E S a. (j). 
iGD ^ ^ ^ iSD ^ ^ 

Multiplying on left v;-ith ({). and integrating 

c' 

over all space 

(S - e.) a. = 2 a. < (t). lY| > 

J D qeD ^ 

= < ([)^.lv| j s P (1.20) 

t,* 


Althoiigh this equation is only a rest'pitement of 
the original eigenvalue problem but it is important in this 
form in that it shows how one can do without both V and , 
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Let us define 


Then 


Tlr> = I ^ 

, ea. (1.20) cen be inxerpreted to imply 


( 1.21 ) 


(Ho + Vf’ 

Compeirison with sp. (1.10) '^ives 


( 1 . 22 ) 


®eff “ "'^eff 


problem of calculating reduces to the one 

f calculating 7gpf , the 'effective interaction’. 


1 . e 
01 


Let us see, how, using eq. (1.2l), can be 

obtained from t. 

= V| vy > = 7(P + 0) 1 'V> 

= t| %> + '/o|iV> 


IIov;, 


Qiy > = S H>j > < $3 I 'll 1 'I>1 > 

i BL 

Z K. > a^ < (f). 1 (!)i > 


iSD 


= Z I (t)i > a. 
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?_ (using eq. (1-20)), 


■i^d 




. 0 


E-H, 


vlH^> = 


0 


Y , 


E-Hq ef f 


> 


Hence , 


Q 


V 


TeffI % > = Vl ^ -eff 

0 


''Vd > 


or i in the operntor notation 


VePf = V + V 


Q 


E-H, 


V 


eff 


(1.23) 


i.e. given a two -nucleon potential Y, the effective 

potential ^Qff is obtained through an integral equation. 

We can see that Y^^^ explicitly depends on energy E which 

it is supposed to reproduce i.e. the effective interaction 

is state dependent. The state dependence may also include 

dependence on local density, centre-of-mass quantum number c 
7 

etc . 


let us compare eq. (1.23) with eq. (1.15). The 
operator Q in eq. (1.15) aliov/s scattering into unoccupied 
particle states only. In eq . (1.23), the operator 0 is 
further restricted to allow the particles to be scattered 
outside the model space only. If one chooses only the 
lowest configuration as the model space, the effective 
interaction Y'^^^ is completely equivalent to the reaction 
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matrix operator G- If model space consistyof more than 
one conf iffuration, to take into account configuration 
mixing, one has to diagonalize the matrix of effective 
Hamiltonian, Then, in eq. (1,15), Q can be thought of 
as made from two parts : one part allows scattering into 
model space and whose effect is treated by diagonalizing 
the reaction matrix in model space. However, the reaction 
mcctrix so formed within the model space should also 
incorporate the effects of tte other part of 0 which 
permits scattering outside the model space only. This 
process of modification of reaction matrix, calculated 
within model space, by happenings outside the model space, 
is called 'renormalization'. One thus, sees th'^t, in 
lowest order, when renormalization effects are ignored, th. 
reaction matrix also represents effective interaction. 

Even otherwise if one exactly treats all the renormali- 
zations, the G matrix is equivalent to the effective 
interaction. The unrenormalized G matrix is also often 
called the 'bare' interaction. Thus, one may use realistic 
interactions for spectroscopy calculations in finite 
nuclei provided one uses a, reaction matrix and makes due 
allowance for renormalization effects. 

Starting with a free intera.ction, Kuo and 

Q 

Brown were the first to calcula.te the lowest-order 
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renormalization of reaction matrix. They obtained 
impressive agreement with experiment al level energies 

‘9 

using the renormalized interaction. Barrett and JCirson 

then calculated the third-order contribution of core- 

polarisation process and they found that it almost 

cancelled the second-order contribution of ’ Kuo and Brown. 
1 0 

Kirson has finally shown that infinite-orders sum of 
perturbation series for all hitherto phvsicallv conceivable 
renormalization processes is zero. i¥e present detailed 
review of renormalization calculations in Chapter II. 

Thus, starting from free N-F potential one m.ay 
get good agreement with experimental data (as happened witL 
Kuo and Brown) but one can not be certain if the reaction 
matrix has been correctly calculated and renormalizations 
properly treated although it is a matter of satisfaction 
that one knows the physical origin of the two-nucleon 
effective potential. In view of these uncertainties, one 
sometimes adopts a reverse approach. One knows a good 
deal about the qualitative character of K-K interaction » On 
the basis of this one sets up a potential e.g. a 
potential of the form 

Vgff = (a + b 

incorporates the spin exchange character of the force. 
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The function f(r) determines the redial dependence of the 
potential and it will also contain the range of the force 
which may be taken to be spin-dependent. The parameters 
(a, b and force range aopearing in f(r) in above example) 
of such a potential are determined by requj.ring that 
experimentally observed properties of the nucleus under 
consideration, e.g. its energy levels. are reproduced. 
Sometimes one also adopts Talmi approach wherein the two— 
body matrix elements themselves are parameterized. Effective 
potentials derived in this way are termed 'phenomenological'. 

I. 2 problems ?0R present work : 

(i) Except for a few unavoidable practical limit- 
a.tions, the 'N’-R force in free space is very well known, 
Hov/ever, it can be realized at the outset that a free 
N-M interaction will not also describe the situation 
inside a bound nucleus for the obvious reason that the 
free R-tT force would be modified by presence of other 
nucleons. In case of phenomenological effective potentials 
the core part is assumed to be inactive and parameters 
are forced to reproduce experimental numbers. In case of 
free interaction we do not have anv free parameter at 

our disposal which might be for>ced to attain some value 
so as to reproduce experimental facts. Therefore, with 
a free N-R interaction between valence nucleons, thpir 
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interaction with core nucleons can not be ignored. 

Further, an infinity of final states is arailable when 
two ps^rticles Bc=(tter in free space. Wif^in bound 
nucleus, they move in quantised orbits,. Yet, while 
initially in a given q;uantum state, they can scatter into 
any of higher unoccupied orbits. In principle, they are 
also infinite in number. Due to our inability to handle 
infinite dimensional matrices we have to define a model 
space* Again, parameters can be forced to get a phenomeno- 
logical effective interaction but with realistic inter- 
actions, some analytical mathematical formalism has to be 
developed to incorporate both core effects and effects of 

outside-of-model-space configurations i.e. the free N-N 

in 

interaction needs to be properly renormalized. Kirson 

has finally shown the inadequacy of present microscopic 

calculations for estimation of renormalization corrections. 

1 2 

Schucan and Weidenmuier believe th^t the failure of 
purely microscopic approach lies in tne facr that in this 
approach the theorv is completely isolated ^^rom 
experime.atal results for bound states, except that the 
free H-F interaction is also most generally derived to 
fit to bound state properties of deuteron, and at no 
stage of calculation do we know where our calculations are 
directing us to. However, this doss not mean that things 
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pre wrong with theorv, of course, there are always 
practicpl limitations, but only that the true renorm- 
alization process has yet to be discovered. So long as 
we do not know that process and mathematical formalism 
to handle it, we can at least deduce in a much simpler 
way the amount of renormalization needed by introducing 
phenomenological and experimental information into calculat- 
ion so that at least the qualitative nature of tie 
renormalizing force may be understood although it is 
rather definite that once we come down to phenomenology, 
there can be more than one ways of achieving the same final 
effect on mathematical numbers and wre do not claim the 
uniqueness of phenomenology we plan to introduce. 

it least one thing is obvious from the results obtained 

by so many people wording with realistic interactions 

that the derivation of reaction matrix itself, starting 

from a fre<= N-IT potential, is one of major causes of 

uncertainties. In veiw. of this, the Sussex matrix 

elements (Slffi) present a favourable situation in that 

the reaction matrix needs not be derived. Previous 
1 S 

calculations have shown that bare SME are not sufficient 
to describe bound state nuclear properties i.e. there is 
a need to renormalize them. However, we do not adopt a 
microscopic approach for calculation of renormalization 
corrections. The SME are available as numbers and the 
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exact form of the free N-N poteutiel which they represent 
is not known. If we were to adopt microscopic approach for 
renormalisation it would be only numbers for relative 
matrix elements which wouild be needed but to decide for 
the phenomenology to be adopted, the behaviour of inter- 
action towards different nuclear properties should be 
known. 4ssuming the potential to be smooth • and a 
combination of central, spin-orbit and tensor parts, we, 
therefore, first decompose the SME into matrix elements 
of component interactions and study the effect of each 
component on low-energy spectra of two -valence-nucleon 
nuclei ^Li, ''®0, "^^Sc , ^^Ca, and ^^Zr. 

Comparison of these spectra .among themselves and of each 
with experimental -one will help in assessing the importance 
of individual components. It will be helpful in deciding 
the phenomenology to be adopted. This way, we will end 
up with an, we can call it semi-phenomenological or semi- 
realistic effective interaction and tbo approach lies 
somewhere between purely microscopic approach to calculation 
of renormalization corrections and a totally phenomeno- 
logical parameterization of effective two-nucleon 
Interaction - let us call it semi-microscopic. We discuss 
this problem in Chapter II. 

(ii) The next problem is that of the isotones 
and ^®C1. Today, we ha/e enough experimental data 
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regarding levels of Paucity of enough experimental 

data has always been realized as the biggest block in the 
way of establishing a reasonable phenomenological effective 
interaction. Moreover, generally even the low-lying 
nuclear levels contain quite complicated components in 
the wave function so that parameterization of a phenomeno- 
logical effective interaction in a model space of simple 
configurations is always doubtful. In view of these 
difficulties, the nucleus presents quite favourable 

a situation since it had beeniealized as back as in 1956 

1 4 

bv Goldstein and Talmi that low-lying states of this 

nucleus are almost pure configurations. Moreover, about 

9-10 low-lying levels of have well established 

energies and- spins and therefore parameterization of a 

phenomenological effective interaction comprising s^of 

central, spin-orbit and tensor dependences can be adventured 

with much less stringent approximations. However, the 

40 

available experimental information for K is gust 

sufficient for determination of parameters of an effective 

interaction consisting of a central, a two-body spin-orbit 

and a tensor part. To check their validity, the most 

58 

immediately calculable thing is the spectrum, of Cl 

15 

through the use of Pandya transformation . We discuss 
this problem in Chapter III. 
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(iii) In the above mentioned two problems we 
attack the problem of establishi;Tg effective interaction 
from two different corners. There are few yet undecided 
factors which decide the nature of a true interaction 

may it be either a phenomenological effective interaction 
or the one derived from a realistic two-nucleon interaction. 
However, the work of Cohen et.al.^^ on pseudonium nuclei 
showed that one could adopt a smaller model space and yet 
obtain an effective interaction by simply an adjustment of 
parameters that the 'pseudo experimental ' data are 

17 

excellently reproduced. On the other- hand, G-upta and Trainor 

proved that this result had no generality and nuclear 

properties dependent on wave functions e.g. electromagnetic 

^ ■ 

transition rates could still|relied upon as sensitive test 
to 'goodness' of an effective interaction. 


12 

Schucan and Weideum u.ller have shown that, a 
priori, there is no reason to expect that the perturbation 
expansion for "^Qff converges. Honconvergence means that 
we are trying to calculate too strong effects in perturb- 
ation theory. We renormalize a free N-H interaction to 
take into account configurations neglected from being 
included in the model space. Nonconvergence would ther. 

imply that such strongly contributing configurations 
should be explicitly included in the model space. This 
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is not the realization of facts for the first time but 
it has always been avoided on the plea of otherwise 
increased computational difficulties. We study this 
problem of truncation of configuration sp^ce for ^^Sc 
nucleus for a randomly chosen two-body central potential, 
however, keeping in mind that the approximate numerical 
results given by this set are not orders of magnitude 
off the experimental results so that the qualitative 
results we plan': to extract do not turn out to be comple- 
tely hypothetical. We first study the separation between 

2 

centroids of T = 1 and T = 0 states arising from (Ofyy' 2 ) 
configuration as a function o^’ the four A|ps strength 
parameters and space dimensionality. This is done with 
a view to assess the relative importance of even and odd 
^ partial waves in nuclear structure calculatiO'ns . This 

study, is expected to be useful since, at times, the use 

18 19 

of a Berber force or just a central, S-state Interaction 

has been tried for calculation o-p two-bodv spectra. This 

study will then tell us the limit in which such an 

approximation can be used that we altogether drop the odd 

part of the interaction. 

There is no a priori reason to believe that all 
configurations within a model space contribute equally and 
a particular configuration within the usually chosen 
model space might be contributing just negligibly in 
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comparison to other configurations within the same 
model space. This study is expected to be useful for 
a more honest calculation of shell model effectiye inter- 
action in the s ense that if some of the usually chosen 
potential parameters and configurations can be established 
to contribute much less than others then they can as well 
be dropped out and an equal number of parameters and 
conf igurat ions can be included in the calculation, which 
would otherwise have to be excluded in view of computational 
difficulties and which can be proved to have more pronounced 
effect - qust without adding to computer time. We study 
this problem in Chapter IV. 



CH4PTER II 


PHEHOMEWOLOCr l C.AJj 5STIMATIOF OP 

TO PRES RUOLEOR-MCLEON iMSRAg PlOI 


II. 1 IFTRODUCTIOR; 

Before one malces nuclear spectroscopy calculations, 
one has to decide for the effective interaction one wants 
to work with. There are two choices open-phenomenological 
effective interaction or an effective interaction derived 
from free nucleon-nucleon (R-R) interaction. Each one 
has its own merits and demerits. 

Suppose we choose some phenomenological form. 

We may not be sure whether we ha\re inclnded all the 
required features of a true R-R potential. Suppose, we 
assume a central force. Then, we are ignoring, say, 
spin-orbit and tensor forces, whose effects, we may not 
know beforehand, might be important. Suppose we include 
them also. Then, we might still be missing some other 
features, say hard core. In principle, we can incorporate 
as many features as we wish while formulating an inter- 
action but we have to remember that each time we include 
a new component, we add at least one more parameter ci,nd 
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tbe limit to number of parameters, which can be determined, 
is set by the amount of available information. We have 
seen in Chapter I that one mostly finds the unknown 
parameters by making fit to experimental excitation ener- 
gies and one, at least, requires as many excitation 
energies as the number of parameters. Such a calculation 
is feasib].e only if the needed energy leyels also have 
well defined spin and parity. Our knowledge of experimental 
energy levels is not rich enough and not many lev^^ls 
have unambiguously established spin and parity, even in 
the low energy region. This means that one can gain only 
a limited description of the effective potential. 

Even if enough experimental data is available, it 
is very important in these calculations that one should 
be able to assign a reasonable configuration to an experi- 
mental level so that it may be meaningful to compare the 
matrix element calculated for that configuration, with 
the level energy. It often happens that, for example 
in two-valence-nucleon nuclei one observes a very low- 

lying O"*", T = 1 excited level which cannot be explained 

20 

on the two-particle shell model « Instead, it has 
been shown to be a core-deformed four particle -two hole 
(4p - 2h) state . With a core-excited state getting 
down so much low in energy, one may legitimately 
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question the purity of other nearby levels in respect of 
shall model configuration. In view of such ambiguities, 
it is no surprise that for same nucleus one often ends 
up with two completel 3 '' different sets of parameters 
which reproduce the energy’" levels equally well. 4n 
example is found in ref. 22 v/here it is shown that different 
sets of parameters give equally good fit to doublet splitt- 
ing in .4 = 14 and 4=16 nuclei. 

To avoid these difficulties in constructing a 
phenomenological form for potential, one may decide to 
adopt the Talmi method. Again, the paucity of enough 
experimental information comes in the v/ay. However, one 
has a way out. One picks up a set of nuclei in the same 
shell e.g. one might select all nuclei with A = 18 to 
20 for determination of effective interaction in 
sd-shell^^. Matrix elements for configurations of more 
than two valence nucleons can be expressed as line=ir 
combinations of two-body matrix elements. One makes the 
assumption that the two-body matrix elements do not 
change in magnitude in neighbouring mass number nuclei 
of the same shell. Thus, one parameterizes the two-body 
matrix elements simultaneously for a series of nuclei. 
However, such calculations have to be restricted to a very 
small model space just because one mav not be able to 
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assign reasonable configurations for ell the energy 
levels even if all of them have well established spin 
and parity. One is thus constrained to assume that 
configurations lying outside the model sp'^ce are either 
ineffective or their effect has been absorbed in the 
calculated values o^ the two-body matrix elements. Unless 
one establishes that effects of left out configurations 
are negligible, one can only infer that, alright, the 
matrix elements represent those of an ef'f’ective interaction 
but it may not be fair to say that numbers so obtained 
are representative of 'true' U-N for>ce. Obviouslv, such 
calculations are strongly space dependent, further, one 
is making a total neglect of many-body effects. One is 
assuming that the two-body interaction in a many-valence 
nucleon conf iguration is not altered by presence of other 
nucleons from a value it is in a tv/o-valence nucleon 
system, and that three-or more-bod.y forces are not present. 
Suppose, we had enough knowledge so that even in this 
approach we were able to find the effective interaction 
separately for different nuclei. Then, if the two- 
nucleon interaction is altered by. presence of otl® r 
nucleons or if more-than-two-bodv forces themselves are 
important, it will be reflected in the calculated values 
of the e.ffective two-body matrix elements. 
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One makes, a simplification over this method. 

Overy two-body matrix element, diagonal, off-diagonal 
or of any configuration, is basically a linear combination 
of relative matrix elements. The relative matrix elements 
themselves stay constant in all two-body matrix elements 
for a given nucleus; geometry of the matrix element only 
alters the coefficients of the aforementioned expansion. 
Therefore, one tries to parameterize the relative matrix 
elements ’ . Within ca given ma;jor shell, only a finite 

number of relative matrix elements appear. This way, one 
requires much less experimental data than in Talmi approach. 

This approach is further simplified by one's 
knowledge about the nuclear fordes. One knows that 
nuclear forces are short ranged. Hence, it may be a 
good approximation to ignore interaction in sta.tes with 
X, > 2 since in higher partial waves the nucleons are not 
allowed to get close enough due to increased centrifugal 
repulsion. One may further assume the existence of 
central forces only in which case, for example, the 

T IZ 'Z 

and ^^2 relative states are degenerate. 

People have gone upto the extent of neglecting even P- and 
D - partial waves. Interaction in P-states can be ignored 
on the plea that forces are much weaker in odd-1 states. 
Neglect of D-state interaction is justified by saying 
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that it v/ill be much weaker than the leading S-state 
term and that once we are all out to determine an 
effective interaction only, let their small effect also 

oe absorbed in the effective interaction perameters. 

5 3 

One should, at least, include the strong 
tensor force but this force only renormali'^es the triplet 
state interaction and so it is also included in the 
triplet-state central interaction. Like this, one is able 
to reduce the number of free parameters drastically but 
only at the expense of losing details about the structure 
of the U-N potential. 

After the relative matrix elements are thus 
determined, one may attempt at interpreting them in 
terms of an effective potential and find the strength and 
range of the potential. 

Above discussion makes it clear that an effective 
interaction, phenomenological atleast, cannot be unambi- 
guously established. In fact, the numbers depend upon 
which configurations one is considering and which 
nuclear properties one wants to reproduce. Since most of 
the attempts are aimed at reproduciug energy levels, wave 
functio.is are never guaranteed and hence no guarantee 
that wave functions will be good enough so as to reproduce 
other nuclear properties also equally well. 
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Due to difficulties in determination and uncerta- 
inties in results obtained for an effective interaction, 
one looks for some alternative aoproach wherein the 
potential appears as a natural conseauence rather than 
a forcible introduction to reproduce experimental data. 

One has a good knowledge of interaction in free space; 
one tries to empl.oy it for calculation of bound state 
properties of nuclei. 

One is renuired to calculate matrix ele'^ients of 
the interaction. Therefore, if the interaction possesses 
a hard core, one should first of all calculate the reaction 
matrix 


Q. = Y + V .2 & (2.1) 

E - 

However, it should be remembered that a reaction matrix 
which is representative of interaction in free space 
will not describe the force inside a bound nucleus 
also. Probably the best example in support of this 
statement is: in free spoce, neutron has a half-life of 
about 12 minutes but inside a bound nucleus it stays 
permanently as neutron! Certainly the interaction of a 
neutron with other particles when inside a nucleus is 
different from the interaction with same particles in 
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free space. To the extent that core could be 

considered as inert, the reaction matrix itself will 

be effective interaction in infinite dimensional 

0 

Hilbert space of valence nucl>- on s. However, the idea of 
infinite dimensional space is not practical and further , 
the concept of effective interaction itself is associated 
with a model space. Let us see what form does eq. (2.1) 
take when the full Hilbert space is partitioned into 
two subspaces. Let, 0=0.^ + ^2 where 0^ allows 
scattering into model space only and Q 2 allows them to 
be scattered into any not-already-occupied state. Assum- 
ing the outside-of-the-model-space effects to be small 
enough so that retention of only those terms where 0^ 
appears onl 3 /- once is a good approximation, the expansion 
on right of eq. (2.1) can be grouped to give 

Qo 

+ G. = G. (2.2) 

E - H 

0 

where G^ is the reaction matrix generated within the 
model space 

^1 

g = Y + Y — G. (2.3) 

' S - ' 

However, a complete expansion gives 
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Qp 

& = G (2.4) 

E - H 

0 

This means that if one wants to work in a 
model space, one will have to iterabe over the reaction 
matrix calculated in the model space. How many higher 
orders will be required in order to get congruence with 
experiment will depend on how good is the description 
rendered bv G^ itself. 

In the present wDrk we have adopted a semi- 
microscopic approach to derive effective interaction. 

What persuaded us to take such a step has been the fact 
as observed by others that even with the present day 
computers that application of more and more sophisticated 
calculational techniques has been possible, we do not yet 
have a satisfactory theory for nuclear structure physics. 
Our attempt has been to see if v/e can understand the 
H-H force in a simpler way. To make a case for ourselves 
(in defence of our simple minded approach over much more 
involved cal-culations of others) it will be in order that 
we review as to .ultimately what has been the fate of an 
almost uncountable nimiber of calculations employing purely 
microscopic approaches. 
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Historically, first attempt of calculating 

reaction matrix for an open shell finite nucleus was 

m»de by Dawson et.al^^ They used the HIT potent ial^^ 

18 

and made calculations for 0 nucleus by diagonalising 
the reliction matrix in sd-shell. The ordering of first 
five states was correctly reproduced and the level 
energies were roughly reproduced. In the lov/est approx- 
imation when the reaction matrix was approximated by the 
interaction itself, it w as clearly seen that on expanding 
the space f rom (Od^/^, to (0d^^2’ ^<^3/2^’ 

the results turned out to be in better agreement with 
experiment. The results further improved vfhen the reaction 
matrix also included 2 'fid) intermediate states. Although 
the gross features one could conclude were that the ground 
state was predicted at a higher energy than experimental 
value and the calculated spectrum vras too compressed, 
yet the results pointed towards a brighter future for 
such calculations provided one mahes a more careful 
calculation of the reaction matrix and assumes a large 
enough model space. 

It is for one's own convenience that, there not be 
a mess of computation that one :-essumes the core as inert. 
Results of Dawson et.al. proved that the reaction matrix 
is not weak enough to retain particles within the model 
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space only. Besides being scattered into outside the 

model space, they could also interact indirectly - via the 

mechanism of core excitation. This process represents a 

long-range correlation which is produced by short-range 

interactions. Classically it can be viewed lihe this: 

One of the two nucleons, particularly onihe surface of 

a large nucleus, interacts with a core-nucleon and excites 

it to some empty orbit , thus creating a particle-hole 

pair. Then, the other valence nucleon interacts with 

the excited nucleon and throws it back in the created 

hole. In terms of diagrams, the process, in which two 

nucleons in the state | ab >, while interacting via 

core-excitation, are scattered into the state |cd >, is 

represented as shown in Figure 2.1a. In this diagram, 

a,b,c and d are particle states within the model space. 

Clearly, the intermediate states in this process are 

those of three particles and one hole (3p-1h). This 

correction to the reaction matrix is denoted by 

It can be noted that to conserve parity, the core-nucleon 

can only be excited by (n = 0,1,2, ) 

single-particle energy. The special cose of n=0 occurs 

when the core is formed by shell closure in the middle 

58 

of a major shell e.g. the pg Fi nucleus. The core 
56 

is Hi with both neutron and proton shells closing 
at Cf y/ 2 ■'^‘^’isnce neutrons being in the IP5/2 




Ffg. 24f. Infinite orders sum of T DA and 
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orbit. The parity reauiremert does not forbid a 
nucleon to be lifted to any of the orbits ^^5/2 

and "1 p -j 2 * 


Let us see, mathematically how the contribution 
of such a physical process can be calculated. In the 
original derivation of reaction matrix, eq,. (2.1), the 
Pauli operator Q was designed to allow scattering into 
particle states only. In fact, the introduction of this 
operator is necessitated by the requirement that the 
particles not be allowed to be scattered into occupied 
states. The idea of core-excitation was not conceived 
at that time. The above decomposition of 0 into 0.| and 
^2 ^nd eventual arriving at eq . (2.4) should be valid if, 
apart from two-particle intermediate states, 0.2 also 
incorporates 'other;- renorm-alisation processes which do 
not violate Pauli principle e.g. the core-excitation 
process. Henceforth we drop the subscript 1 on G-^ with 
an understanding that the reaction matrix is being 
calculated in a model space. Then, the lowest-order 
correction, second-order in 3-, to the bare 3 matrix will 
be 


0 


3 


3p1li 


3 


5p1h 


3 


where e is the appropriate energy denominator and 



57 


ph 

27 

It was first of all pointed out by Bertsch 
that the core-polarisation process has the useful 
property of lowering energy levels of low J and raising 
those of high J. It is just the desired ctfectrthe 
calculated ground state will be in better -u; reement with 
the experiment ant the compressed spectrum will spread 
Importance of core-excitatio i process is indicative of 
importance of dynamic effects of the core after the 
static effects have been eliminated by factoring out the 
core energy* 


Similarly, the correction arising from two- 
particle intermediate states (ligure 2.1b) is calculated 


as 



G 



G 


Qoc = ^ |2p > < 2p I 

2p 

In Figure 2.1b, both p^ and Pg are particle 

states outside the model space. Ov/ing to the intermediate- 

state energy occurring in the denominator, it is expedted 

that, except for the particular case where intermediate 

state excitations into same major shell are possible 
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like in M, intermediate states with s ingle particle 
excitation energy 2’tnCO will be most important. 



We have considered the process where one core- 


nucleon is raised through 2T-\6J energy. However, ere rgy- 
wise it should be equally probable that two core nucleons 
from just below the Fermi surface are lifted, each by 
I'KiO, into the model space. In this process, intermediate 
states are of four particle-two hole (4p2h) character. 

This correction is represe.ited as in Figure 2.1c and 
denoted by 0 - 2 ]^ • 

Thus, the total, renormalized interaction becomes 


G + G-^ 


3p1h 


G2p + G 


2h 


'Kuo and Brown^ applied these ideas to structure 
18 18 

calculations of 0 and P nuclei starting from Hamada- 

28 

Johnston (HJ ) potential • They first calculated the 
spectra with bare G interaction. The results were nowhere 
in agreement with experimeat thus confirming the necessity 
of renormalizations. However, they included only the 
G^p^h correction. They did not have to go beyond lowest 
order correction (second-order in G). In the second-order 
itself, the agreement with experimental spectra was 
satisfactory. It was natural to thinh that third - and 
higher-order corrections all cancel amoung themselves. 


Kuo and Brown neglected "the G 2 p correction on the 
grounds that it is essentially a second-order ladder 
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diagram and hence its effects are already included in 
the calculation of reaction ma.trix^ 

However, energy levels are a crude test of the 
effective interaction and they, infact, failed to reproduce 
E2 transition rate from second O’*" state to first 2"*” 
state. These states are usually believed to be deformed, 
containing 4p-2h core-excited component. However, they 
were able to reproduce reasonably well, energies of these 
levels as two particle states with 3P'^k'. corrections. The 
plea that they gave for neglect of ^2I-' correction was 
that the 4p-2h deformed states are highly collective and 
hence are much lowered in energy and so it will not be 
proper to calculate the G- 2 ]^ correction in perturbation 
theory taking a 2 energy denominator. This neglect 
thus showed up in the transition rate. 

The Kuo-Brown G- matrices were recalculated by 
29 

Kuo using the same intermediate states assumptions 
but treating the short-range repulsion of the interaction 
and the state dependence of the second-order tensor term 
more carefully. His results for G- and were 

significantly different from those of Kuo and Brown 
and this destroyed the good agreement with experiment 
they had earlier. Kuo, therefore, incorporated the G2p 
and ^ 2h corrections also and this improved the fit. 
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Since the first successful celculation of 
core- polarisation correction by Kuo and Brown, there 
have been a large number of calculat ions^^"^'^ starting 
from a free H-K interaction. These calculations have been 
performed with either 3- + ov G + + G2p ^2h’ 

They all substantiate the original findings of Bertsch that 
the bare interaction does not give satisfactory agreement 
with experiment and even philosophically it is wrong to 
make comparison with experiment at this stage. The 
ground state energy is not correctly reproduced and the 
ca.lculated spectrum is too compressed. However, inclusion 
of core-polarisation lowers the calculated ground state 
and also causes an overall spreading in the sprectrum thus 
bringing it in better agreement with e xperiment , 

Thus, although reasonable agreements with experi- 
nent were obtained but from a theoretical viewpoint it 
was a disturbing feature that the theory did not provide 
a unioue prescription for calculation of reaction matrix 
and for number of terms to be kept. For ex.arople, Kuo and 
Brown got impressive agreement with G + G^p.|^ but when 
Kuo did an improved calculation of G, he had also to 
include G 2 p and (x 2 j^ corrections to restore the agreement. 

Even Kuo and Brown had realised that it might not 
be sufficient to include just one-bubble term in the treatment 



of core-excitation. Instead, one should calculate it 
as a series of forward going bubbles (known as TDA.) 
shown in Figure 2.1d or, as a se 3 oies of both, forward 
and backward going bubbles (known as RP\) shown in 
Figure 2.1 e. 

Thus, the total interaction should be as shown 
in Figure 2.1 f. v^here the hatched circle on left neans 
a summation of series of diagrams on right. I. series of 
particle-hole pairs created one after the otter constitutes 
what is knovm as a 'phonon'. The reason whv the inter- 
action amopg valence nucleons should be renormalized by 
a phonon and not by a single bubble is that some of 
these phonons may be highly collective and thus much 

lowered in energy and then t aking 2itOJ as energy denominator 

8 

would be much undermining its contribution. They argued, 
however, that some of these partic].e-hole states will be 
spurious and they should come at zero energy. Further 
there is every chance that several of tte se phonons will 
lie much higher in energi’’. Hence, they thought it 
legitimate to i.nclude just one bubble correction. 

is more and more core-polarisation calculations 
were made, it became clear that the first-order corrections 
(second-order in G-) had a large effect on spectra and 
in some cases the second-order correction was even larger 
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than the bare Interaction G- itself. One, then, becomes 
suspicious whether or not are the third-and higher-order 
corrections actually negligible. 

in attempt to incorporate the strings of bubbles 
introduced above was first made by Kuo'^ . He constructed 
the core-po lari ^sat ion phonons in RPl and included their 
effect by coupling them to the valence particles in 
i=18 nuclei. He observed that low-lying spectra calcul- 
ated ?7ith G + & + ^ 3 p-]h were quite similar, 

flence, he concluded that to calculate the core-polarisation 
effects involving exchange of positive parity phonons it 
is sufficiently accurate to calculate the effect in 
first-order theory itself i.e. include just the 
correction although calculation with phonons is desirable. 

An arithmetical error in Kuo's calculations was 

•7 ^ 

discovered by Osnes and '-Varke^ who repeated the RPA 
calculations for A = 18 nuclei and concluded that signifi- 
cantly different results for renormalization of the two- 
nucleon effective interaction are obtained depending on 
whether the core-excitations induced by the valence 
nucleons are represented by single unperturbed particle- 
hole excitations or by the full RPA solutions for the core. 

Since then a trend has set up for calculation of 
core excitation corrections to effective interaction 
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between valence nucleons in TDK and RPA. Zamick"^ 
concluded th^t the core-polarisation calculated in 
RPl is much ].ar^er than the one in second-order pertur- 
bation theory. Ellis and Siegel explicitly summed to 
second-and all-orders the perturbation theor*'?' diagrams 
of the TI)4 and RPA for renormalization of ef’^ective inter- 
action between two sd-shell valence nucleons. They 
demonstrated, the eauivalence of infinite sum. to the 
method of construction of phonons and observed large 
enhancements compared to the first-order perturbation 

correction. This equivalence has been demonstrated 

39 

again by Goode and Eirson . Starting with Eallio-Eollt- 
veit (KE) interaction^*^, Siegel and Zamick^^ have observed 
the same enhancement in respect of effective charge for 
closed shell + one nucleon nuclei. In fact, for certain 

38 

choices of G matrix elements, the enhancement is infinite 

42 

Eirson and Zamick have named this process 'propagator 

renormalizat.ton ’ . .All these calculations end up with. 

the conclusion that strong coll.ective enhancement occurs 

when these series are summed up and that, in particular, 

the ground states are depressed far below the experimental 

energies. It, then, became clear that some compensating 

effect was needed to off-set 'this negati's^'e enhancement; 
43 

Osnes et.al. have repeated these calculations for 
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Valence nucleons in fp-shell. They also find that the 
difference betVireen RPT and one-bubble descriptions is 
much larger for fp-sbell valence nucleons than for 
sd-shell ones. i possible explanation for this difference 
is that the core is much bigrer than the 0 core 

thus allowing a larger number of intermediate states. 

Bruecknerhs idea behind introduction of reaction 
matrix Mtas that after all singularities are treated out, 
the resultinLg G- matrix will be weak enough that meaningful 
calculations could be performed in a model space or at 
most one would req.ui'^e one or two lowest-order corrections. 
However, the inadeQuacy of bare G- interaction and large less 
of second— order corrections had necessitated investigation 
of tliird-and higher-order corrections. Therefore, at the 
same time that above mentioned HP'l corrections vvere 
started, Barrett and Klrson undertook the calculation of 
third-order correction for = 18 nuclei using Kuo 
reaction matrices derivefl from rJ potential. In a 
preliminary note^"^ i^hey observed th«t there is no evident 
term- by- terra convergence v/hen the core-polarisation is 
treated ia perturbation theory- Instead, an analogy 
with Betlie'’s nuclear matter results'^^ suggests that the 
core-polarisation is so strong that it cannot be treated 
by perturbation theory at all but must be handled by a 
method which sums perturbation theory to all orders. 
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They, then, reported^ the results of their calcuitions 

to their-order and also of a few selected terms to 

fourth-order which contribute to tho pertuibation 

series for 0'*', T = 1 effective interaction between valence 

nucleons in ‘i = 18 nuclei. They found the total third- 

order correction to be as big as or bigger than the 

second-order contribution. They, thus, concluded that 

there is no apparent convergence in the perturbation 

expansion for in powers of G. However, they did 

observe that the contribution of 4p - 2h c orrection in 

third-order was an order of magnitude smaller than the 

corresponding number in second-order. It v/as, thus, safe 

to conclude that no collective enhancement is produced 

by this tvpe of correction. This conclusion is consistent 

with that of Goode^^ and Rajewski and Ilirson^'^ who have 

summed this effect to all orders in G and found little 

48 

or no enhancement. Barrett and Srison also performed 
calculation for tho simplest third-order number conserving 
set of diagrams which represent the renormalization of the 
bare react j.on matrix for A. = 18 nuclei. The sum of terms 
within this set was found to be small, being 10 per cent 
of G or less. 

There is double-counting involved when Ggp 
correction is included in third-order. Although the exact 
amount of this double-counting is uncertain, it has been 
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shown by Kirson^^ that t-ho inclusion of two-particle 
ladders in an expansion in^erms of the Kuo G matrices 
creates a double counting which can be as big as 100 per cent. 
For this reason, Barrett and Kirson had not included 
the G2p correction in their third-order calculations. 
Barrett^'' has attempted to determine the size of this 
double counting by comparing the Kuo G matrix ele'^ents 
with the Barrett et.al.^^ (B®l) reaction matrix elements. 

He finds that the amouit of overcountins^ depends strongly 
on the size of the gap between occupied and unoccupied 
states in the single particle enccgy spectrum in the 
calculation of the BHM reaction matrices, being largest 
for a small gap. In spite of the oyercounting problem 
which occurs, it has been suggested that two— particle 
ladder terms should ba included in the expansion for 
computed with the Kuo G matrices. It nas also been 
argued^ ^ that the third-order RPl diagrams should be 
included eTsnthough these dia^grams contain excitations 
of 4-y^(UOand ere not a pune contribution. In view 

of these suggestions, Barrett^^ recalculated the tnird- 
order correction with two-p^^rticle ladder terms and RPl 
diagrams included. The results weakened his earlier 
conclusion regarding nonconvergence of the perturbation 
expansion for although he still maintained that 

the prospects of improyed convergence were questionable 
in view of oyercounxing problem. 
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In the same model space (sd-shell), G-oode and 

54 

Koltun have shown that the (JT) - weighted average of 

the perturbation theory terms is repulsive in third- as 

well as in second-order. Therefore, the conclusion of 

Barrett and Kirson that second-and third-order terms 

cancel for O"^, T = 1 levels cannot be generalized to other 

•JT states. Goode has performed calculations for A = 18 

nuclei for all JT states in sd-shell. He concludes that 

with the possible exception of O"*”, T = 1 levels, the 

calculated spectra are quite similar whether calculated 

in second-or in third-order perturbation theory although 

the spectra calculated with calculated upto second- 

or upto third-order terms are much different from the 

ones with bare reaction matrix. Grosslv, it can be 

said that third-order effects on the T = 1 states tend 

to be small compared to corresponding second-order 

O'f'fects. Howeyer, this is not quite the case with the 

T=0 states where oyidence for nonconvergent behaviour is 

found, probably due to manifestation of the tensor part 

o-p the F-N interaction. The T = 0 spectra calculated 

in first, second-and third-order perturbation theory 

are quite different. In particular, the most dramatic 

+ -h 

third-0 ‘?der effects occur for the lowest 1 and 2 , T = 0 
states each of which lies about 1 MeV higher when third- 
order correction is aoplied than when only second-order 
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effects are included. It appCcsrs th=it the two-lying 
experimental 1 = 0;- spectrum of is better represented 
in second-order perturbation theory than in third-order. 

Although Barrett and Eirson had proved a term-by- 
term nonconvergence of perturbation series in G-, yet 
they had suggested that it might be possible to rearrange 
the perturbation series into physically significant groups 
of diagrams so that each group converges order-by-order in 
G. However, one must then find some way to sum. all there 
groups of diagrams. Otherwise, as such, simple perturb- 
ation theory does not seem to be a reliable method 
for calculating the effective tv70-nucleon interaction 
when currently standard calculation techniques are used. 

42 

It Was suggested by Kirson and lamick that for 
calculating '^Qff between valence nucleons, one should also 
renormalize the interaction between the valence nucleons and 
the core. This was re-^erred to as vertex renormalization 
and it was suggested that it should also be summed to all 
orders. They calculated the modifier vertex to second- 
order in G, included it in the RPA and found that the 
resulting vertex renormalization tended to cancel the 
propagator renormalization. They also found that the 
multipole structure of the renormalization was more 
complicated than that given by a pairing or a pairing- 
plus-quadrupole interaction. 
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There is no reason why the pqrticle-hole inter- 
action should not be renormalised by the core as well. 
This process was called 'screening' of the particle-hole 
Interaction. Core renormalization of the particle-hole 
interaction with various potentials has been considered 
in refs. 56-61. Blomqvist and Kuo^ employed HJ inter- 
action for calculation of negative parity particle-hole 

40 43 

excitations of Ca and Ca v/ith and without including 
the core-oolarisation. lov^^est 3~states were found to be 
very sensitive to the core-polarisation. The RPA v/ith 
the bare G- interaction overestimates the collectivity 
of these 3-states. The inclusion of self-screening 
diagram largely reduced the collectivity and the final 
result was nuite similar to the TDl solution with bare 
G interaction. It must be kept in mind that or 

are big cores and allow many more particle-hole excit- 
ations than does a smaller core say ^^0. However, if 
one can e stablish the same fact-thai 4 renormalized RP.A ' 
solution is as good as 'bare TDl' one - for other cases 
as well, one has obtained a much simpler way of dealing 
with core-excitation phonons-TD\ with bare particle-hole 
interaction. Since the screening of the particle-hole 
interaction so greatly affects the collectivity of the 
core states, it must, in turn, be expected to be 
important for core-polarisation of the particle-particle 

- 4. 4- • 63 

interaction 



50 


1 0 

Kirson , then, undertook to sum the RPA series, 

the vertex renornializqtion series and the self-screening 

diagrams to all orders in G and found th^t the net 

result of these three infinite series taken together was 

approximately zero. The TD4 and. RP4 series greatly 

depress the spectrum. The Inclusion of screening 

damps the RP4 back to TDi. Pinallv, the inclusion .•f 

the vertex or black box renormalization produces a 

spectrum which looks like that for the bare G. Klrson^^ 

has rece.ntly incorporated one more effect neglected in 

his earlier calculations, namely, the phonon-rmediated 

screening of the vertex connecting zero-particle, zero- 

hole’ states of the core to two-particle, two-hole states. 

He concludes that core-polarisation, when self-consistently 

treated is too weak to account for significant renorma- 

1 6 

lizatlons required for the 0 core, whether applied to 

particle-particle •r hole-hole effective interaction 

or to E2 effective charges. Calculations shnilar to 

those of Kirson (however using different forces) were 

65 

recently performed by Jop;ko and Sprung for the fp-shell 
and similar conclusions arrived at. They conclude that 
the renormalizing effects found in the TD4 and particularlv 
the RPA were strongly damped by screenin^^ the particle- 
hole interaction. The results obtained for the screened 
TDA and RPA are auite similar and intermediate between 
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those obtained for the correction corresponding 

to describing the core states as unperturbed particl-e- 
hole excitation;! and the bare TDi. Kuo and Osnes ended 
up v/ith the s ame conclusioii for effective charge and 
also for effective interaction^"^ for sd-and fp-shell 
valence nucleons. 

Kirson's conclusion-the sum of three series 
taken together almost vanishes-is depressing since one 
knovrs that one does need non-e;ero renormalization. 

However, there are two points which offer hope. One is 
the G- matrix calculations of Barrett et.al.^ which 
indicate that the G- matrix can be made much more attract- 
ive by varying the intermediate states' spectrum, especially 
by decreasing the gap between the occupied and unoccupied 
states. An alternatiye method for computing G- matrix 
exactly and simultaneously diagonalizing it for two 

68 

nucleons outside an inert core was developed by lawson 

who solved the coupled iiitegro-dif f erential equation 

describing such a shell model problem. He also studied 

the dependence of G on the g ap between the occupied and 

unoccupied states and obtained results similar to those 

69 

found bv BHB/I and by McCarthy 

1 useful aid in calculating le that of 

shell model diagonalization in a model space of reasonable 
size (like loludice et.al. ' ). Such!f'''5a!^a^'tii^.x|^@f^ '^el'y 

i rghlTRAt LIBR^ I 

1 Ace. No. A 4 188 2] 
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lengthy and difficult to perform and are probablv 

impossible to carry out for heavy nuclei. But, in cases 

where such calculations can be perfo med, they mav 

indicate the correct magnitude of ^Qff Using these 

results, one may be able to find a simple expression 

or model for ? which can reproduce results similar 

to those of shell model diagonali’iiation. This method 

is, in essence, non-perturbative and involves diagonal- 

70 

izing bare G- inanenlarged space-Lo ludice et.al. 
performed a shell model calculation of the energy levels 
of 1=18 nuclei in the configuration space of two particles 
plus three particles and one hole using Kuo G matrices. 




I2p> 


+ 


2 

3p1h 


^3p1h 


3p1h> 


These results, when compared with first-order perturb- 
ation theory are found to be generally similar. This 

71 

calculation was repeated by Lo ludice et.al. however, 

using a different method of calculation. It turns out 

that the linked cluster folded expansion converges in G 

in a reasonably rapid way. However, the fact that they 

underestimate the effective charges and B(E2) values 

suggests that considerably more core-polarisation is 

required than they predict it to be. A similar attempt 

72 

has recently been made for 1=18 nuclei bv including all 
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the original Kuo and Brown calculations, na'nely, the 

core-polarisation acts like a pairing interaction and. 

depresses in energ*^ the low-lving states of low J. 

Barrett and Kirson^ also c.alcul^ted 

large third- order terms contributing to with 
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Clement and Bara-nger ^ matrix elements and the SMS and 

found that the results for the different potentials were 

sirall^=r in that the third-order contributions v\rere 

large in comparison v/ith the second-order contributions. 

The results of these calculations would indicate that t he 

convergence problems in trying to compute ^Qff by 

oi'dinar^.?- perturbation theory (i.e. the relative sizes 

of terms contributing to ^Qff) are independent of the 

form of the F-N potential as long as the potential fits 

the F-F scattering phase shifts data and bound state 

properties of the deuteron reasonabl 3 r well. It was, 

therefore, also felt that the differences in the off-shell 

properties of different potentials are .not important. 

Many of the G matrices which appear in second and 

higher-order terms are of i-the-energy-shell. A study 

of off-shell effects of this kind has been carried out 

77 

by Pradhan et.al. who find it is possible for these 
effects to be nuite large-fisore than 1 MeV for some low- 
lying states. A-gain, there is question of the physical 
validity of the phase shift equivalent G matrix elements 
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v/hich produce lerge off-shell ef'j^ects'^® . In any 

case, the influence of off-shell effects on higher-order 

terras needs to be studied in more detail- . Another 

difficultv concerning calculation of reaction matrix is 

the choice of the single- particle basis e.g. '^^’oods- Saxon 

IQ RO 

wa'^-T'e functions have be-'-n used' ^ for initial and 
final- states between which the matrix elements of G- are 
calculated and it was found that the G matrix elements 
for 1=18 nuclei are appreciablv more repulsire for a 
Woods-Saxon basis than for a harmonic oscillator basis 
causing upward shifts in the nergy spectra by as much as 
2 MeV. "i^et another difficultly in the calculation of react- 
ion matrix is the proper treatment of the tensor force 
which is not weak enough and tends to scatter nucleons 
into intermediate states of medium energy so that the 
G matrices which result f-rom potentials which contain 

the tensor force are still strong. Calculations by 
8 1 

Barrett indicate that it is the contribution of the 

tensor force to the G matrix which gives rise to the 

strong collective excitations of the nuclear core and 

that if the tensor forces are weaker, the contributions 

of these terms in the perturbation expansion for 

would be less important. Similar conclusions regarding 

68 

the spectra of A=6 nuclei were arrived at by Lawson 
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It should be noted thqt the calculations by doode 
and KirsonJ^ do not prove that tho perturbation series 
for converges. Their calculations also do not 

constitute a surrmation O"^ this series since they do only 
partial summations. They only shov/ that when certain 
subsets 0-^ diagrams are sumiiied to all orders in G, one 
obtains a finite result-. One can argue tha u this procedure? 
being based on limited physical insight, might overlook 
i-inportant series of terms. Therefore, resul..s of this 
type, although certainly instructive, can hardly be 
considered conclusive. Schucan and eidenm J. Il^r have 
investigated the perturbation expansion j.or of .ective 
operators and algsbraically studied its convt.rg^-nce 
properties. They conclude that, a priori, there is no 
reason to expect that the perturbation expansion for 
T (or any effective operator) converges; it is likely 
to diverge or at least converge poorly. However, they 
also show that it is possible to r---rr •■ng... tne seri..K. 
algebraically so that convergence is improved. 


Unless a complete absence of raany-body forces i; 


established, the next logical step 
pair interaction exactly should be 
particle interaction exactly? then 
interaction, and so on - a cluster 


after treating the 
treating the three- 
the four-particle 
ap nr OP oh. Such an 
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approach, has had considerable success in the theory 
of infinite nuclear matter®^’^^’^"^. Clearly, the 
convergence of the cluster series vyill be related to the 
ratio of the 'interaction volume ' of a nucleon to its 
'occupation volume' l.e. tc tne quantitr (force range/ 
nucleon separation)^. For the short raogo part v^, 
this is of the order of 10 per cent. Thus, cluster tech— 
niciues are anpropriate for the treatcoent of v^ although 

85 

use of reaction matrix is Inevitable. in estimation 
of three-body forces in nuclear matter suggests that 
whereas the lowest-order three-body effect is small, 
second— order effects are larger with the net result of 
about 2.5 MeV per particle additional binding energy at 
saturation density. In an analvsis®^ o-f” effective inter- 
otion in p-shell nuclei, the importance of Feingold's 
three- body vector interaction arisins from second— orde^. 
effects of the tensor force has been demonstrated, "'/arhe 
and G-unye^^ have calculated binding energies and r.m.s. 
radii of nuclei in the mass region 16 < A. < 208 in the 
Jastrow cluster exp'=nsion method by including the second- 
order correction to energy using the heid soft core 
potential^^. They conclude th^t the contribution of the 
three-body cluster correction to the two-body correlation 
energy is at the most 9 per cent. Quesne^^ has generalized 
the method of Talmi for extracting two-body interactions 
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from experiment to include all possible many-body 

forces compatible with the chosen configurations. It 

is shown that the interactions of the order higher than 

two are important in 0 Osnes^^ has found 

that the contributions of effecti'/e three-body forces to 

the birLd'’ng energies oi calcium isotopes c'^lculated 

with G- matrix elements deriTed from HJ potential are 

92 

small. Eisenstein and Kirson have shown that if a 
phenomenological shell model fit to energy levels of 
Ofy^^p -nuclei is performed by using both 'effective' 
two-body as well as three-body forces, then the contri- 
bution due to the effective 1iiree-body forces is not 
negligible. In other words, the microscopic approach 
invoking core-polarisation reno-^malization will have to 
take into account these effective three-body forces 
as well, fact is th-'t the problem of nuclear effective 
interaction is a three-body problem-two valence nucleons 
and a core. In view of the fact that the ordinary per- 


turbation series for defined in terms of two-nucleon 

interaction and its reno ’'’mall sat ion by core excitation, 

95 

does not converge, Lovas and Vegh have attacked the 


problem from a three-body point of view. Other 
94-96 

Crclculations regarding the sise of many-body forces 


in the three-and more-body problem indicate that they 
can be signific.'^nt and hence cannot be neglected. Employing 
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procedure first sugsestad bj Brondow , B-^rrett et.al. 

constructed the Q-fective interaction in A=18 and ,\=19 

nuclei bv performing a direct truncation into a smaller 

model space usin?; th^ rcsultc of a large shell model 

Calculation. Iti Cases where th^ truncation is severe 

i.e, when there is small overlap between th.' large space 

and small space v'ave functions, they find that a large 

effective three-bodAr term must be added to the effective 

two-body interaction obtained for A=18 nuclei in-order 

to reproduce the results for A=19 nuclei. On the other 

90 

hand, the calcul-'tions of Quesne indicate th-'t for more 
than three particles, there is a partial cancellation between 
tho three-body contribution from all other higher-order 
terms, so that the oyerall many-body effect is reduced. 

Thus, there are many eff-ects which heye not been 

adoauately considered or could be treated more accurately 

in the calculation of G, "^nd energy denominator 

and which could possibly lead to reasonable agreement 

with expe'^lment within the framework of the present theory 

and available computational technioues. Schucan and 
1 2 

VToidenm u Her have proposed that the nonoonvergence 
of the perturbation series follows from thej fact that 
we have no control oyer which experimental results will 
be best reproduced by the theoretical calculations, for 
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example, the etates produced by niuy not correspond 

to the lowest experimental states i.e. there may be 

'intruder' experimental states arpona the calculated 

states. It then follows That Wt-j will, probably, be 

forced to ititroduce phenomenological anri/or experimental 

information into the microscopic theury to guide the 

calculations of the effective operators away ■from the 

physical effects v/hich are causing the expansion to 

di.verge. This point of vi.ew has also been advanced bv 

08-100 

Harvey and Khanna" . The introduction of such infor- 

mation is no more a downgrading of the microscopic 
theory than the introduction of G in place of V in order 
to treat the short range correlations. It is sim.ply 
a matter of letting nature guide the mathematics so that 
the mathematics will get to the physical answer in the 
moat economical manner. 


Without getting involved in the mathematics, let 

us see what nualitative features v/e have learned about 

renormalipiation corrections, in particular the core- 

27 

excitation. Even Bertsch , at the time of Introducifig 
the idea of oore-pola--isation, had demonstrated that the 
core-polaris'-^tion process tends to simulate a pairing 
force. \lthough the quadrupole component of the core- 
polarisation is the most important, the other even 
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multipoles contribute enough to give the overall force 

==■ pairing character. The nrgaement 'wes given originally 
3 1 

by Kuo . He considered the diogram (figure 2 .. la) 

wher'3 the particle-hole excited pair is coupled to (J'', T' 

Each vertex in the figure must be a scalar. There'^ore, 

if the particle and the hole are coupled to (J'', T'Os 

th.^m thti corresponding tensor component of the effective 

interaction ie picked out. This is becau.se the two 

valence nucleons can be visualised as interacting via 

wavy lines (the G- interaction) with a (J'', T'') coupled 

particlc-hole pair scueezed in between. Thus, for this 

diagram, the core-polarisation correction from a given 

J’’, T'' behaves as an irreducible tensor force of 

31 

rank J'',T''. Kuo has made calculations for nickel 

1 p ^^2 model space and finds that 
the J'' = 2, T'* = 0 component provides the largest 
contribution. It then means that next to pairing, the 
moat important component of the renormalizing force is 
the nuadrupole-quadrupol interaction between the valence 
nucleons. Mathematically it means a P 2 However, 

the arguement is weakened by the fact that above diagram 
represents only the direct term in the correction. In 
fact, one should include the exchange term also as shown 
in Figure 2.2b. In this diagram, we have only one vertex. 


isotopes in llP'^/2 ^^5/2 
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Hence, if the p-article enh the hole ere coupled to 
(J'', T'')> then only this component of the force 
enters the upper interaction but all allowed tensors 
enter the lower interaction. This me^ns that the (J'', T'' 
renormalizing component should not be solely viewed as 
one of correspond ing multipolarity. The correspondence 
is only approximate. 

Freed et.al. have e xploited the above arguements. 
They use Tabakin interaction as residual interaction in 
lead region thus avoiding handling the singularities. 

The single particle wave functions in t'vo-particle 
intermediate states are taken as plane waves. Their cont- 
ribution is calculated in second Born approximation. 

However, the core-excitation effects are calculated 
phenomenologically to simulate the qualitative features 
discussed above. This approach is thus intermediate 
between purely phenomenological and purely microscopic 
approaches. It is argued that the interaction of the 
valence nucleons through the particle-hole excitations 
of the core, to some approxim'^tion, can be incorporated 
through a first-order Pg foi'cs. Similarly, discrepancy 
between calculated and experimental ground state 0"^ 
levels in even-even nuclei may be removed by addition 
of a pairing force. The procedure for fitting the 
strength parameters of these renomalizing components 
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was that first the quadrurole force streriiffth was so 
adjusted that the first 2"^ level w-^s reproduced. Then 
the pairing force strength was so adjust od th^t in the 
presence of* quadrupole force of strength as determined 
above, the ground state energy is correctly repro- 
duced . 


In the present work we take the SM.T as effective 
interaction in first order. However, we do not calculate 

the second-and higher-order corrections microscopically. 

101-104 

Instead, like ?reed et.al. , we have also estimated 

the renormalization corrections on the basis of qualit- 
ative features of the renormalizing force as discussed 
above. Obviously we have had to introduce both phenomeno- 
logical and experimental information. However, our attempt 
has been to use least amount of experimental information. 
This particular choice of residual Interaction greatly 
conveniences the calculations since the numerical values 
of relative matrix elements are tabulated'^. We take the 
advantage of this in th-^t we can simul ^taneously perform 
calculations on several nuclei in different mass regions 
and thus attempt to explore the properties of efrective 
interaction, say, as a function of mass number. Although 
the specific choice of SI4S is purely for convenience, 
as discussed above, the results care expected to be general 
and be true for other interactions as well e.g, the 
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1 05 

ones with hard core like Yele , HJ etc. or the ones 
incorporating nonlocality like the one of Tpbqkin. To 
avoid ourselves being led into a mess of oomput etional 
difficulties, we consider nuclei with two valence 
nucleons only. We hope stopping at two nucleons to be 
sufficient because we are not considering manv-body 
forces and if only two-body forces are considered, even 
for a man'-^-valence .nucleon system, the interaction 
energy is a linear sum O"^ two— valence nucleon inter.actions. 
However, wc consider such nuclei where both T = 0 and 
T = 1 interactions Cvon be used. Thus, we studv the 
^Li, '’|p, ^qC-o, ^^Sc, IgITiand nuclei. We hope 

that if ?.'o can successfully establish the effective 
interaction for these nuclei, it will be usable for 
spectroscopic study of those manv-valence nucleons nuclei 
at least where the valenco nucleons are in the same 
orbit '-.s the lowest •='Ctlve orbit in these nuclei e.g. 
the ^0 40^^ nuclei etc. In the present 

work our emphasis has boon, only on to see that the two- 
particle interaction energy in the ground state and the 
excitation energies of at least a few low-lying levels 
are reproduced as best as possible although we under- 
stand that the task does not finish here. Evan if we 
can arrive at an interaction with which the calculated 
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level energies ogree oompletely with experinent , this 
does not guarantee thot our wave functions are also good ; 
thev might fail to reproduce other properties of levels 
e.g. transition rates. 


II. 2 METHOD OF CALCULITIOH: 


The SME are available as nianbers in a form 
directly usable in nucle'^r structure calculations. However, 
the analytical form of the free potential which they 
represent is unknown and our experience with other known 
free H-N potentials tells us that it will be complicated 
enough. We avoid introduction of complicated dependences 
like explicit dependence on momentum, density etc. Instead 
let us assume that the interaction is a scalar constructed 
from combination of spin and orbital degrees of freedom 
of a nucleon. Since the spin operators are tensors of 
rank one, we can obtain from <5^ and ^2 only a 
scalar ( vector e.g. ( 5^ x or 

a tensor of rank two. Hence, the total interaction, is 
a sura of central, spin-orbit nnd tensor parts characterized 
formalin by the rank k = 0,1, and 2 respectively in the 
scalar product ( ) where operates in 

coordinate space only and S^' '^operates in spin space only. 
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The 


total interactiorL 


if5 then obtained as 


V=E y(k) = J (yftO.jCk)) 

k 

Tht- SME are ^"^allable as matris elements of 

r between Yerioue channels of relative states In^^S+D 

However, by employing: -avoo. n j. -it. • j. • 

■ •• '^figular momeatum algebra, it is 

possible to deduce matniwr -i n. j? tt/i \ j:- 4 ., 

^rix elements of Y(k) from those 

of Y using 


<n' 1 ’ S ,l| V(k) I nLs ^ > 



( 2 k+l ) ‘ 


aZ (-)'^^tS+i'( 23 , + i) 

D' 


1' s ;i 




k 


<n'^'Sj' Ivjnf^S 

(2.6) 


> 


Eurtuer, except tensor force channels, 

tlie SME a’^o tabul^^tad onpy for n = n’, n' -f 1. However, 
in a harmonic oscillator basis, other matrix elements 
with (n-n’) > 1 mav be obtained from 


<n'|V|n> - |ii(n'-n-l) < n'-l|?ln> + '/"(n-i-'l ) (n+^ + 3/2) 


X.' < n’-llvln+1> + Y'lrTn+X+y757 <n< _i |vln-1> 


1) (n* - 1/r~-^ < n’ - 2 |vi n > 

’'^3 , ) / ( 2 . 7 ) 

I ^ >] flTTn’ + ^' + 'l/^J 


+ < n' -1 
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The potential is aesumel to be momentum independent 
so th'^t the last term in'/olvine commutator of r^ end V 
vanishes . 

'I complete knowledge O'f an effective interaction 
constructed as a sum of various components involves 
a comprehensive studv o-p the effects of* the constituent 
components on various nuclear properties and on nuclei 
in Various mass regions. If we adopt a totally microscopic 
approach, given a free two-nucleon interaction, in 
principle, the effective interaction is ob-'.ained as sum 
of infinite series. For a piecewise studv of effective 
interaction, a properly renormalised free M-ll force will 
be an ideal subg'^ct. However, if we want to avoid the 
use of microscopic approach for calculation of renormali- 
siatitn corrections (as discussed at length earlier in 
this chapter, the final outcome of a purelv microscopic 
apnronch indicates its futilitv) we shoul.d h^/e some 
idea as to the nature of the renormali3''‘tions needed. 

We know, philosophically it is wron^ to compare with 
experiment the results obtained with a free H-N interaction 
since we know there a need of reno rm.alia.ing the force, 
nevertheless, a study of the spectra with free H-H force 
and its central, spin-orbit and tensor components and 
their comparison among themselves and with experiment 
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v^ill enable -qs to draw conclusions about rel'-->tive 
inportaace of those components in nuclear structure 
oalcul-^t ions will thus provide us with a clue as to 
which part needs nior>e attention be paid to and thus guess 
ab-JUu thij needed renormalizations. 


T'hc SME are tabulated as a function of th:- 
oscillator length pe . , / f 


parameter b ( = ■ Using virial 

theorem, b can be related to the r.m.s 
nucleus 


radius of the 


b = 


N < r'" > 




( 2 . 8 ) 


where is the degeneracy factor i.e. the number of 

nucleons in an (n 1) orbit. Thus, for a completely filled 

.*a 

(n.^d) orbit, = (21 + 1). The number U has to be 
decided on physical grounds. If <r^> is the charge r.m.s. 
radius, U will be the total number of protons in the 
nucleus. If < ^ matter r.m.s. r-^^dius, N will 

be the total number of nucleons in the nucleus. The 
experimental values of r.m.s. radii have been taken from 
ref. 106 .- 


Once the b value is fixed, the appropriate 
rels'tivo matrix elements, picked from tha tables in ref. 4, 
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decomposed into cGntr.-^l, spin-orbit end tensor 
components by employing eq . (2.6). 

The required single particle energies are d educed 
from th-' spectra of neighbouring odd mass nuclei. 

The relative matrix elements are tabulated unto 
G— Wave. However, a study of the numbers reveals that 
the relative matrix elements decrease rapidly with 
increasing.^ . Therefore, we feel it sufficient to 
include upto 7=2 only for calculation of low-lying 
energy levels. 

Harmonic oscillator wave functions 
have been t ake-i a s single particle basis states. 

Calculations ha'Te been performed with the 
following interactions ; 

(1) Pull interaction, (ii) central part, 

(iii) spin-orbit part, (iv) tensor part, (v) central + 
spin-orbit parts and (vi) central + tensor parts. 

The two -body matrix elements in jj -coupling 
are calculated using the formula (1.29) given in appendix A. 

The respective model spaces for the di^^ferent 
nuclei have been chosen within single major shell only. 
Therefore., only diagonal-iu-n matrix elements of V, except 
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for the tensor force channel where n and n’ differ 
by 1 , ere needed. 

Experimentallv , excitation energies of levels 
arc Icnown. To make a comparison with calculated inter- 
action energies feasible, we first calculate the inter- 
action energy in thi ground state between the last tvi/o 
particles in a nucleus from known total nuclear 

bindlnff energies (B.H.) as follows. If both the particles 
outside the closed shell core are neutrons, the inter- 
action energy is 


B.E. ( ) + B.E. ( ) - 2 B.E. ( ^^2^) 

and if one is a neutron and the other a proton, the 
required number is 


B.E.( ^"^^X ) + B.E. (/,Y) - B.E. ) - B.E.( 


A 


■Z-1 


'\+ 1 - 

7 ,‘ 


A+1 

Z-1 


I) 


In calculating these numbers using these formulae, 
numerical values of the binding energies nre to be 
taken with a negative sign. The experimental binding 
energies have been taken from ref. 107. 
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Below we discuss configuration sp«ce, single 
perticle energies and b values for different nuclei. 


i ) Op-Shell : 

There is ambiguity regarding the value of 
^op ~ ®0p ^ ^ value for ^Li nucleus, 

M /2 ' ^3/2 

The electric charge ra-^ius tor ^ Li is approximately 2.5 

fm. Correspondingly, b, as defined abore, calcul.ates 

1 8 

out to 1.84635 fni ^¥hercas the value of b for "0 turns 
out to be 1.8466 fm. f’urther, the b value increases with 
mass number. As such, it appears that the calculated 
value 1,84638 fra is too high for a light nucleus like 
"Li. Also, the Stanford data' on elastic electron 
ncatterin.cr from ^Li show that the charge form factor for 

^Tii could not be interpreted either ir^erms of the 
conventional simple harmonic v»rell shell model as used 
for other p-shell nuclei or in teems o."‘ the modified 
haimionic well shell model in woich s ano. p nucleons 
are assumed to move in different harmonio wells. Ho 
acceptable fits to the data could be found .'.rom the 
cluster model approach with either of the two models viz. 
a + d and ^He + % structures, However, a fit to the 
data could be made with the phenomenological expression 



73 


for the charge form factor of which yields an 
r.m.s. radius to be 2.54 +0.5 fm. If this is so, 
then harmonic oscillator wave functions are not a suit- 
able comolete set that one could have chosen for this 
nucleus. Instead, one may look for the wave functions 
generat-^d in some other single particle potential. One 
calculation"''^^ yields the r.m.s. radius of ^Li as 2.52 
fm with a parameterized single particle local potential 
using two different values of the same parameter r^ of 
range for the s and p shells. 

110 

Lauritsen and A-.'jzenberg - Selove have not 

been able to assign an unambiguous value for AO- The 

Opwp level in ^He and ^Li appears much diffused over 

a v/ide energy region. A.ssimakopoulos et.al. have 

used a value 5 *2 MeY for A®. To get over this ambiguity, 

calculations have been made parameterizing the value 

of Talmi type phenomenological calculations of 

Cohen and Kurath"'^^ and Amit and Katz^^^ give A 0 to be 

1 1 4 

0.64 MeY and - 4.74 MeY respectivelv while Kumar 
reports a value of around 4 MeV for A.O. Dawson and 
ITalecka"’"'^, who calculated the energy levels of A =6 nuclei 
using the realistic hard core interactions, find that 
^0 MeY gives best agreement with experimental 

spectrum. Halbert et.al. however, obtain best fit 
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Z\ S - 3.2 MeV usin--^ the HJ potential. 


117 

It is possible to prove 
energy splitting between the 0 p ^^2 
particle levels is given by 


analytically that the 
and '7p-^y'2 


= I ^27 ( ^Pq) + 3Y (^P^) - 5V(^P2)J 

where the V's are the diagonal-in-n matrix elements of 
the potential in the respective relative states. Por 
Op shell, the energy index P (= 2n-j + 1^ + Png + I 2 ) is 2. 

Hence, in relative P - states, only the value zero is 
allowed for the relative total Quantum number n. In 
Tanle 2,1 we have tabulated the relevant relative matrix 
elements for b = 1.5 fm, 1-6 fra and 1.7 fra. 

TABLE 2.1 

s 

Biagonal-in-n (=0) SME in relative P-states as relevant to °Ii. 


b \ state 
(fm)\ 

5p 

^0 

(MeV) 

5p 

^1 

(MeV) 

3p 

^2 

(MeV) 

1.5 

- 2.57 

3.14 

- 1.85 

1.6 

- 2.46 

2.54 

-1.48 

1.7 

- 2.27 

2.08 

- 1.18 
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The value of As tu-rns out to be 5.0? ¥eY , 3.79 MeV 
and 2.85 Me 7 for b = 1.5 fn, 1.6 fm and 1.7 fm respect- 
ively. The basic philosophy behind use of realistic 
interactions is tbs elimination of free parameters and 
spectroscopy calculations starting from first principles 
as far as possible, ’’'^e see that the SIIC0 give the value 
A 6 ill range of values iust discussed as used by 
other authors for a b value in the range 1.5 fm - 1.7 fm. 
We, there "ore, use A 8 =3.72 MeV as used by Joshi"^^® and 
make calculations for three b values - 1.5 fm, 1.6 fm and 
1.7 fm. This f'^ay expect to be able to find the b 
value anpropriate to this nucleus. However, the energy 
of ground state A level will depend on the amount of 
configuration mixing which, in turn, depends on A 8. 
Hence, more precisely, the value o" b arrived at would 
be valid for A 6 == 3.72 MeV only. 

4t the time these calculation^ were performed 

only an old co!f‘pilat ion"^ of experimental levels of ^li 

119 

v/as available. They have been compiled again ' recently. 
The new energy levels differ little from old values. 
However, the j”^ of second level with T = 1 has been 
unambiguously established to be 2 . 
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ii) isOd - shell 

18 18 

For 0 and F nuclei, the configuration 
space is 

Od5/2 ^^ 3/2 


18 17 
?or 0, we can use single particle levels of 0. 

1 8 

However, for P it will be mors appropriate to tahe the 

17 17 

average of single particle energies in 0 and P 

. .120 
nuclei 

17o (MeV) 17p (MeV) 


5/2^ 

0 

5 / 2 * 

0 

1/2+ 

0.87081 

1 / 2 * 

0.49533 

3/2+ 

5.08300 

5 / 2 * 

5.10300 


Hence, 


the single particle 


energies for the 


two nuclei are 


Od5/2 

1Si/2 

0d3/2 


(MeV) 

0 

0.87081 


Od^/2 

^®l/2 

Od3/2 


(MeV) 

0 

0.68307 


5.08300 


5.09300 
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The experiniental spectra have been taken from 
ref, 121. 

\7e hqve used b = 1.S fm for both and 

i i i ) Oflp - s hell 

42 42 

For Sc and Ca nuclei the configuration 
space is 


0f7/2 

IP3/2 

1 P 1/2 

Of 3/2 

The single pa'i^ticle 

41 122 

energies in Ga 

and '^"’sc 

nuclei are 




^^Ca 

(Me7) 


(MeV) 

7/2" 

0 

7/2' 

0 

3/2" 

2.1 

3 / 2 " 

1.9 

1 / 2 " 

3.9 

1 / 2 ' 

4.1 

5/2" 

6.5 

5/2" 

6 . 4 
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Hence, the required values are 
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^2ca 

«Sc 


(Me'7). 

(MeV) 

Ofy/2 

0 

0 

IP3/2 

2. 1 

2.0 

1Pl/2 

5., 9 

4.0 

Ofj/g 

6.5 

6.45 


42 

The charge r.m.s. radius for Ca is 3.53 fm. 

Hence, b = 2.03795 fm. ¥e use b = 2.0 fni for both 

42 ^ , 42 „ 

^ Ca and ^"Sc. 

The experimental level spectra have been taken 
from ref. 124. 

cr o 

For ^ Ni nuc^-eus the configuration space is 

IP5/2 ^■^5/2 ^^1/2 

57 1 25 

The single particle energies in Hi are 



(MeV) 

3/2 ’ 

0 

(5/2)- 

0.76 

(1/2)' 

1.08 
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The ideiitificpition of 0 ^ 5/2 '^P'i /2 single 

perticle levels is doubtful. However, we use the follow- 

i “1 1 

ing values as also used by Joshi ' and bv Jain 



(MeV) 

1 P 3/2 

0 

Ofj/g 

0.78 

V1/2 

1.08 


58 

The charge r.m.s. radius for Ni is about 5-8 

fm and b turns out to be 2.0525 fm. However, for the 

whole nucleus we should use matter r.m.s. radius ratber 

58 

than charge r.m.s. radius. In Hi nucleus the croton 
shells close at level while neutrons fill 1P^/2 

level. Hence: the matter r.m-. s. radius is expected to be 

larger tlidn Charge Km.s. radius. We tahe b = 2.1 fm for 

58.T, 

Hr. 

The experimental energy levels have been taken from 
ref. 127 . 

iv) 2s1d0g - shell 
92 

Por Zr, active orbits are 

1 d 5/2 ^^1/2 


113/2 0 ^ 7/2 
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Joshi excludes the 0 grj ^2 level and assumes 
the single particle sequence as 


(Me?) 


OJ 

in 

0 

2Si/2 

1.22 


2.07 


1 22 

We call it set I. Cohen et.al. have given the single 

89 

particle energies of neutron levels in 53Sr nucleus as 



(Me?) 

115/2 

0 

2s ^/2 

1.138 

113/2 

2.045 

OSq /2 

2.205 


If we assume that addition of two protons to 
®^Sr does not change the neutron single particle levels, 
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t 


91 

we cen use these values for Zr as well. We c^^ll it 
1 28 

set II. Cohen has also given the single particle 
91 

energies for Zr as 


(MeW) 


ld. 5/2 

0 

2 s 1/2 

1.55 

CM 

2.70 

057/2 

2.70 


We call it set III. YIe mahe calculations with all the 
three sets. This study night enable us to comment on the 
importance of the OgY/p Qi^ergy levels by 

comparing results of set I with those of set II and set III. 
Furthermore, a rel-^tive study of spectra with set II and 
set III may also reveal the relati'^e importance of the 
1 d ^/2 *=^nd 0grj^2 levels because in set III, these two 

single particle levels are degenerate. The r.m.s. charge 

88 

radius for is about 4.1 fm which corresponds 

92 

to b = 2,1284 fm. We take b = 2,2 fm for Zr. The 
experimental energy levels have been taken from ref. 129 , 
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II. 3 RESULTS Uro DISCUSSION: 

Ne first present the results of calculation with 
bare SME. Results for T = 0 levels are presented in 
Figures 2. 3 - 2.6 while those for the T = 1 levels are 
presented inFigures2.7 - 2.14. Results for ^^Zr 
nucleus presented in Figure 2,14 are those calculated 

1 s 

with set II. The single particle energies for "F and 
"1 8 4 " 2 4 2 

'0 and for Sc and '^Ca have been taken to be different 

This will affect the absolute level positions and the 

exact amount of configuration mixing and hence nuclear 

properties calculated with the wave functions so obtained. 

Apart from this, the general qualitative features of 

T = 1 spectra of "'^F and "^^0 and of ^^Sc and are 

expected to be, and in fact they are observed to be, 

6 18 

similar. Hence, we wi3.1 discuss 1=0 spectra of Li, F 
and and T =1 spectra of ^^Ca, ^®Ni and ^^Zr. 

*18 expected, the calculated level energies are 
not in agreement with experiment , maximum discrepancy 
among lowest few levels being for the ground stats. We 
also conclude the oldest knovm result - central force 
gives maximum binding energy. The level energies come 
almost wholly from the central part and, in fact, in 
all the cases studied, the ground state calculated with 
central part alone is in best agreement with experiment 






c+so 



Energy (Me 













Energy (MeV) 
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and, for ground states, the central part gives larger 

binding than the total interaction meaning thereby that 

addition of spin-orbit and tensor parts causes an effective 

repulsion. It can be verified from the figures that, 

1 8 

except for the P, T = 0 levels, even the lowest state 
turns out to be unbound when calculated with the tensor 
part alone. 

Besides level energies, the level splittings 

are also not correctly reproduced, the agreement being 

worse for T = 0 levels, let us, for the moment, concent- 

2 

rate only on the levels arising from the lowest ,j con- 
figuration and consider only the total interaction and 

6 . 18 

its central part. Whereas for Li and 'P nuclei the 
calculated spectra are compressed as compared to experi- 
mental ones, not only is the ground state of Sc 
not correctly reproduced but it is also predicted to 
lie much higher in energy. However, .in view of the 
separation between the bunches of (l”^, 7"^) and (3 > 5 ) 
experimental levels, the same fact might also be inter- 
preted as predicting that the interaction is too weak to 
lower down the level to its right ordered position. 
Correct ground state of ^^Sc is not achieved with 
the total interaction, its components or their combinations. 
The spin-orbit and the tensor parts predict very little 
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£?plittin^ among various J-levels arising from a given 
tv;o-particle configuration and the calculated levels 
with these parts are closely bunched near the single 
particle excitation energy O''"’ 'the two-particle confi- 
guration in question. However, the tensor part gives 


I’^rger splitting than fhe spin-orbit part. This is in 
confirmity with the earlier: observation ’ that the 
two-bodv spin-orbit force is mainlv exhausted in provid- 
ing single-particle spin-orbit splittings. Prom Figure 
2,5 one gets an impression that the central part of the 
T = 0 interaction causes smeller level splitting than 
does the tensor part and this point is more substantiated 
when one makes a similar comparison for T = 1 interaction 


(Figure 2.10) where the observation is not found to be 
true. However, it can only be called accidental as can 
be obseryed from a study of Figure 2.15 where we have 

plotted the puoe (Odj/g)^ 

seperatelT. Here we regain the felth ttot the oentrel 
part causes much larger splittlm than the spln-orblt 
or the tensor part . The obserTed bunching Is fortuitous 
and is caused due to the fact that when dlagonalitatlon^ 
is made, the f level dos® change but the f and 3 
levels are pulled down. the lowest 3 " configuration 

levels, the splittings due to the spin-orbit and tensor 
parts are of nearly equal magnitude but the tensor part 
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b3Comes more important for higher energy configurations. 

It is observed that the spin-orbit and the 
tensor parts do not always split a two-particle confi-. 
guration into different J's in the same direction as 
is done by the central part. Thus, addition of these 
parts to the central part often results in change in the . 
ordering the levels. The addition of tensor to the 
central part results in a drastic increase in the splitt- 
ing of the 2| and I 2 levels of ^Li. The shuffling of 

levels as a result of addition of tensor interaction 

42 

can be seen from spectra of Sc. In this case, although 
the central and tensor parts predict the same level 
ordering, the ordering of and 5'| levels with central 
+ tensor interactions is reversed. The 3-| and 5 
levels of are predicted to be reversely ordered with 
the central part than with the total interaction. 

However, for lowest few levels, the said shuffling 
occurs for T “ 0 levels only thereby -demonstrating the 
relative importance of the tensor force for T = 0 states. 
This study is similar to one made by Young who has 
carried out shell model calculations with purely 
phenomenological central and central + tensor inter- 
actions. 

% 

Although the ground state <1 as calculated with 
tensor part does not agree for T = 0 levels, also (except 
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18 + 

for F), the lowest O"^, T = 1 level, in all cases, 
is predicted at a positive energy value well above 
other levels arising from the lower 1 conf iguration. 

This explains why is it that the difference between the 
energies of the lowest - state , calculated with the total 
interaction and with the central part respectively, for 
T = 1 states is much larger than the corresponding 
difference for T = 0 states. The tensor part appears 

1 s 

to be playing an important role for P nucleus. Except 

for the elevation of the lowest O"*” level, other J levs, s 

of T = 1 from lowest j configuration, calculated with 

the tensor part, are more closely bunched than the 

zero isospin levels from the same configuration. This 

arises because of the fact that we have included the 

3 3 

strong off-diagonal tensor force (^S.| - I).|) that 

operates in T=0 states only. Otherwise, the tensor 
conpononts of the pure ^^.,^2,3 channels of T = 0 and 

q n channels of T = 1 are much weaker, ks calculatedj 
the lowest state in turns out to be T = 1 with 
central part; 3 '*', T = 0 with central + spin - orbit 
parts and 1"^ » T — 0 with central tensor parts. ThUo, 
presence of tensor interaction is necessary to achieve 
the correct ground state J for this nucleus. The 
numbers are given in Table 2.2. 
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TABLE 2.2 


Ener'fties 9 Tid of lowest levels of respective 

1 ft 

isospins in F, as calculated with various inter- 


actions. 


T = 0 


Interaction 


Energy 

(MeV) 


Energy 

(MeV) 


Total 


- 2.76901 


•2.64716 


Central i 1 


- 2.90082 


■3.04721 


Central + 5 + 

spin-orbit 


-2.93859 


-2.89633 


Central + 
tensor 


-2.76922 


-2.76369 
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Overall, it can. be said that the agreement for 
T = 1 levels is better than for 1=0 ones. 

4. study of eigenvectors obtained with various 

parts of the force reveals that the central part alone 

causes appreciable mixing of configurations. Often, 

t^’ere is shuffling of energy levels from ordering tha,t 

would otherwise be expected from a knowledge of diagonal 

18 + ' 

matrix elements alone. For example, in F, the 

^2 

state lies at - 1.32 Me? and comes from (1s^y'2^ 
configuration. The 1^ state lies at - 0.57 Me? and 
comes from (0d^^2) configuration. These are the 
values of diagonal matrix elements. However, after 
diagonalization, the 1^ st.ate shifts to - 2.9 Me? 
energy and is predominantly a (Od^y' 2 ) configuration 
(48 per cent contribution) with large admixture of 
(lai/ 2 )^ conf iguration (30 per cent contribution). 

The it state occurring at - O .71 Me? energy receives 

^ 2 • • • 

68 per cent contribution from (ls^y' 2 ^ cohf igur.ation and 

29 per cent contribution from (Od^y' 2 ) configuration. 

The matrix elements of the spin-orbit part are, in 

general, quite small and the configurations are nearly 

pur^?. The tensor part of the interaction alone causes 

appreciable mixing but less than that due to the central 

part. 
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Let us discuss the role of different sets of 

92 

single particle energies in the fr nucleus. The 
results are shown in Figure 2.16. We have plotted only 
the spectra calculated with the total interaction and 
with the central part. However, it should be clearly 
bo^-nie that not too serious comparison of the level 
positions among the three sets should be made since they 
have been calculated with d ifferent sets of single 
particle energies. However, one thing is clearly seen 
that the set II gives the best agreement with experimental 
ground state. The fact that 6 ^ level calculated with 
set II falls lower than the one calculated wath set I 
can be attributed to va^^ious reasons. Let the single 
particle energies have any va,lue , but inclusion of a new 
configuration in the configuration space will invariably 
bring the olj’ level down. Even if we had done calculations 
with set II with Og^/g configuration chopped off, the 
calculated ground state would still be expected to bo 
lo'v/er than ths? one with sot I since the single particle 
levels of set II are more closely spaced, and hence 
more configuration mixing, than those of set I. Similarly* 
the fact that 0 ^ level with set III lies in between 
those with set I and set II can be explained on the 
basis that single particle levels of set III are more 
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widely spaced thsn those of set I and set II. The 
relative importance of the siiigle particle levels 1 d^y '2 
and Ogrj ^2 onn be asserted from a study of the wave 

functions obtained with set II and set III respectively. 

92 

In the ■ Zr nucleus, within the model space considered, 
we expect only four 0'*' levels. It turns out tha'f' both v/ith 
set 11 and set III, when calculated with the central 
part, the 0^ level receives largest contribution from 
( 0 gyy 2 ^ configuration while the (ld^^ 2 ) configuration 
constitutes most of the 0^ level. This should be viewed 
together with the fact that in set II the 0gy^2 
lies above the 1d^^2 The actual numbers are given in 

Table 2.3. This observation can be attributed to the 
pairing force component of the SME which acts more 
strongly in states of higher angular momentum. This 
effect is observed onl^’^ witn the central part of the 
force and not with the spin-orbit and tensor perts. How- 
ever, it is acceptable in view of the fact that spin-orbit 
and tensor forces^^hre comparatively longer range forces 
whereas pairing force is very short ranged. Sinailt^r situ- 
ation also occurs for total interaction, however, only 
for set III. First observation means that if only 
central forces are operative, neutrons prefer to fill 
the OS'j/2 first -because of the typical pairing 

force. However, inclusion of spin-orbit and tensor parts 
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2.3 

Rel^ti/e intensities' of confi- 

gurations in the structure of oj and ot levels of '^^Zr 
as c'='lculated with the total interaction and its central 
part respectively. 


Interaction 

Level 

Set 

(115/2)^ 
(% ) 

II 

(OSy ^2 ^ 

(3: ) 

Set 

ddj/gjZ 

(“4 ) 

III 

(Og’y ^2 ) 

(X ) 



95 

0.02 

1 

97 

Total 

\ 

i 

j 

0.01 

99 

95 

2 

j 

i 


19 

77 

8 

89 

Cental | 

i 

i 

i 

1 


77 

21 

i 

89 

9 
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weakens the doiPinating e.ffect the pairing component 
to the extent that in presence of all the three compo- 
nents , the neutrons prefer to fill orbits in order of 
oingle p'=irticle energies. However, if two orbits of 
different angular momenta but s^me single particle 
onorgv ore encountered, the neutrons still prefer to 
' ill the one of iiigher angular momentum first. This, 
thus, demonstrates the net dominance of tho central part 
over others. 


A studv of Table 2.3 reveals that calculat- 
ions with central part for set II yield more configur- 
ation mixing between the 0^ and 0^ levels ever^though the 
1 dj /2 and particle levels are degenerrote for 

sot III. It is to be attributed to different values of 

matrix elements of the potential in the Cld^/ 2 ) ^ and 
2 0 
( 0 ^ 7 / 2 ^0"^ configurations. The exact numbers for the 

(ldj/ 2 )^_^ and (0grj^2)^^ conf igurat ions are - 0,17415 Me? 

and -0.17861 MeY respectively for the total interaction 
and -0.55563 MeY and - 0.99038 MeY respectively for the 
centpral part. The unperturbed matrix elements, including 
single particle energies, for set II, and set III are 
shown in Table 2.4. The corresponding perturbed levels 


are also shown. For central part, the unperturbed 

P 

^ level lies lower than the unperturbed (ld^y' 2 )' 



Perturbed (E) ■end unperturbed (E^) energies of 0^ and 

j\- 0 2 

0^ Isvels of Zr celculeted with totel interaction 
•and its central part. 






lUnfiGrtal^ed Energies 

: (^n) 

®o^'’‘^5/2^ 

T" 

j Perturbed Ener- 
I gies (E) 

2(0^) 

iter- 

ition 


i 



. P 

(MeV) 

0^ ^ 
^4 

(MeT) 

-E(0+) 

(MeY) 


i 

i 

(B'ley) 

(MeV) 

(MeV ) 


Set 

1 

II 

3.91585 

4.25139 

-0.51554 

4.19182 

4.29737 

^ 0.10555 

• t 

Sot 

Illj 

1 

5.22585 

5.22159 

0.00446 

5.27525 

5.45133 

0.17610 

1 

Set 

II 1 

3.55457 

3.41962 

0.15475 

i 

5.58766 

3.81607 

0.42841 

jntral 

Set 

iir 

4.86437 

i 

i 

1 

t 

4,4096 2 

i 

i 

I 

0.45475 

i 

' i 

4.59721 

5.06341 

0.66620 i 
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level for both set II and set III. As a result of 
configuration interaction, the 0^ - 0^ splitting 
further spreads up. Samething also happens for set III. 

We may be using anv set of single particle energies but 
the off-diagonal matrix elements, which cause the 
configurations to mix, always stay the same. The fact 
that for central part, set II is a case of stronger 
configuration mixing can now be explained on the basis 
that the two unperturbed 0*^ levels are much closer for 
set II than for set III. However, calculations with total 
interaction employing set III yield less mixing than even 
the calculations with central part and set III eventhough 
the splitting between the unperturbed O"*" levels in this 
case is very small. This may be attributed to the fact 
that the off-diagonal matrix elements for the total 
interaction will be di-^ferent from the ones for the 
central part. The difference between the matrix elements ; 
of total interaction in (1d^^2^Q+ O’*" 

configurations is so sm-^.l! that the largeness of single 
particle energy of the OS’j /2 dominates ^ 

over this difference and the unperturbed (0g,^^2^ 0"^ ! 

level itself lies above the (11^/2^ 0"^ level. Config- | 

uration interaction sends them further apart and it turns oui 

i 

to be a case of almost pure configurations. This study | 

1 

-also proves that the closeness of two unperturbed levels | 



is not sufficient to ensure large conf igurc^ion mixing 
in wsve function. Off-diegons.1 m<=itrix elements else 
eouelly metter. However, for « given set of of f-diegon-^1 
Hieti-ix elements, the closer the unperturbed energies, 
the stronger is the conf igur--^tion mixing. 

Energy l.-jvels o^ ^li calculated v/ith various inter- 
actions have been plotted in figures 2. d, 2.4, 2.7 and 2,8 
-as a function of b. It is seen that the snin-orbit 
part stays almost constant, largest vari^^tions occur in 
the central part. The tensor part also varies appreciably 
but in between the central and the spin-orbit variations, 
however, only for T = 0 interaction, for T = 1 interaction, 
the variation with b, of tensor part, is almost as much 
as that of the spin-*rbit part. >lt any b value, addition 
of tensor to central part causes, a large depression in 
2 ^, 1 = 0 level and the 2 ^ level c'^lculated with central 
+ tensor interactions is also the one which varies most 
sharply with b vo.lue. 

YsTe now come to the renormalization part. 7/e 
have seen that the central fore'e constitutes the ma;ior 
and also the most important component of the two-body 
force. Therefore, it is anticipated that changes in 
the matrix elements of this part might bring in the desired 
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result. In qny cese, the spectra with pure spin-orbit 
tensor forces ai^e nowhere near the truth and it will 
be a futile o,ttenipt to change these parts only. Mathemati- 
cally, it might to possiule to change Tn'^trix elements of 
these foreSvS only an^i got a new set of matrix elements 
which could give the right 'numbers' but physically it 
.a^ppoars difficult to conceive th'^t these two components in 
bare interaction have so less importance anrl ths renormaliz- 
ation should change these parts alone. On the other hand, 
we have other more strong reasons to believe th 3 ,t it might 
be sufficient to change the central part alone. In any case, 
it is safe to neglect the spin-orbit part altogether. 

Spectra with central + spin-orbit and central+tensor inter- 
actions respectively are in general resemblance with the 
spectra with the pure control part both, as regards level 
energies as well as level orderings. However, .for finer 
details, we have seen the importance of the tensor force, 
p.articularly for I == 0 interaction, llso, foe example, che 
tensor .force, in second-order, has bean loiown to contri- 
bate about 150 Me? to 0 binding energy. 

Let us first discuss a possib3.e way of finding a 
suitable renormalizing force for 1=0 SME. Let it be a 
T=0ora!E=1 state, particularly for ground state, it 
is expected to be predominantly a relative S-state. Eor 
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T = 0 it becomes state ant for T = 1 it is Sq 
state. Even experimentallv , the lareest inaccuracy 
inlieterrainin'^ the matrix elemen'"s fr::!! scattering data 

3 

occurs in the state. The coupling parameter 
3 3 

E between the ' S.j and the channels is also poorly 

3 3 3 

known. This fact makes both and S^- -- T)^ matrix 

elop'ents ' values uncertain. However, the tensor force 

operates in spin triplet states only and, thus, effectively 

renormalizes the triplet S - and D - state central 

interaction. The relative matrix elements decrease 

rapidly with . Neglecting odd state interactions for 

the moment (which are expected to be much weaker), one 

might then conclude that the tensor force also effectively 

- 

renormalizes the S.| state interaction only and that 

3 

the central part is made almost wholly of the S.| state 

3 

interaction. It is then expected th'^t if only matrix 

elements are arbitrarily varied to get a fit with experi- 
ment, the net change in these matrix elements will 
imbibe in itself several effects, namelv, the inaccuracy 
in the bare interaction matrix elements themselves, 

the renormalization of ^S.j part itself and the effective 
renormalization of interaction as a result of 

possible renormalization of the tensor force in the 
coupled channel. 


I 
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Therefore, our procedure has been th'=t we make 
th'=‘ S^- relatiwe state matrix elements (1 + a) times 

their bare interaction value and calculate the spectra 
again with the total interaction so obtained, each time 
increasing the value of a and thus pinpoint a value of 
a such tii'-.t a particular level in question is in complete 
agreement with experiment. There is no a priori reason 
to believe that a single value of a will reproduce all 
the levels. However, if the values of a as obtained for 
different levels are not widely different, a least square 
fit can be attempted to obtain a compromise value of 
a spch that the considered levels are best, although not 
exactly, reproduced. Y/e have minimised the function 



Z 

i 



, calc 



exp 


) 


2 


where the siimmation over i runs over the levels consi- 
dered. The T-- ,-1 's are the level energies calculated 
as a function o-p a and E. ^^„'s are the correeponding 

-i- y 0Xp 

experimental energies. 

We now discuss the T = 1 interaction. The 

level of Li as calculated with total interaction^ for 

any b value, is not in agreement with experiment. Similar 

18 4-2 58 . 92 

situation occurs for Ca, Hi and ' Zr nuclei. 
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However 5 there is a difference. In these latter nuclei, 
the lowest O'*" level occurs as ground state whereas for 
°Li it is an excited state. The calculated 0^ level in 
these nuclei is alwavs above the exDe-r'imeatal ground 
state indicating an insufficient attraction . On the other 
hand, for ^Li , the calculated 0^ level turns out to be 
below the expe^'ninental one meaning that a repulsive 
re norrpalization is needed. Is discussed, the Sq part 
0 ^ the interaction plays a dominant role for lov^-lving 
T = 1 sta-^es. In a state where the spins are paired and 

the relative orbital angular momentum o^he particles is 

! 

zero, the force keeping them in this state is anticipated 
to be a relatively short ranged one. The exact potential 
that the SME represent is unknown but with a hard core 
interaction, the major part of attraction is cut off since 
it is absorbed in compensating for the hard core. Thus 
effectively, a hard core realistic interaction is shorter 
in range than an equivalent well beliercd smooth potential, 
like a phenoraenologic^^l potential with Gaussian shape. 

The observation regarding the 0^ level suggests that 
the realistic interaction, let it he sm or others which 
have been worked upon n number of times, leeks sufficient 
pairing. If this interpretation is right, it can be 
checked by adding a pairing force to the bare sm and 
seeing if it is possible to get correct ground state 
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energy at least. 4s a primary calculation we have 
used a constant strength pairing force defined bi?" 

^ j-] 2 2 • ^3^4 ^ 

= - - hi f'hl'TTK'iJ'^T) i h- i ^2.9) 

2 ' ^ ^ ^ d -j 2_ ^ ^ 

The results ^re presented in Table 2.5. These calculations 
have been done for ”^^0’ ^®Ni and ^^Zr nuclei 

only. We see it ^ possible to get the correct ground 

T'lBLE 2.5 

Pairing force strenffths |Grj for respective nuclei such 
+ 

that the 0.j level energy is correctly" reproduced. 


Nucleus 


hi 

(MeV) 

18o 


0.254i 



0.2480 

. 58^1 


0.0864 



0.0680 
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st^t.e energy by =5 suitj^ble choice O’^’ the strength [G-j . 

The vslence nucleons in ^Li are e neutron end a proton 

end an idealized pairing force as the one defined by 

eq* (2,0). cannot be used for neutron-proton pairing. Hor 

the s»ne reason, these calculations have not been tried 

1 R AP 

■^or T = 1 states of ? and Sc. For the particular 

6 + 

case of Li, the 0^ level demands a repulsive renorma- 
lization rather than pairing force. In these calculat ions , 
the first excited 0 state cannot be reproduced with the 
same pairing force strength which reproduces the 
ground state energy. For excited 0'*’ levels the calculated 
ene^^gy level varies very slowly -as a function of the 
pairing force strength and if at all we want to reproduce 
this level too, mathematically it might be possible to get 
i’he right number but a physically unacceptable large 
value of jG-j would be required. This was perceived from 
the range of values of jG| during the calculation for 
obtaining the correct ground sta'^e energy and hence was 
not tried. However, these results are not unexpected 
cither, in view of the fact that the pairing theory is 
panticularly successful for ground states.' 

Let us again concentrate only on the levels 

2 

arising from the lowest .i configuration. • Making a 

IX 

comparison with corresponding energy levels of same J , 
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ono observes th^t the theory predicts a much compressed 
wSpectrum as comp=ired to experiment. Thus, the agreement of 
energies of' J O"*" levels is better than that of the 
0^ level. Hovvever, it invariably turns out that as a 
result of compression, the levels cf highest J compatible 

p 

with the said 1 configuration (and wl: ich also happens 
to lie hlvhest in energy) always falls below the corres- 
ponding experimental level. This means that whereas 
one Invariably needs an attractive reno-arnali'sing force 
to get correct odj” level, one needs an attracti\7'3 or 
repulsive renormalization for higher levels, but, as 

discussed, always a repulsive renormalization for 

42 

highest J at least, for example, in Ca, on comparing 

th*'> results for total interaction with experimental 

ones, the O"^ level is predicted to be too high but the 

discrepancy is much smaller for 2'| level, rhe 4-^ level 

is in excellent agreement v\rith experiment and the 6^ 

level is again in disagreement but below? the experimental 

58 

level. Similarly, the ground state configur=it ion of Hi 
is (1p^/2^^* level is above while 2 ^ level is 

below the corresponding experimental levels. 

On.0 xmr^ 6 d i ^ "fc 6 iro.pxov'oinBn.'t oyst "th-O pnix*in,S 
force model for the renormalizing force, that has been 
shown above to be successful for reproduction of correct 
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iSrroand state energy, is the replacement of the idealized 
pairing force by a 6-function force 

V = Vq 6 (r^ “ ^2 ^ (2, in) 

The justification ■'^or the use of a 6-nu;iction force is 
that for the lov;est O"*" level it beh-^yes exactly .as 
does th-> pairing force. However, it also has nonvanishing. 
m-^trix olvCneits in J O"^ states and it is, therefore, 
ant ic ip-‘t ed th^t suitable adjustments of strength of the 
6-function force, v;hen added to the bare SIE, might 
compensate for th.‘ renormalizing effects. However, it 
is well established that the 6-function force gives 
very little splitting in J ^ levels while lowering 
the ground state O"'' level and, ther^re , it can be seen 
beforehand that the strength of 6-f unction force that 
givv’s correct ground state energy will not necessarily 
reproduce the higher levels also. Matrix elements of 
til''! 6-function potential are calculated using formula 
(i.28) given in appendix i. The results a,.e presented 
in Table 2.6.''inlv results for the lowest j^ configuration 
levels are presented and, in fact, these were also the 
only levels for which some sensible value of could 
be found. Otherv/ise, other excited sta.tes were found 
to vary very slowly with VQ-aid for large they almost 
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TfiBLni; 2.6 


I 

II 

III 

TV 

Y 

YI 

Nucleus 

J 

Bxp. 

(MeV) 

B=>ra 
(Me 7) 

(MeV-fm^) 

(MeY) 

18o 

0 

-3.9044 

-2.5865 

-95.35 

-3.90014 


2 

-1.9223 

-1 .6254 

-312.50 

-1.69464 


4 

-0,3515 

-0.6887 


-0.70959 


0 

-0.2727 

-0.1402 

- 25.55 

- 0 .65336 

^4,-. 

0 

-3.109 

-1.3635 

-190.95 

- 3.10105 


2 

-1.589 

-0.8172 

-504.00 

-1.05356 


4 

-0.359 

-0.3879 


- 0.52468 


6- 

0.081 

-0.2437 


-0.28357 


0 

-1.936 

-1.5723 

-193.80 

-1.9383 


2 

-0.482 

-0.7867 


-0.79042 


92 ^ 

-I. 434 OO -1.0278 

-73.875 

-1.43368 

2 

-0.49954 - 0.441 

- 85 . 866 

-0.49117 

4 

0.061600 - 0.2452 


-0.25267 
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1 , 

becnme ^sympljofcic to some value. As anticipated, the 

0^ level is reproducible for all the nuclei. Hovrever, 

not all the levels, even, those arising from the lo'/vest 

.2 

0 CO nf igur-ucion , are reproducible even with different 

for di'farent levels. In all the four nuclei 
considered, only those levels, other than the 0^ level, 
vere reproducible for which calculated energy was higher 
than exoerirDGntal one, and, except for ''>li, only the 2^ 
levc'l was ■^ound to &.11 in this category, for others it 
could have been possible -to get the right number by 
'assuming a positive V^. However, that would be un- 
physical. .Por this reason, the agreement of the said 

level, when calculated for that gives correct 0^ 
■■.•‘lergv, Improves (■■■•Ithough exact congruence with 2^ 

I "cperimcntal level energy is not achieved, since that 
r'KiUiirec a di-^fernnt from the one required for the 
Icv.-'l). The agreement for other levels becomes worse, 
since, in principle, they require a positive . 
.Incidentally , the Og level of ”*^0 occuring at -0.2727 
it-'V energy is obtainable for a much weaher than the 
ons roauired to get correct O’! level energy. However, 
the exact shell model configuration of experimental 0 ^ 
level is uncertain and it has been felt that it contains 
4p - 2h de.formed component which should be explicitly 
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included in the model space rather than being 
treated in perturbation theory. Thus, although we ha/e 
been able to get the right number, yet, for this 
reason, we reserve that our v/ave function for this 
level may not represent its real structure. 

We, thus, sec that, to some extent, it is possible 
to achieve an improved description of level energies by 
trving to simulate renormalisation corrections via a 
6-function force. If we stop here, we endup with a 
J-dependent effective interaction-thr.t too for very 
few selected levels. The idea of a repulsive 6-function 
force for bound states seems to be more a product of 
fantasy rather than of ,any physical meaning. The addition 
of 6-function force only ranormali'ses the short range 
P'=rt of the T = 1 SME and we would like to attribute 
the observc.tion, that addition of 6-function force only 
is not sufficient for higher levels which need an 
effective repulsion, to the fact that the long range 
part of T = 1 S?® also needs be prop'jrly renormalized. 
There is conceptual difficulty in assuming a repulsive 
6-function force but while constructing a model for long 
range renormalization we may remember the caution of 
not introducing a function where such difficulty arises 
again. 
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When we calculate the two-body m'^trix elements 
(TBICS) of 6-function potential empolyin.'^ the method of 
Slater integrals, although all allowed values of k in 
the expansion 




E 

k 




f. 


k 


( 2.11 ) 


contribute but, finally, owing to the special features 
of the 6-function force the geometry of expressions 

occuring in such a way, the TBl/Di) is proportional to 
only i.e. effectively only k = 0 component contributes. 
On the other hand, for a non-zero ranged force, all k 
values contribute differently, a compact expression invol 
ving only one is not derivable and the relative 
contributions of different k's depend upo-n the range of 
force. Roughlv, the highest k value making appreciable 
contribution varies inversely as the range of the 
force, let us, for a moment, make multipole expansion 
of a finite range Gaussian force 


.-rPa^ 


= 1 - 


r^ 

T 

a 


+ 


2a 


( 2 . 12 ) 


let us work in the limit a In principle , r-the 

two-nucleon separation - also goes to infinity but 
let us work in a model where wave function fall off for 
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l^rge internuclGon separation 

r = l - ?2 I ■ 

Defining . 

' R = -I ( ?, + ?2 ) 

as the centre-of-mass coor'llnate , 

= 2(r^ + rg ) - 4 (2.13) 

TJow, on right of eq. (2.12), first term is 
a constant and ca,uses no splitting. The second term 
can be expanded using eq. (2.13). On right of eq. (2.13), 
first term represents an average oscillator field and 
the second term represents the motion of the centre-of- 
mass as a whole. Hence, the second term O'f' eq. (2.12) 
also does not give any level splitting. The third term 
can be expa,nded as 



where 9 is the angle between r^ and r 2 > 

The’ first two terms raodifv the oscillator field . In 
effect, they introduce anharmonicity. The last term 
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has odd p^^rity and does not turn out to be Talld 
representation of ]\T-¥ interaction in physical problems 
where parity of a two-nucleon st^tte is conserved. 

Within a major shell, th’-^ third term gives the same 
contribution in all two-particle states. It is the 
fourth term th-^t h'-^s a P 2 (Cos 9) type angular dependence, 
conserving parity and which might be expected to produce 
level splittings. 

Therefore, in the long range limit of a, that 
the expansion on right side of eq. (2.12) could be 
terminated at the third term, we see that a long range 
d-il force, apart from providing the averaas harmonic 
oscillator field in which the nucleons m.ove, can be 
approximately simulated by a dependence of type 

^1 ^2 ^2 0 ) 

We have derived this particul'^r form for Gaussian 
shape. However, we expect that, for anv general shape 
of the potential, the angular dependence can be given 
by P 2 (Cos 0 ). This discussion also explains the fact- 
the average nuclear field is provided mostly bv the 
long range pf\rt of the interaction. Were it not 
for the largeness of a, it would not have boen possible 
to terminate the series (2.12) so soon and thsn the 
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depsndance on and r2 O'f angle-independent terras 

! ■ 

would be much more complicated. 

Thus, in go-ieral, we anticipate th'^.t a long 
range force can be appr.G'jcimate'ly represent -d by 

f (p^, ) P2 (Cos e) 

where f (r^ , r2) is some radial dependence. However, it 
is customary to take 

f (r^, rj) = 

In this approximation, calculations become extremely 
simple. The justification for such a simple choice is 
that if the potential 

V P2 (Cos 0) 

is included as a part of two— nucleon potential, it 
ganeratc.s the deformed field needed to explain quadrupole 
deformation of deformed nuclei. 

The pairing-plus-quadrupole approach has been 
used”'^^ extensively in nuclear structure calculations. 
In these calculations, this sum is taken « s the total 
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rcsidupil interaction. However, in our calculations, 
we are only taking it -as the renormalizing force, 
the pairing -iud the quadrupole p'^rts renormalizing the 
short '^nd tho long range p^^rts of the free force 
respectively, for the pairing part we haae used 
6-function force instead of the idealized constant 
strength pairlnv force. The choice of pairing and 
quadrupole forces has been for convenience and also due 
to our knowledge that these forces have earlier been 
successfully used. Otherwise, if the philosophy is 
only to add a short and a long range part whose strength 
parameters are to be determined from a fit to experimaa'" ' , 
data, in principle, any phenomenological form for thesc- 
p.arts should work. Also, such an approach would be 
expected to be more successful it the pairing-plus-quadrap'; 
approach is not very successful for many excited levelfn 
This is anticipated bec=iuse the short range part will 
largely renormalize the higher k components of the free 
interaction while tho long range part will mainlv contribui 
to low k components. Th'^t way, a combin'''tion of them 
would give a complete multipol-- expansion of the renor- 
malizing force. However, it h-s not been tried sinco : 

calculations are anticipated to be more involved and if i 

I 

several strength and range parameters are included, | 

that would spoil our case of formulating an effective | 
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inter-iction in whose derivation o minimum number of 
experiraentnl levels ere used. For this reason we hnve 
restricteci upto pairing nnd qu- d-^npole forces only. 

Although wo h^"*' given o different explonotion 
for use of qusdrupol' force os e component of renormaliz- 
ing force, v;e could hnv2 anticipated at least for the 
angular dependence to be of the type P2 (Cos 9) from 
an analysis of second-order perturbation theory 
results as done earlier in this chapter. 

Matrix elements of the quadrupole force are 
calcula.ted using aq. (A. 24) as given in appendix A. 
However, it can be seen that inclusion of this new 
renormalizing force increases the number of unknown 
paq-’ameters to three (the 6-function force strength and 
the quadrupole force strengths in spin-singlet and triplet 
states respectively) from one and thus the basic point, 
that the theory should contain no or minimum number 
of parameters, behind the use of realistic interactions 
se'-ms to bo losing, ’since , like this, we can go on 
adding physically accxaptable' renormalizing forces one 
after the another • and then one can question the use of 
such a method at all - compared to such a procedure, 
the use of a completely phenomenological force will 
also be equally .justified. To avoid this objection, we 
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adopt a procedure wherein we effectively use only one 
experimental level for determination of tho three 
par-<ni 3 ters, although the computer time is incrc'^sod in 
this process. We use only th'; ground state energy 
fo"’^ vdeterminn tion co the p‘^r‘^''^eters . Pirst a rougii 
guess of the upper limit of parameters is obt^i-ned in 
such way that, other two parameters remaining zero, 
the two-body matrix eleme it for 0^ level itself becomes 
equal to the experimental level energy. This is tte 
upper limit value for this parameter since even at this 
value, as the configuration mixing calculation is done, 
the calculated ground state will certainly move belov; 
its unperturb 2d value (i.e. the experimental value) and 
therefore, if, starting from zero, the strength for this 
particular parameter is increased in small steps, even 
if the other two parameters do not contribute at all, 
we will be able to pin-point a value where the calcul- 
ated ground state energy, v/ith this parmicter only 
nonvanishing, ag^^eos with the experimental ground state 

energy. 

In any case, the 6- fu action force has to be kept 
attractive but the quadrupole force strength- can be taken 
to be attractive or repulsive. lu fact, when all the 
three -parameters are simultaneously varied, there can be 
a number of sets of values that give the same description 
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of the ground state energy. Now, it is here that, if 
required, we c=^at.^te one or both spin components of 
the quadrupole force tc be rep’ilsxT-’. The quadrupole 
force does contribute in O'*’ states also. Hence, 
introduction of a repulsive component in the effective 
force requires the attractive part to be strongf^r in 
order to get correct ground state energv. How, the 
6-function force contributes to J 7^ o'*' states also. 
Therefore, eventually, a limit to the variation of 
parameters is set, both on positive and negative side 
beyond which agreement of higher levels cannnot further 
improve , 


After pin-pointing parameter values for ground 
state energy, excited states are calculated ^fith each 
set. It then turns out that not all the sets can 
describe the excited states equally well- It would be 
desirable that a single set of parameters describes all 
the excited states. However, if it is not possible, 
we can, at le'^st, find a set, different for different 
levels, that gives best agreement for one level at a 
time. Thus, effectively we will end up with a J-dependent 
effective interaction. 

Since the T = 1 levels of ^Li require a repulsive 
renormalization, we cannot use a 6-function force. 
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experimental levels as compared with results for b = 1.5fni 
and 1.6 fm but the 1'^ level is in worst agreement which 
me'^ns that a comparatively much larger renormalization 
is needed for the level. The results are tabulated 
in Table 2.7 ani plotted in Figure 2.17. 

One thing can be observed - except for b = 1-7fiii, 
for b = 1.5 fm and 1.6 fm, the value of a for the 2 '^ 
level is almost twice that for the 3"^ level. It may also 
be noted that the 2 '^ level comes from the (0p^^2 '^ 13 / 2 ^ 
configuration while the 3*^ level comes from ttie ground 
state configuration (Op^y' 2 ) • have earlier made 

a mention that the low-lying states are predominantly 
relative S -states. This observation makes us suspect 
that may be for the 2"^ level, effects of higher 
partial waves cannot be entirely simulated in the 
matrix elements and that may be this approach for renormali- 
zation is more successful only for levels arising from 
the lowest configuration. Even otherwise, we know 
that a good fitting is possible only if the 

numbers to be fitted are not very much different from one 
another. Therefore, we have repeated the calculations for 
36 ^ fitting without including the 2”^ level. Largest 

O 

change in % from this omission occurs for b = 1.5 fm 
where the best fit value of a becomes positive after 



1 58 


-h OJ 


+ 




4“ 

OJ 


o 

s 




OJ 

m 

tr\ 

* 

O' 

I 


CD 

\o 

fON 

♦ 

f 


OJ 

CD 

cn 


C\J 

LO 

IPs 

•• 

vx) 

I 


• {:> 
{x! S 


!>- 

!>- 

UO 


VO 

CD 

CO 

4 

uo 


a) > 

u o 

C: rS 


CO 

LO 

a\ 

o 

00 

cr 

o 

« 

o 


o 


CNJ 




on 

on 

to- 


co 


Lf\ 

I 


-P 

•H 

Cf-U 

4.D 

W I 
0" 
CD 


O 

d 

x:! 

-P 

•P 


OJ 

> 

OJ ^ 0 


ur\ 


d 


VD 


m 


ir\ 

P- 

OJ 

\D 



cc 


KT 



0- 


CO 


00 



LA 


CO 


o 



cr 


ur\ 


T— 



VD 

v- 

r — 




cr 

OJ 

• 

» 

T-” 

# 

!>■ 

cr- 

• 

Lr\ 

o 

CD 


0- 

cr 

LA 

i 

[ 

cr 

I 

O 

— 

i 



cr 


VO 

CJ 




p- 


o 

o 




O 


• 

4 




4 


tA 

o 




o 






GO 

I 




!>■ 


P" 


00 



LA 

if\ 



LA 







O'' 



A 

-p- 

CJ 





LA 

• 

# 





♦ 

tA 

O 


CA 



A 


p- 

r- 

VD 

cr. 

Lr\ 

I 


-p 

•H 

Ch^ 

-p 

m 

0 

pn 


0- 


m 

VD 

p~ 

LO 


d 

> rd 
0 

‘=^.' ^ 
fs:^, *r-f 


Cvi 


OJ 


>* 

0 




00 


Contd. 



Gontd. Table 


138a 





T— 



CD 



VD 


VO 



C'- 



i:-“ 





T— 


CO 



lA 



VD 


Lr\ 

c 


CO 



O 



• 





• 



c 

Q" 


00 

'T- 

A 


l 

I 

o 

1 

A 

A 

1 

CM 



I>- 


t> 

GO 


A 



VD 


00 

O 


-r- 



o 


A 

T- 


VO 





♦ 

« 


» 



o 



o 


c 


l>- 


T-* 



00 



o 


VO 



A 



VsO 





"T- 


VX) 

o 


K> 






T~ 


A” 



A 


« 

• 


« 



» 


C\! 

i 

o 

\ 


CM 

1 



CM 

1 



^■ 








CM 








X) 








00 


CD 





CVi 

CM 


VO 



00 


LA 

CD 


r- 



T— 



VO 


00 



CM 


tA 

•T” 


A“ 



T“ 


♦ 

* 


« 



* 


LA 

o 


A 



MO 

1 



^ 

w e -p 
oj S -H 
is 


> 

CM , CD 


t> rSd 
-P ^ -P 

CQ fs^Li ’H 

rr- 


CV . 0) 



-4.86631 -1.63 -4.06 1.44235 


138b 


N'M 

CT^ 

— 

MO 

•T — 

-:;h 

CM 

O 


iS^ 

• 

♦ 

T— 

— 


a? 


CM 

t>- 


^_ 

u 


LTN 

r- 


OT' 

o 


r- 

» 


• 

o 

oo 



*r~ 

1 





cc 



v~ 



• 



o 


cr^ 


* CM 

00 



00 


ir\ 

0) 


LTn 

o 


CO 

« 


m 


♦ • 

c 


I 


00 

OD 

00 

• t 


MO 


'M" 

m 



1 

"M- 

MO 

p 


O 


o 


KO 

LTn 

CM 



* 

• 


# 

o 

o 


o 


oo 

LTv 

cJ 

t 


CO > 
0" feEj 


'sO 

VX) 
% — 

in 

CM 

O 


CM 




O 

cc 

C' 

\4D 

\ 


•r"! 

CHf- 

CO 

(1J> 

po 


o 


CM 


CM 


> 

o 

F&., -H 




t>* 

CD 


in 


VO 


00 

CM 

VD 

I 


-P 

•p b 

> 

-p cu -P 
m s -H 

CDw > 

P. 






00 


r>- 




141 


the oraission. This should be viewed together wlhli 
the f'^ct th^t, with 2^ level included in tine fit, 
the 1"^ ground siato renormnlior-jd Ijvol is in even 
worse cgreoinont with experimental level theii wlnet 
ngreeraent the b'^re interaction itself C'in give. Omiss- 
ion of the 2"*" lovci from fit causes the agreement 
to improve on renormalization. S.t all the tlnree b 
V'^lues, improves when the 2 level is not included 

in the fit. drestic improvement occurs for b = 1.5 fm. 
for b = 1.7 fm, the 2^ level calculated with b^'re 
interaction is not very much different from experimental 

level and also the value of a for this level nt th.is b 

2 

value is positive. Hence, improvement . in %. for b = 1.7 
fm is not as much as it is for b = 1,5 fm. 

Jhe I 2 lovel is also explained better efter 
modification 0 ? the matrix elements. However, it 

varies very slowly with oc anfi it could be perceived tliat 
a very large value of a would be required if at all this 
lcv:l docs not as^ptoticallv approach a constant value 
ns was also t>u. Indication during the range of a values 
considered. This h-s bean the criterion througHout that 
calculation of a for all those levels v/as not tried for 
which it could be observed that a physically inconceivably 
large value of a would be reauired. Ve can only reason 
that it is not justified to assume that the renormalization 
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of the force in various channels can be effectively 
dumped into the matrix elements of the relative 

state alone fcr all the states of a nucleus. 

The spectrum of ^Li is best reproduced for b = 1.7 fm. 

Ho?/ever, it does not necessarily mean that it is also 

the true b value for this necleus. At the most we can 

say that, consistent with 7^^ - 6^ = 3,72MeV, 

^Pl/2 ^^5/2 

b = 1.7 fm is the most reasonable value of b for this 
nucleus. We have no^t made calculations for some still 
higher value of b, say b = 1.8 fm, but it can be seen 
from, the figures that if the levels continue to behave 
with b in the same regular fashion (which is expected 
also in view of the fact that the relative matrix elements 
smoothly with b), the agreement of levels calculated 
with .bare interaction will become worse (particularly 
for the ground state) although it might be possible to 
arrive at a more satisfactory best fit value of a. We 
suppose this ambiguity will continue to prevail unless 
the b value and the single particle energy of the OP-j/2 
level is uniquely established. 

One can reasonably questio,n that if three known 
values are different, a one-parameter fit to two of them 
is bound to be better than the fit to all of them. It 
would then weaken our possible conclusion, that a 





reason.'sble renormalization c^n b^established in terms 

/ 

3 

of matrix elements alone, particularly for the 

2 

levels arising from the lowest 3 configuration, which 
appears to be emerging by the observation of a better 
■;’it to 1"^ and 3"^ levels of on exclusion of the 2 ^ 

level. 


This doub+ can be checked by performing similar 

18 4-2 2 

calculations on 1 and Sc nuclei - the (0d^^2) 

2 

(0fY^2) configurations give respectively three and four 

• i R 

levels. The results for P are tabulated in Table 2,8 and 

42 

plotted in Figure 2.18. Similarly, numbers for Sc are 

shown in Table 2.9 end these are plotted in Figure 2.19. 

It is seen that the agreement does improve a lot when 

2 

only levels from the lowest 0 configuration are considered 

in the fitting procedure eventhough it is also possible 

to find out reasonable values of a for one or two other 

levels as well. The value of a for the 4-'*" level of 

+ 4- -h + 

differs much from the a values for 1 , 3 , 2 and 5 

2 + 
levels and therefore, 9/ is very much improved when 4 

2 

level is excluded from the fit. The improves 

further when the level is also excluded from the fit. 

2 

Hov\^ever, this time the change in the value of X is not 


substantial. Similar results are found in case of 


42 


Sc 


nucleus also. 
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+ 10 ^ . 

The observ-'^t ion -when 2 level of is 

excluded from the fit, the fit does not impT’Ove much- 
desorves further ennlysis. This level crises predominantly 
from the (ls^y '2 '^'^ 5 / 2 ^ co ofigur^tion. The single 
perticle energy of this conf iguretio’:! is ohout 0.7 Me7. 
Results improve in ^'^Sc nlso v/hen the 2'*' level, which 
receives l-^rgest .contribution from (IP 3/2 ^^ 7 / 2 ^ configur- 
ation whose sinrle particle energy is 2 M 6 V, is excluded 

42 

from the fit. Howeyer, the improvement in Sc is not 

-IQ ^ 

es much as it is in r on exclusion of 4 level which 

arises from ( 0 d^y 2 ^*^ 5 / 2 ^ conf igur-^^tion whose single 

particle energy is about 5 MeY. The 2^ level of ^Li 
can also be sirailarlv discussed as it arises from 

( 0 pi /2 *^P 3 / 2 ^ configuration with a single particle energy 

of 3.72 MeY. These all results mean that it is sufficient 

to change only the ^S.j matrix elements to get a renorm- 

.2 

alized force for levels arising from the lowest j 
configuration. This approach is also successful, to 
some extent, for levels arising from excited configurations 
with low single particle excitation energy. Further, 
mathematically it might be possible to find a reasonable 
yalue of a for excited configurations with higher 
excitation energies. Howeyer, with such leyels included, 
it is not possible to arrive at a best fit value of a 
which would describe other low-lying leyels as well as 
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c-^.n be bone by excluding them. From Tables 2.8 and 2.9 

we Can see that this difficulty arises because the 'value 

of a for such levels is much di'fferent from the ones for 
2 18 

the lowest j configuration. For F, it can be seen 

tint the value of a for 4-"^ level is much laT’ger than 

those required for other levels. This might mean that 

the relative dominance of the relative states in 

excited configurations is not as much as it is for tho 

ground state configuration. However, then, we cannot 

+ '’ 1-2 

explain how come is it that the 2 level of ' Sc can 
be explained by a much smaller value of a than the others 
required for the ground state configuration levels. 

We now discuss the results of calculation of 
pairlng--plus-quadrupole renormalization of T = 1 Sl'ffi. 

The results are presented in Figures2.2n - 2.23. In ell 
these figures, modified energy levels, as a result of 
pure 6-function force renormalization, are shown for 
even J's only. The odd - J levels for whom the 6- function 
force does not contribute ava not shown. Us explained, 
there were a large number of sets of strength variables 
possible each of w^hich could give the correct ground 
state energy. However, in these figures, results are 
plotted only for those parameter sets, each of which could 
very well describe at least one excited level. The 
numerical results are given in Tables 2.10 - 2.13. 
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4 study of Figure 2.20 rcveols that the overall 
agreemeut of predicted spectrur with experimeutal one for 

^0 is achieved with p'-^'r-'^meter set (-0-^, 0, -0.03053'16) 

at S.ho. 5 of Table 2.10 i.e. spectrum Y. Similarly , in 
4? 

“Ca , spectrum II gives an overall best agreement. The 
case of ^ Hi is the one of worst agreement that with all 
efforts even the 2^ le'^a?! could not be elevated beyond 

-0.72447 MeY (spectrum I) while it great Iv disturbed 
all other levels. In spectra II, III, lY and YI, an 
improved agreement of excited 2"*” experimental levels 
other than the 2(5” one is seen but we can not say if 

i 

physics— wise also we one correct since we ^re not ma.king 
a one-to-one correspondence between the calculated and 
experimental 2^ levels. Instead, the said agreement is 
established if the 2^ calculated level is ignored from 
the comparison and the 2^ experimental level is compared 
wi+h 2^ calculated level. If the 2^ experimental level 
is really r. true two-particle configuration and if our 
comparison turns out to be correct on calculation of 
other properties of this level, this would mean that an 
experimental 2"^ level corresponding to 2^ calculated 
level is yet to be observed. However, we would like to 
take the blame on ourselves. For this nucleus, we have 
not been able to achieve even the correct renormalization 
for the 2 !^ level (in other nuclei at least this has 
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been possible), woull ratlaer like to say that our 
renormalization approach has failed for this nucleus. 

.For few excited levels we hare be-.n able to get satis- 
factory numbers but •■^’q make a reservation to s-'’y if the 
predicted parameters also repr- sent the -truth, ^^.1 though 
the 4*^ level is b--st reproduced in spectrum Y, we can 

say that spectrimi 71 gives the more satisfactory 

92 

agreement than any othev sot. The Zr nucleus also 

+ 

does not present a much happy situation. The 2^ level 
is best predicted in spectrum I but only at the cost of 
all other levels. Instead, at least correct ordering 
is predicted in spectrum II by Just a negligible sacri- 
fice in the fit of 2^ level.. The spectrum IV describes 
^he 22 level best but otherwise th.3 situation is 
unrealistic - even the predicted 2^ level goes below 
the ground state (not shown in figure). Overall, the 
spectrum V describes the levels best. 

II. 4 COnCLUSIOV: 

Thus, mathematically at least, we have been able fo 
represent the renormalizations required for few low-lying 

3 

levels in tems of the matrix elements of the 
realative state, 4 two-body matrix element is the final 
outcome of an interpl.ay among several channels of 
interaction. This statement also holds for the renormali- 
zing two-body matrix elements. Their representation by 
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rel^itiiv'o ]nqi;rix ©IsrrieTi’fcs alone is an ovens iinplif’ico.'tioii 
of things. As such, it seems difficult at this stage 
if some detailed stiidy of the renormalisnt ion could be 
made in the presentl}^- adopted simple minded approach. 

In conclusion we sev that our approach turned out 

to be much less successful for T=1 levels than for 1=0 

levels, for T=0 levels vvs calculated a f or • different 

levels separately and then found the best fit value of 

a. Thus WG nsec several levels for establishing the 

effective interaction. However, it can be seen that, at- 

2 

least for loveest j configuration levels, the values of 
a as calculated for different levels are auite close 
and therefore, even if we had not used several levels, 
predicted excited levels, even calcul'^tad with a as for 
ground state, would have been better described than with 
bare interaction. In case of T = 1 levels, the attempt 
of determining three parameters by using onlcr one 
experimental level turned out to be the root cause of 
troublvCS. We only tried to find parameters so that the 
ground state level energy is correctly reproduced. 
However, in this process, it was never guaranteed that 
other levels would also be reproduced. Therefore, we had 
to scan the whole range of possible values of parameters 
in small steps with a hope that/we might strike a parameter 



169 


set which, is else suitable for description of higher 
levels. Ill first run, when the p-^r^meters were varied 
in crude stops, to see if such epproech is going to 
yield some results ct ell, such close ngi'sement wos not 
found because of which, in the finnl run, parameters were 
varied in finer steps. As per methom'^t ic=>l results, for 
1=1 levels, we have been able to predict the renormali- 
zation in terras of potential parameters for few of the 

58 

low-lying states. Except for Ifi, wo hove been able to 
reproduce at least the 2^ level corroctlv, however, 
using only one level - the ground S'tate O"*”. 



CHAPTER III 


A PH E HO ME '•■TOIjO&I CAL EPPECTI'IC liWERACTlOl^l POR 

III . 1 I:MTR0DTJCTI0N’ : 

Of the two lauclei and ^yCl, the fact that 

energy levels of one could he derived if those of the 

other were known was one of the successes of n 1-coupling 

nuclear shell model in its infancy. In their ground 
40 58 

states, the K and Cl nuclei have respectively three 
and one protons in the 0*13/2 Both of them have 

a single neutron in Ofy/2 low-lying energy 

levels of these nuclei are expected to arise as a result 
of interaction between the Od^/2 <orbit protons and the 
lone neutron. Through the use of coefficients of 
fractional parentage, matrix elements of twr-particle 

rz 

interaction in the configuration (0d^^2 *^^7/2^ 

expressed as linear combinations of matrix elements of 
the same potential in the two-particle conf iguration 
(Odj/2 Ofrjr/2)* means that if same radial wave - 

functions are assumed for the two nuclei (which amounts 
to assiming same oscillator length parameter b for them), 
energy levels of one can be derived if those of the 
other are known. The lowest four levels in both these 
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nuclei have J = 2 , 3 s 4 and 5 . Por Cl, these 

can be considered as corning from coupling the angular 

momentum 3/2 o-p the last p^roton to the angular momentum 

7/2 of the odd neutron. If j .j - coupling is valid for 

the valenc'- nucleons, the throe protons in the 015/2 

in compliance with the Pauli princiole, can give a total 

•angular momentum of 3/2 only and then the lowest four 

levels of '"K. can also be considered as a result of 

coupling of the angular momentum 3/2 to 7 / 2 . A.ssuming 

these couplings, with a knowledge of excitation energies 

in nucleus, it should be, then, possible to deduce 

excitation energies in ^^01 nucleus and their comparison 

38 

with experimental spectrum of Cl v/ill confirm the 

.14 

above assumptions • G-oldstein and Talmi found th-^t 
the said transformation worked excellently. 

Mb' 

The Od-5/2 can at most acco' modate four 

particles and thus, the three protons in the Od^/g o^rbit 
are equivalent to a vacancy - the hole - in the sajne 
orbit. Thus, the problem of handling three Od^/2 
protons raid one 0^7/2 neutron is mathematically equivalent 

15 

to one of one hole and one particle. Pandya showed 
that the matrix element of interaction between one 
hole-one particle states can be expressed as matrix 
elements of same interaction between two-particle states. 
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The relation he erriTed elj 
<(3h^ Op) JT |Vl ( 2 ^] Dp,) JT >' 

( 1/2 1/2 T 'j 

1^1/2 1/2 T'^ 

<0 0 J' |V| D/ d., J' T' > (3.1) 

h' P ^ P 

,He else verified this relation for th ' c.-^se of «nd 

^^01 nuclei. 

Since then, the problem of these two related 

nuclei has been discussed by di-"‘ferent authors with 

30 

different interactions. Kuo and Erown c--^lcul!=ted the 

spectra with reaction matrix derived from HJ potential. 

59 

Dieperinlc et.al. used Tabakin interaction and the 

modified surface delta interaction (MSDI) and worked both 

1 35 

in TDA and RP4. hieperink and Brussaar^ have tried 

to parametrize the effective interaction in Talmi approach. 

Besides those specific nuclei, attempts havo also been 

made of determining the ■ d^y 2 - ^7/2 interaction 

matrix elements from studies in neighbouring nuclei. 

Mo inester and Alford ^ have determined the ~ ^ 7/2 

interaction using multipole sum rule methods. Sartoris 
57 

and Pjpmiek have shown that the discrepancy between 


= - 2 (2J' +1)(2T’+1) 

J'Tf 


^h ^p 


-1h» J 
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P'^r tide-hole matrix elements of nuclear intei’cctions 
e.s obtained by using e realistic interaction and "^s derived 
from experiment is partly but not coyiipletel3;" removed by 
considering the effects of core-polarisatj-on. They had 
proved this statement by using HJ and KK interactions 
and with the ^'^7/2 conf igurat ions as example. 

However, the paper was written in the year 1967 and 
today we knov/ that microscopic core— polarisation calculat- 
ions are not satisfactory. Hrne has discussed the 

problem of nuclei in the range "^S — Ca under the 
assumption of an inert S core -^nd, inste^’d of using 
Pandya transformation, ho has calculated the matrix 
elements for the 6)d^y'2 terms 

of *^^7/2 configuration particle-particle ma,trix 

elements which were then determined in Talmi approach. 

In the past few years^ several measurements of 
energy levels of have been mode and now we have a 

better knowledge of the energies and spin-parities of 
quite a few low-lying levels of this nucleus. In the 
previous chapter we h^^ve pointed out that the paucity of 
experimental information is the m^^jor obstacle in the way of 
determination of phenomenological effective interaction 
for any nucleus. As such, with the improved knowledge 
about the levels, prospects of establishing a reasonable 
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effective interaction tod^y are brighter . 

Johnson end Kennett have determined level 
energies ano deauced spin-par itie^s of levels in 
ITink -ind Schiffer have recently performed (d,p) 
experiments on ^'Cl and with 12 MeY deuterons. 

In the previous Chripter we have seen that 
introduction of both phenomenological and experiiiEiital 
inform^^tion into microscopic calculations improves the 
results and' also makes them economic in respect of cow'ljoie^j, 
time. However, a comparison of the results for T = C 
and T = 1 states shows it clearly that the improveirevot 
depends upon the amount of information put in. Basical^ 
a realistic interaction is also derived through pheno- 
menology. However, ths amount of available experimentohl 
data is so vast that in construction of the phenomenelogXc<x(l 
form of the free M-H intersection, /fe cen incorporate as 
many depet\dences as we wish, How, H is only speculation, 
on our part th-^t bv '•ddition of seemingly justifiable 
corrections vre can modify it properly so as to be usablvs 
for c-^lculation of oound state properties- Such an 
approacb h?.s ultimately been shown to have failed and, 
asy^iscussed, the semi-microscopic approach has also not 
met enough success unless sufficient experimental 
information is fed into the calculation. Then it see;-- 
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to us th-st this semi-microscopic approach, a hybrid, 
is no better than the purely phenomenological approach 
to nuclear structure except for the fact that in the 
former, at least a. part of the effective interaction 
is knovm to have some well-founded origin. 

In the previous chapter we hare used Sussex 
matrix elements (SME) under the assumption of their being 
composed of central, spin-orbit and tensor components. 
This interaction ?tlso, gives gross similar results as other 
realistic interactions do and a.lso , at least the study 
of spectra with the central, spin-orbit and tensor 
components did not lead us to any conclusion which might 
have made us feel doubtful as to the basic constitution 
of the interaction as regards its components. Thus, we 
feel th-’t the true free II- ST interaction, as represented by 
the SME, might bo possessing several complicated 
dependences as evidenced by existence of such terms in 
other realistic interactions, yet, effectively it can 
be understood as consisting of central, spin-orbit and 
tensor components. 

With aboye considerations in mind we believe that 
if parametrisation of a phenomenological particle-particle 
interaction incorporating the central, spin-orbit and 
tensor dependences could be possible, it would be as 
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good as a semi-realistic interaction. it least it 
would be designed to reproduce experimental energies. 

For that matter, even realistic or somi-realistic inter- 
actions also do not gu^^rantee wave functions. i purely 
microscopic appro-'^ch does not lead us to the true numbers 
and once we start introducing phenoraenology , there can 
be more than one ways, each of them being re^isonable 
and sensible in its own place but different forms of 
addition will result in different structures of ways 
functions . 

For the purpose of parametrizing a phenomenologic-^' 
effectiye interaction, we suppose, the - ^®C1 nuclei 
present a fayourable example. Now-, sufficient experimer . 
information about low lying leyels of is at hand 

so that an effectiye interaction, constructed from centra' . 
spin - orbit and tensor dependences of the 11-17 potential, 
can be tried to be parametrized. Enough experimental in- 
formation is not available about ' Cl levels. Hence, the 
effective interaction p-^rameters cannot be determined 

•Zp 

from known levels in ' 01. However, as pointed out earlier, 
in eq-. (3^1), the hole-particle Interaction is same 
as the particle-particle interaction Y for which ^he 
two-particle matrix elements, as appearing on the right 
of eq. (3w1), are calculated- Thus, once the parameteri- 
zation is done for levels, the energy levels 'of 
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are automatically calculated. ¥e take this relation 
between the two sets of energy levels r.s the test of the 
effective interaction parameters "'a "^infi the best set 
of parameters for and apply it to energy levels 

calculation in ^^Gl. 

III. 2 IffllHOD OP. CAICUIilTO'I: 

Bass and Wechsung"'^^ have tried to deduce the 
amount of configuration mixing from a study of electro- 
magnetic transitions between negative parity states in 
4-0 

'K and they find that the structure of low lying levels 
can be reasonably described in terms of relatively simple 
configurations. Their results are shown in Table 3.1. 

It, thus, becomes clear that, except for the 
excited 3~ level, the wave functions are almost pure - 
a result perceived iij Goldstein and Talmi much earlier. 

It presents a favourable situation Lur parametrizing a 
phenomenological effective interaction - at least the 
P'^rameters obtained will not bear the blame of unduly 
absorbing large configuration mixing effects of configur- 
ations lyi'ig outside the model sp^ce. Although a model 
space allowing the hole to be anywhere in sd-shell and the 
particle to be anywhere in fp-shell would be much desirable 
yet, in view of available information about wave-functions. 
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TLBLE 3.1 


Rel^tiTG Intensities of component configu'c-^tions 
40 

in the structure of ^ K lereJ s derived by Bass and 
1 40 

Wechsuiig 


Level 

energv 

(MeV) 


Intensities 


j’^ ( Od^y '2 0 ^ 7 / 2 ) ( ^^^3/2 "*^3/2 ^'^7/2^ 


"("in per cent ) 
1 ... . .. .-1 


0 

4" 

100 


0.03 

3' 

98 

2 

0..8 

2" 

99 

1 

0.892 

5" 

100 


2.048 

2" 

3 

97 

2.070 

3" 

3 

65 

2.104 

r 


100 

2.626 

0" 


100 
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we feel it would be reasonable enough to choose 
the smaller model space : 

hole states : 0d^y'2> 

particle states: OfY/2’ ”^^ 3/2 

In principle, unless we show in our case also that 
the effect of left out coiif igurations is really small, 
our assumption of this truncated space will rema,in 
questionable. However, upto the extent of nonaccountable 
effects of left out configurations on spectra, calculated 
with effective interaction parameters derived in a trunc- 
ated model space, we assume that the effective interactic 
will explain low-energy levels. 

Goode and Zamick"’^*' have proposed two sets 
of single particle energies in this mass region 
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SET I 

SET II 


(Me7) 

(MeY) 

Od 5/2 

-11.7 

-9.9 

IS1/2 

- 9-7 

-7.7 

Od 3/2 

- 7-2 

-5.4 

0f7/2 

0.0 

0.0 

IP3/2 

1-9 

1.9 

IP1/2 

4.25 

4.25 

Of 5/2 

6.25 

6.25 

®o' 7/7 ®oh/2 

7.2 

5.4 

Odj/2 " 1 ^ 1/2 

2.5 

2.3 

They observe that 

the spectra of 

and ^’’ca are 

insensitive to cba 

nges of 1 MeY 

or so in their single 

particle energies. 

-Although the 

results in "^^Ca are 


improved by decreasing the gap b3tvveen Ofyy '2 ^*^3/2 

39 41 

orbits from 7-2 to 5'-4 MeV, the results in K and Ca 
deviate even more from the experirreiital values-. Dieperink 

{T Q l.n 

et.al.^'^ have used a different set for Ca nucleus 
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T\BLS 3*2 


Expsrira-int.'^.l sitii^tion reg^^rding e'-xorgy levels 
40 

in K as per different available compilations. 


Ref. 138 

Ex 

(MeY) 


Ref. 142 

Ex 

(MeY) 


Ref. 139 

(MeY) 

4 

0 

4" 

0 

4" 

0 

3 - 

0.0296 

3" 

0.0296 

3" 

0.05 

2" 

0.8 

2" 

0.8001 

2" 

0,801 

5" 

0.892 

5" 

0.8916 

5" 

0.891 

2"(3") 

2.0473 

2" 

2,0474 

2" 

2.048 

3" 

2.0693 

3“ 

2.0697 

3" 

2.072 

r(2",3" 

) 2.1043 

(1)" 

2.1036 

(1)" 

2.105 

(2", 3") 

2. 4188 

4" 

2.3976 

4" 

2.397 

3 " 

2.4575 

2" 

2.4191 

2" 

2.420 

17 2", 374") 

2.6265 

(0)" 

2.6258 

(0)- 

2.628 


1 " 2.7301 

r(273“) 2.7562 
2"(3") 2.7862 


Contd 


Contd. T'^ble 3.2 
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2"(3") 

2 " 

3“ 

r 

3 "( 2 -) 

1 " 

r 

1 " 

3“ 

r(2",3' 


Ref, 138 

2.8077 
3.1283 
3.3671 
3. A3 37 
3.4857 
3.6636 
4.0198 
4. 1040 
4.2536 
4.4628 
4.5378 
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TA'BL'S 3.3 



Adopte 

d level 

energies 

and prospective 

domin^^nt 


configurations 

in the 

structures of re 

spective 


levels 

• 




Level 

number 

j-K 

Ex 

(MeV ) 



1 



0 

0 d^y2 

Of 7/2 

2 



0.0296 

0d^y2 

®^7/2 

3 


2i 

0.8 

Od ^/2 

Of 7/2 

4 



0.892 

0 

1 

ro 

Of 7/2 

5 


^2 

2.0473 

1 

0 

1 p ^/2 

6 



2.0693 

Od ^/2 

IP 3/2 

7 



2.1043 

0 d;^2 

IP 3/2 

8 


^3 

2.4575 

1 s ;)2 

017/2 

9 


^3 

2.4191 



10 


4-2 

2.3976 

1 Si ^2 

017/2 

11 



2.626 

CM 

1 

0 

IP 3/2 
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function, The assignment to levels 1-4 is unique. 

Excitation of the particle from to the 1p^y^2 

requires 1.1 l/eV cnsrg''^ while 2-3 r:,-2 are needed for 
exciting a hole frqm th- Od^/^ to tne 1s ^ /2 ortit. ¥e 
do anticipate a 1 level around ? i'iet energy as resulting 
from excitation of the p-'rticle to 1pj/2 while the 

hole stays in o^'bit. On the basis of this we 

— 1 

assign configurations to levels 5-7. The 04^/2 "^^3/2 

configuration will give a 0“ level also. On the basis 

of consideration of single particle (hole) energies, 

other 0“ levels arising from 1 s^y '2 ’'P-]/2 ^ *^‘^5/2 

Of^/ 2 ) configurations are expected to lie much higher 

in energy. Yfe assign this conf iguration (0<3.^/2 ”'^3/2^ 

to level no. 11. In the ascending order of unperturbed 

energies, next comes^he configuration 1s^^2 *^^ 7 / 2 * 

assign this configuration to levels 8 and 10. However, 

this assignment is not expected to be very satisfactory 

- 1 

since the unperturbed energy itself ofthe ^^ 7/2 

configuration is 2.5 Ke7 while each of the 3^ ^2 

levels appears at an excitation energy less th^n 2.5 MeV 
which is possible only if the particle-hole interaction 
in this configuration is strongly attractive so that 
the unperturbed m='trix elements of Hamiltonian for both 
these J states fall' quite lower than 2,5 MeT since 
configuration mixing will send them up again. Om^-may be 



186 


that there ire other 3 ~ and 4 ~ unidentified levels 
which are causing configuration mixing. Otherwise j 
except for th‘^ oossibility of some coTnolicated .confi- 
guration level getting down lov^ in energy which might be 
a possible c-^use of this conf igur'''tian mixing, there 
®re no other 3 '^nd 4” levels within the model space 

considered. Iven otherwise, the unperturbed energies 

— 1 -1 

of ”^^ 3/2 "^^1/2 *^^7/2 

1.9 MeY and 2.5 MeY respectively - quite close. As such, 
we anticipate approximation of assigning pure configur- 
ations to 32 ^h.d 3 ^ levels to be somewhat restricted. 

In compi^rison to the level energy, the 2 ^ level 

appears at much too low an excitation energy to be 
-ccounted for by some simple configuration like (Od^y'2 ^^1/2^’ 
(1s'^^2 '^^3/2^ *'^‘^5/2 0^7/2^- include it 

in the fitting calculation. As more and more dependences 7 ' 
in the effective potential are introduced, larger 
number of experimental levels 'ir^;. rv.:quirc'd. Attempt 
has been made to avoid use, as far as possible, of levels 
of uncertain configurations. 


For , 

1^1. 

Hence 

, using 

V 

1/2 

1/2 "j 

t 

^ = 1/2 


1/2 

1/2, 

I 
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^.nd 


n 


1 


1/2 

1/2 


1/2 1 

1/2 J 




= 1/6 


the matrix eleme-its of H'^’^'ilton.iaii be'^'-veer one particle- 
one hole configurat ions '^re given by 


<(3h'' 3p) JliH| (%l -1 ,) J1> = 


.-1 . 


0 + -p M .hi, hp.jp, 


- i Z (2J'+1) 

J' 


Dj, Dp 


J 


] 


1 


2 (2J'.+ 1) 
J' 


J . 

(. ^h ' '^p ' 

J'j 

jh h 

b 

[ h ’ ^p ' 

j’j 


V <1h»V'0|V|Di,Dp,J' 0 > 


< Di,,.1pj'l!vhi,Dj5J’1> 


(3.2) 

where experimontal values for single particle (hole) 
energies Bp (S^, ) are to be used. In the ground state 
conf iguretion there Is a hole in 0 ^ 3/2 ^ 

particle in OfY /2 experimental excitation 

energies are known as relative to the 4~ ground state 
level arising from this configuration. The gap between 
the Ofyy'2 *^^ 3/2 common to all levels and, 

hence, simpl?/ irrelevant for excitation energies 
relative to ground state. Like this, effectively the 
single particle energies needed in configurations 

(0d"/2 ■*^3/2^ 

respectively. 
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?irst we take the residual inter-action as the 
p mx* 0 c © Tiiu 3 ? 0 X f 0 12 ^ c ©' 


V^Cr) 


2T+1 , 2S+1 
= Z P 

T,S 




OS 


(3.3) 


where ’ 2S+1p pro:iection operator that 

projects out states of isospin T an'^ spin S, im is 
the strength in the particul=<r (P^S) chamiel of inter- 
action ana rQg is raage in the spin state S i.e. 
we assume an iso spin-independent ranffe. Matrix elements 

I* Q 

of V^Cr), as a function of \ (= — “ calculated 

V 2 b 

using the formula (A. 30) given in appendix /i. fhis 
is done in view of simple relation -^or Taloii integrals 
ihterms of X 


i >? \ P + 5/2 


1 +x 


(3.4) 


40 

r.m.s. charge radius of Ca is 3.5 
so that the corresponding value of b comes out to be 
2.013635 fm. We take b = 2.0 fm for Thus, after 

X is determined, the range can be calculated from 


■OS 


r? X, 


(3.5) 


Thus, there are six' parameters specifyijQ^g the 
central force. However, we have been able to reduce 
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the effective number of unknown parameters to four. 

This has been possible since we have an idea of the 
upper and lower limits for the nuclear force range rQ. 

We let Xg (for both 3 = 0 and S = 1 ) vary Independently 
between 0.3 and 1.5 in steps of 0.1 in first run. Thus, 
for a given ( Ao> number of unknown para- 

meters is reduced to fou’". In principle, four unknowns 
can be determined from a knowledge of four experimental 
values. However, in the present case we can find four 
.A-ijg parameters for e-^ch ( Aq, set. Each set of 

the Ijig parameters obtained in this way will describe 
the four experimental values considered equally well 
and it will not be possible to d istinguish the true set 
which will also be good for prediction of higher levels. 
To avoid it, for emery ( A q, minimise the 

function 



E 


i,Calc 




E. 

- , exp . 


( 3 ^ 6 ) 


for five excitation energies. The axe the 

calcul-.ted excitation energies corresponding to the 
respective experimental excitation energies 
Eor a given ( Ao» set, %^calc'^ linear 

functions of all the four A/p 3 's. The unknown strengths 

2 

are determined by imposing mathematical restriction of ^ 
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being minimal v/iih respect to each Arp 3 . 

4n experimental energy level carries in itself 
the effect of so many off-di-^gonal ra-^.trix elements of 
the Hamiltonian matrix. If the E^miltoni-n matrix is 
of dimensionality greater than two, it is not possible 
to write down an explicit analytic expression for per- 
turbed energy eigenvalues in te.rms of off-diagonal 
and unperturbed diagonal ma + rix eler^ents. tre, there'f^’ore , 
have to resort to approximate methods. We assume that 
in first approximation the excitation energy of a level 
can be approximated by the difference of diagon-^l matrix 
element of Hamiltonian in that particular state and 
the diagonal matrix element of Hamiltonian for the 
ground state 4” level, in the limit of pure configurations. 
Fortunately, for this particular nucleus, it happens to 
be a very good approximation except for the 2.0693 MeV 3 2 
level , (Cf. results of Bass and Wechsung^ . In this 
manner it was found th'^t the value of Aq around 
0.3 - 0.4 and thatof around 1.1 - 1.2 gave the best 

agreement for the five excitation energies under consi- 
deration. 4fter this, Aq was varied between 0.3 
and 0.4 and A-j between 1.1 and 1.2, both in steps 
of 0.01, and the above procedure repeated for determining 
the strength parameters. Thus, the best set of 

2 

strengths was chosen which gave least value of . 
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However, it was found th^t, for v^'rietion of 0.01 or 
so in ^min sensitive enough although 

the strengths predicted, for such two successive sets 
of Values of Ag, were different. Thus, it became 
obvious that small changes in the ranee can be adjusted 
through 9 variation in strengths, for the same reason, 
no fine'-’ v=u’intions in Xg were tried. 

However, as v/e have said, true energy levels 
require an exact di.agonalis.ation to be made and also 
otherwise there is no way to distinguish v/hich one of 

2 

the two or more sets, which all give the s-^me minimum ^ , 

is better, furthermore, a set of parameters may give 

2 

a least value of ^ in case of pure configurations. 

p 

However, -V obtained on diagonalieation with this 

2 

parameter set, mav not be the minimum obtainable -yL . 

Still, in view of the earlier observation that this 

nucleus is a case of pure enough con-figura>t ions we felt 

2 

that a parameter set that gives minimum yc diagonal- 

ization would not be very much different from the one 
that gives minimum for pure configurations. 

2 

Therefore, besides th-^ sets which gave minimima yC for 
pure configurations, wo also picked up other few which 
gave a yi^ value not very much different from 
The Hamiltonian matrices are then set up for all J states 

0 
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i — i 

in the model space ('ls^y'2 ^*^3/2 *^^7/2 ^^3/2^ 

each pareineter set seperetely and diagonalized- The 
2 

value of is calcul^^ted again for the S'^'me five 

excitation energies, this time with c-^lculated perturbed 
energy eigenvalues, with e'^ch pav'T-retef set. The 

2 

different parameter sets, although mav give s.ame 
for pure configurations, will contribute differently 
to off-dioLgonal matrix ele'^enfcs. This difference will 
show up in the value of ^ obtained for perturbed ' 

2 

eigenvalues. The set which gives a minimum value of ^ 
is thus taken as the best representative of effective 
interaction parameters for the five excitation energies. 
However, derivation of a set of parameters, which give 
the best description of the variables used in determining 
them, is not all that is expected of an effective 
interaction. Instead, besides the known quantities 
used in determining the paramete es , the interaction 
should also satisfactorily reproduce other quantities - 
In the present case this means that besides the "Pive 
excitation energies employed in parameterizing the 
effective interaction, the p^^rameter set should also 
satisfactorily predict the excitation energies of higher 
levels. It then turns out that the parameter set which 
p* ive s a least does not a.lso describe other 

excited levels equally well as can be don€ by scrificing 
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2 

the velue of by a small amount, finally, 

select only th-^t set which gives an overall best agreement 
with levels' excitation eaergies. 


4,fter the work with the central force is over, 
we proceed to determine the effective interaction 
parameters for a two-body potential of th'^ form - a 
central force + a two-body spin-orbit force. The spin- 


orbit force 


V 


TjS 


(r) 


2T+1,3 ^ 

S P ^ . S e 

T 


(3.7) 


contributes to spin - triplet states only. Therefore, ^ 
the contra! force parameters already determined will 
be no good for the central part of this combin-^tion 
since -addition of the spin-orbit part will modify the 
spin-triplet central effective interaction. Then, co 
restore the fit, the central part will also have to be 
modified. We assume iso spin (T) dependence of the 
spin-orbit force strength S,j,. The matrix elements of 
this force are calculated using formula (A.* 36) given in 
appendix 

Besides strengths, we have no reason to believe 
that the S = 0 and S = 1 central force range parameters 
Xq will also not change as a result of 

addition of the spin-orbit part. Instead, the modificatio]^^ 
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in central effective interaction miffht as well be 
reflected both, in ranve.as well as in stren.^th para- 
m.aters. However, it is rather firmlv established that 
the central part constitutes the na.ior and most important 
part of the two-bodv e-!"fectlve interaction with spin- 
orbit and tensor parts providing finer details of the 
spectra. This 'He have also seen in Chapter II in 
connection with our study of relative importance of 
central, spin-orbit and tensor parts of Sussex: matrix 
elements. 4s such, we assume th-^.t ).q and will 

lie between 0,3 end 0.4 =<nd between 1.1 and 1,2 
respectively and, as earlier, we let them vary in steps 
of 0 . 01 . The inclusion of spin-orbit part increases the 
number of unknown parameters from four to six and 
there.fore two more experimental energy levels are needed. 
The spin-orbit force is operative in spin-triplet states 
only and hence, the range of this part is also expected 
to be round about the range of central part in triplet 
state. Therefore, for the spin-orbit part, we 
let X , defined through 

= 2 iTo' X (5-8) 

vary from 1.1 to 1,355 in steps of 0.015. 'Por each 
set of Aq , and X ’ strength parameters 

are calculated through the least square fitting procedirre 
as describod abnye. 
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The Same procedure was repeated for par*omet— 
rising the e’^fectiv'o potential of the form central + 
tensor. The form of the tensor force is 


r 


Y 


Ten. 


= E 


2T+1 , 3 




3(^ .r)(o^.r) 


( Qj • ^9 ) 


•(rVr^_ ) 


m/ 


(3.9) 

where iso spin (T) dependence of the strength is 
assumed. Two-bodv matrix elements ot the tensor force 
are calculated usinr the formula (A. 40) giren in appendix A. 
The range parameter ' as defined by 


^Ten. 


= 2 f~2 ' 


(3.10) 


was varied from 1,1 to 1.205 in steps of 0.015. 

Still higher values of ' were not tried so as to 

2 

save computer time and the trend of ^ for pure 

conf igurations was identical with that for inclusion 

of the spin-orbit p=^rt and it was realised that with 

2 

such low sensitivity of % to ^ extension of it 

to still higher range of values will not bring any 
significant change in the final results on diagonali- 
c^ation. 

We need five excitation energies for paramet- . 
rization of pure central, force, ^^'e take the 4.^ 
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(ground st-te), 3 “ (0.0296 MeY), 2:j' (0 . 8 MeY) , 

5 (0.892 MoV), 1 “ (2.1043 MeV) and 0':j' (2.626 MeV) levels. 
For the Inclusion of spin-orbit ami tensor parts we 
consider the 2~ (2.0473 MeY) and 3^ (2.0693 lleV) 

levels also. 

.Yow we proceed to parametrize the effective 
potential in the form central + spin-orbit + tensor. 

Thus, there are eight strength parameters reauiring a 
knovifledge . of nine excitation energies, that is, ten level 
energies including the ground state. However, only 
eight excitation energies were ava.ilable (the 4^ 
level at 2.3976 MeY v/as established later). The '5^ level 
at 2.4575 MeY ¥/as -the only other, hitherto unused, level 
which could be used v/ith some reasonable degree of 
definiteness as regards the main component in the wave 
function. Therefore, v?e had to resort to one more 
approximation. We assumed spin-isospin-dcpendenca of 
strength and spin-dependence o-^ range of central part 
(four parameters). Wo also assumed isospin-dependence 
of strength of spin-orbit part (two parameters more). 

In analogy with the spin-orbit part, it would be 
desirable to include isospin-dependence of strength of 
the tensor force p«rt. However, we assumed that, on 
the basis of our previous knowledge that exact fixation 
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o.f' strength, and range is rather arbitrary and it 

is rather a ^oint contribution O'f the two which is 
parametrised, the iso spin-dependence of the strength 
of the tensor force can be absorbed by talcing ranges 
in the two isospin wstates to be different and assuming 
that the strength in both the isospin states is same. 

Lihe this the number of strength parameters was effectively 
reduced to seven although letting the range parameters 
for the two isospin state’s vary independently increased 
the computer time. However, the seven parameters could 
now be determined in a least square fitting procedure fr,5in 
a knowledge of eight excitation energies, lie used the 3^ 
level for fitting, besides the eight others mentioned above. 

Thus, instead of the tensor potential as expressed 
in the form in eq. (3.9) we use in the form 


^ 2T+1 ,3 ^ Q f 

^Ten- j ^ 


3( Sij.r )( ) 


I 


-- --p) 


/ 2 / 2 
-(r /r 


Ten„T 




(3.11 ) 


and correspondingly 


^Ten. ^ ^ ^ (3.12) 

The central force range parameters As 
varied between 0.3 and 0.4 and’ between 1.1 and 1.12 in 
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steps of 0.01 for S = 0 and S= 1 respectively. The 

range parameter fo-r the spin-orbit part 'vas restricted 

to three values 1.1, 1.115 and 1.15. In the first run, 

for A' assigned value 1.1, 1.115, 

1.13 and 1.145 in succession and it was found th«t 

was very sensitive towards the change of A'-j^ from 

1.13 to 1.145, becoming much wo se "for X''' “ 1.145. 

1 

Similar behaviour was observed in case of X''' although 

the pC ^ sensitive. There'^ore, from /\q~ 0.31 

onwards, X"' (both for T = O.and T = l) were varied 
/ T 

independently over three values 1.1, 1.115 and 1.13* 
Fixation of A.] ii^ 'the much smaller range (1.1 - 1.12) 
was intuitively done in viev/ of the earlier observation 
that with central + spin— orbit and central + tensor forms. 


WPS predicted to be 1.1. 

4-0 

Twelve energy/" levels for K are expected to 

arise v/ith in the chosen model space. Out of these, 

we have used nine for determination of parameters. Is 

discussed, a serious attempt of com.parison of calculated 

42 with experimental 4-2 level should not be made. In 

the energy region, around 4.5 Me/, where we expect 

1 

the I 2 e,nd 2^ levels arrising from the (■ts^y ^2 ''P3/2^ 

configuration to lie, the spin assignments are uncertain 
and there are several levels suspected to be 1 and 2 . 
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Hence 5 , a reasonable comparison of predicted energies 
of these states is not possible. 

Hoa'e'rer, ir the p-'ocess of establishing effective 
interaction for the tv/o-partid e matrix elements 

in the same model sp^ce as tbe one considered for 
particle-hole states are antomatically calculated and 
the parameters can, thus, be verified ggainst energy 
levels of ^^Cl, predicted in the model space 

proton : 

neutron: ®^7/2’ ”*^3/2 

38 

The experimental energy levels of Cl are 

Energy 

(MeV) 

2 ~ 0 

5" 0.67127 

3- 0.75526 

4“ 1 , 50887 

3" 1.61721 

3" 2.74310 

The lowest four levels are precisely the ones 
that would be expected .from the ®^7/2 
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However, appeerence of tv/o low-excited 3~ levels makes 

it much .doubtful if the simple two-particle model will 

work for them as it is expected '^o h- for the lov\rest 

four ®^7/2 ^o^'^igu.r^tion levels. The 3” level 

appears at an excitation energy of. 0,36195 MeV above 

the 3-j level wherces an excited 3” state is expected 

from the Od^y '2 ”*^ 3/2 configuration for which the single 

particle excitation energy above the lowest Od^y '2 ^^ 7/2 

37 

configurati 'n is not known since the levels of ' S and 
37 

01 are not established. However, taking same single 

particle energies as those taken for the single 

particle excitation energy of this conf igijration is 1.9 

MeV. On the basis of this we expect the excited 3~ states 

to hav.-^ complicated structure rather than .just two-particle 

states. We shall only see hov/ well a^re the lowest four 
38 

levels of 01 predicted with an effective interaction 

40 

parameterized for K. 

T O 

The ^01 nucleus is one with neutron excess 
and hence the w,avo functions do not have cure isospin. 
Instead, v/e write a total angularm omentum (J) coupled 
proton-neutron state as |a(p)b(n)J > where the proton 
is in single particle configuration a, and the neutron 
is in the single particle confi,guration b. This state 
is obtained from states of, well defined isospins through 
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=i(p) b (n) J > = 



u 


b J T + jabJT 


1 

°J 


Therefore, a proton - neutron interaction matrix 
element '.vill be 


<a(p)b(n)J |7| e(p)d(n)J> = -J jjCabJTi 7j cdJT 

+ < abJT \y\ cdJT >,p^Q 

Therefore, particle-particle matrix elemeats of 

T = 0 and T = 1 states (which have already been calculated 

40 

duririe; calculations for K) for the appropriate config- 

TO 

urations in the model space as considered for ^ Cl are 
picked up and the Haiailtonian matrices for different J 
states constructed by taking averages of corresponding 
matrix elements in T=0 and T=1 states which are then 
diagonalized to get the energy levels of Cl. 


III. 3 RSSULTS im PISCTJSSION: 

As discussed, six experimental levels were used 
for p'^rameterization of the pure central force and there 
remained three more levels which were to be predicted 
with the parameters which gave best fit (i.e. minimum ^ ) 
for the six levels used. Therefbre, it was for central 
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P'T.rt onlv th'^t the said six levels were best repro- 
duced for the (/^o»A^) set nt volues (0.39, 1.11) 
vdiile the overall egreeroent with experimental spectrum 
including the other excited states also v^as obtained 
for tip-' volues (0.38, 1.14). Otherwise, for paramet- 
ri'i^ation in tbo form, central + spin-orbit and central + 
toiisor, oiplit levels were already used for proramotri- 
■'3ation and there \\as not much point in seeking the fit 
for tln-i 3^ lev.,)l for reasons discussed above. As such, 
for pa ramatrization of the effective potential in these 
two forms, the particular set that gave an overall fit 
with the spectrum was identical with the one which gave 
a minimum for the eight levels used. For the 

par«me triaation in the form central + spin-orbit + tensor, 
all the- known nine levels 'were used up in parametrizing 
the pot'/ntj.al and hence, for this form also, there was 

2 

orio unique sot of p^r-imeters which gave minimum value of 'y^ . 

The '.iff active interaction parameters are shown 
in Table 3.4. The c'lcul ^ted, pure and configuration 
mixed I'iv.'ls for contra! potential as calculated with 
both the parameter sets are given in Table 3.5. Similar 
results for effective potentials in the other forms are 
presented in Table 3.6. The value of obtained in 

each case is also indicated. The same results are plotted 
in Figures 3,1 aM 3.2. The wave functions as obtained 



TA.BLE 3.4 


S.^o. 1 

Interaction 

Central j 

Centi'’'^! +' 

C-'ntrnl 
+ tensor 

Central + 
Spin-orbit 
+ tensor 


Parameters 

I 

II 

• i 

spin - 

orbit 

1 

Ao 

0.39 

0.38 

0.39 

0-'3 

0.39 

2 

Ai 

1.11 

1.14 

1.10 

1.10 

1*10 

3 

o 

o 

? 

-154.9 - 

218.7 

-694.9 

1337.50 

590.80 

4 

o 

o 

- 19.3 - 

20.0 

- 30.2 

19. i 

6,10 

5 

i^QdleY) 

-90.1 - 

89.7 

- 39.0 

-944.5 

-342.50 

6 


- 15.1 - 

13-5 

-4.2 

- 50.8 

- 35.90 

7 

X 



1.235 


1 .13 

8 

SqCMcV) 



-1.1 


-1.90 

9 

(MeV) 



-8.3 


11.40 

10 

A" 




1.115 


1 1 

S' (Mo7) 

0 




-1 .500 


12 

S! (Me7) 

1 

A'o' 

t f f 

Al 




5.500 


13 

14 





1.10 

1.13 

15 

(McV) , 





-0.3 

16 

Sj' (Mo7) 





0.3 

17 

roo (fm) 

1.103 

1.0748 1.105 

0.84851 

1.103 

18 

i’q^ (fm) 

3.1395 

3.2243 3.1113 

3.11130 

3.1113 

19 

^LS " 



3.4932 


3.1962 

20 

^Ten. n 




■ 3.1537 


21 

^Ten. ,0 (4rn) 




3.1113 

22 

^Ten, , 1 I* 





3.1962 
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with etLective potentials o-f" the foriPs central, 
central + spin-orbit, central + tensor anl central + 
spin - orbit + tensor are tabulated in Tables 3.7 - 3.10 
respoctivelv. The results quoted for pure central 
potential are those for the (^0’ -^1^ (0.38, 1.14) 

i.e. the one that- gives an overall agree'^ent with 
experimental spectrum. 

It is clearly seen from the Tables 3*5 and 3.6 
that addition of tensor force siibstantia.lly worsens the 
value of for levels predicted under the assumption 

of pure configurations. The energy of 3(J' state is 
0.02b6 MeV and, therefore, even small deviations in 
this level energy contribute largely to . The 

largest effect that coaf iguration mixing does to this 
nucleus is to reverse the ordering of 4.[ aud 3^ levels 
thus predicting the 3)J level to be the ground state 
which is not true. It is because of this factor that 

bocomes much worse on diagonalioation. It is only 
on the addition of the tensor force that the correct 
ordering of tlie 4),' and 3)J levels is achieved and the 
value of for configuration mixed levels is also 

the best for the central + tensor form. In view of 
this it can be said that presence of tensor force is 
a must as a part of effective interaction for this 



TABLE 3.7 


2H 


Coriipon.fci'its in the wave -functions for different 
Icnrels fitted with central force. 


predic:t&d 
excit-^tion 
e ne tp;j 
(KeV) 


5^/2 Of 7/2, 

I 

I 

®^3/2 *”^3/2 

1s^^ 9^7 

*“1 

1s.j ”^^3/2 

2.14804 

r 


0.99371 


0.11197 

6.03039 

1" 


0.11197 


-0.99371 

0.72810 

2~ 

-0.97716 

0.18853 


0.09801 

2.01898 

2~ 

■ 0.18011 

0.97963 


-0.08876 

4.99446 

2“ 

0.11275 

0.06908 


0.99122 

-0.04488 

3~ 

0.98605 

•0. 10762 

0.11975 


2.68470 

3" 

0.04635 

0.90223 

0. 1-2876 


3.78248 

3” 

0. 15419 

0.41762 

■0.89545 


0 

4" 

0.07335 


0.99731 


4.37134 

4" - 

0.99731 


0.07335 
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TA-ELIC 3.8 


Components in the wave functions for different 
levels fitted for cehtral + spin - orbit form of the 
effective potential. 


Predicted 
excit-^^tion 
energy 
(Me '7) I 

■'1 

1 

1 

Od^ Oi Y 

• 2 

1 

Od^^ IP 3 

^ ^ I 

IB^' Of 

1 s^^ IP 3 

2 

2.0'2948 

r 


0.Q8711 


0. 16003 

7.50552 

r 


0 . I 6 OO 3 


-0.98711 

0.69725 

2" 

-0.96142 

0.27505 


0.00606 

2. 16279 

2" 

0.27208 

0.95592 


-0. 12650 

5.22926 

2” 

0.04058 

0.11997 


0.99195 

-0.08 f60 

5" 

-0.97455 

0.22200 

-0.03680 


1.98059 

5" 

-0.19856 

-0.92511 

-0.32364 


4.15917 

5" 

0. 10589 

O. 5 O 8 O 4 

-0.94546 


0 

4“ 

-0.99851 


0.05807 



0.99831 


4.5162? 


4 


0.05807 
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table 3.9 


Components In the wave functions for different 
levels fitted for central + tensor form of the e'Pfective 


pote’T 

tial. 

Predicted ^ 

excitation 
anerpjy 
(MeV) 

1 

1 

Od.J Of^ 

3 ( 

7 2 

! 

-.■1 

Od^ 1p^ 

1s^^ Of.^ 

^ 2 

''P 3/2 

2.07357 

r 

- 

-0.9954 


0.0954 

5.23519 

r 


0.0954 


0.9954 

0.B0051 

2" 

-0.9704 

0. 1946 


0.1433 

1 .91236 

2" 

0.2089 

0.9736 


0.0925 

5.10337 

2“ 

0.1215 

-0.1197 


0.9354 

0.075B5 

3" 

0.9S33 

-0.1291 

0. 1286 


2.47890 

3" 

0.0324 

0.8183 

0.5739 


2.83760 

3“ 

0.1793 

0.5601 

-0.8088 


0 

4" 

-0.9863 


0.1649 


4.83220 

4“ 

0.1649 


0.9863 
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TAB LE 3»10 


CompoTients in the vvave function for different 
levels fitted for centr--=il + spin - ordit + tensor form 
of tl.G ef'fective potential. 


Predicted 

excitation 

enerf^y 

(MeV) 


Od,’’ Of^ 

0 / 

Oh/2 ■'h 

2 

1s^^ Of^ 

^ 2 

1 Si^ IP3/2 

■5 

1 .8828 

■ r 


-0.9602 


0.2793 

3.4736!5 

r 


0.2793 


0.9602 

0.71581 

2“ 

-0.9821 

0.0355 


0.1851 

2.36065 

2" 

0.0452 

0,9978 


0.0431 

4.80505 

2" 

0. 1829 

-0.0556 


0.9816 

0.00132 

3 

0.9B35 

-0.0341 

0.1776 


2.35687 

3" 

0.04S5 

0.9958 

-0.0774 


2.96056 

3“ 

0.1742 

-0.0847 

-0.9811 


0 

4“ • 

-0.9924 


0.1228 


4.50107 

4“ 

0.1228 


0.9924 
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nucleus, "to get correct ordering of the ^'^ 3/2 *^^7/2 

configuration levels e.t lecot, although correct excitation 
energy of tn.o level is not reproducible. 

It can be seen th'-'t the vSpin-orbit and tensor 

forces ^ct in opposite di^'ections for the 3 '^ level. 

Addition of spin-orbit to the cen.tr-'^l part further lowers 

the 3 ^ loV'.,il. On the other hand, addition of tensor 

P'-^rt to the central part pushes up this level to the 

extent that correct ordering of the and 3^ levels 

is achieved v/ith. the 3^ level appearing at 0.07385 Me7, 

A.ddition of both spin--orbit and tensor parts to the central 

part causes an elevation in the 3^ level position, from 

a V'^.lue obtained otherwise with the pure central force, 

which is much smaller than that obtained with the central 

+ tensor form. Thus, with the central + spin-orbit + tensor 

form, the 3"^ level is predicted only at 0.00132 Me7 

and the ovor'^11 agreement with exper Ln;untal spectrum is 

also worst for this form. The discrepancy may largely 

be attributed to the fact that the 3, level at 2.4575 MeV 

3 

had to be considered for determination of parameters 
and its appearance at a comparatively low energy makes 
the comparison of the matrix elemeut 

^^^1/2 ^^7/2 l^effl "^^ 1/2 ^^7/2 ^ 

itself with the experimental energy doubtful. 
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It is anticipated that the 3 levels involve 
a substantial amount of configuration mixing (as .judged 
from experiment'll le'/el energies and single particle 
energies being used) due to which it should not be a valid 
assumption th-’-t the theoretical excitation energies, 
as c'-'lcul'-'ted from unperturbed singl.e particle energies 
"^nd diagonal matrix elements of interaction, be directly 
equated to the experimental excitation energies, This 
is best illustrated by studying the wave functions for 
the 3 states as obtained with the central + tensor 
effective potential. The 3 'j|' level is relatively pure 
but the 3 p 3 ” levels are strongly mixed with each 

other. Such large configuration mixing is observed in 
none 0 " the other levels. Since the effective potential 
of the form central + spin-orbit + tensor has been 
paramotri.'aed under the assumption of pure configuration 

for tho 3 ~ level, theref'ore the wave function for the 

_ ^ -1 
3 level is mor 3 purf; in rouspect of the *^^7/2 

conf Lgura fcion than the corresponding structure for the 
effective potential of the form central + tensor. However, 
the 3 ^ level for tho centr.al + spin - orbit + tensor 
effective potential is in worse agreement with experiment 
than the 3 ” level obtained with the centra,! + tensor 
effective pote.ntial and we have seen that the spectrum is 
also overall worst reproduced in former case. 



220 


Enoirgy lovcXs of* ^^01 a,s cr.1c ul^tisd v/ith dlffsi'sn.’fc 
effective potent ir^ls e.s determined for K ^.re presented 
in Tflble 3.11. The results quoted for centrel effecive 
potential '='re tloe obt'^inod for the ^'^0’ ^ 1 ^ 

(0.38, 1.14). It is cleo,rlY seen that the 3“ level 
position is in worst ngreemont - the effective potential 
of any fotm predicts 3” the ground, state J instead of 

2"” as exparimontally observed. The 3 predicted level 
is closest to scro for the effective potential of the 
form central + spin-orbit. However, wo may recall that 
this potential was also the one which gave the worst 
agreement for 3-| level in 40jr. The central + tensor 
potential predicted correct ordering of 4^ and 3^ levels 
in but for ^®C1, this form predicts .a deepest 3 

level. These observations regarding the behaviour of 
3“ levels are iii 3 onf irmity with those of Finh and 

I 

Schiffer^^^ who find that an improved value of (01^/2 
two-p.article matrix element in ^'^Cl predicts a worse 
agr.,-ement for 3" level in Howovor, except for 

the 3" level, the ordering of the lowest 2 , 5 and 4 
levels is correctly reproduced and the excitation energies 
of 5' and 4" levels are in quite an impressive agreement 
with corresponding expo '’imeiilal numbers. Nothing much ohould 
be expected of the form central + spin -orbit + tensor, for 
re- sons discussed above. However, the form central + tensor 
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TABLE, 3.11 

Lovel 3xcit--^tiorL energies of Cl predicted for 
different forne of ttie effective notenti'^l with perametei'S 
ns those given in Table 3.4. for eny form of tte potential, 
set I of la '/els is for pure particle-particle configar^tions 
and set II gives level energies on diagonalization. To 
facilitate comparison with experiment, all excitation 
onorgies have been calculated relative to the 2"^ level. 


Expt , 

(MeV) 


I 

(Me?) 

Gentnai 
' II 
(Mo? ) 

Central + 
I 

(Me?) 

spin-orbit 

II 

(Me?) 

0 

2" 

0 

0 

0 

0 

0.67127 

5“ 

0.69494 

0.71470 

0.69512 

0.72939 

0.75526 

3" 

-2.72652 

-2.78045 

-0.83476 

-0.81810 

1 . 30887 

4“ 

1.32250 

1.34226 

1.32177 

1.35604 


0" 

-0.08350 

-0.06375 

1.11642 

1-. 15069 


r 

2. 10979 

2.12955 

3.08537 

3.11964 


2" 

2.86838 

2.90790 

3.05396 

3.12250 


3“ 

0.30610 

0.39955 ■ 

1-. 41411 

1.46599 







^ ^ ^ ... . 


Contd. Table 3.11 
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Contd. Tf>ble 3-11 


Expt . 

(MeV) 


Central 

I 

+ to:inor 

TI 

Cririijr'^ l+spin- 

I 

-orbit+tensor 

ir 

0.0000 

2" 

0 

0 

0 

0 

0.67127 

5 ~ 

0.51963 

0.54318 

0.37851 

0.37966 

0.75526 

3" 

-6.64042 

-7-05091 

-5.73825 

-6,32996 

1.308P7 

4" 

1.14031 

1.16386 

1.01757 

1.01882 


0" 

1 . 46004 

1.48359 

-0.85892 

-0.85767 


r 

0.09179 

0.11534 

0.23312 

0.23437 


2" 

3.02191 

3-06905 

2,58219 

2.58469 


3" 

0.51689 

0.98048 

-0.11107 

0.48314 
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•ilso does not reproduce the excitation energy of the 
5 level as well as is done by ce itrnl end centre! + spin- 
orbit forms. Best agreement, excepting the 3*" level, 
with experimental spectrum is found with the central 
effective potential. 

III. 4 CONCLUSIONS: 

Above discussion makes it clear that for the 
40 38 

Problem of K and Cl nuclei, it is possible to 

40 

par'^metrize effective potentials for the K nucleus in 

such .a way that besides those of energy levels of 

38 

'Cl are also satisfactorily reproduced. The maximum 

discrepancy is observed in 3~lev'els. The tensor force 

A.C) 

plays an important role _in K nucleus particularly for 

^ O ft 

the 3"^ level. The level energies in K ■•^nd Cl -are 
correlated in such -a way that an improvement in levels of 
causes those of ^^Cl to become 'worse. The results are 
worst for the form central + spin - orbit + tensor, the 
possible reason for which has also been discussed above. 

In viev>/ of this it can be said that tho p'^.rameters of the 
effective potential of the form central + spin - orbit + 
tensor are a forcible imposition of m^^thematics on 
physics - the oldest question on validity of parameters 
of a phenomenological effective interaction. 
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feeling thnt nuclear level energies are relatively 
insensitive to details of effective interaction - closed 
shell + two nucleons nuclei are the mos t widely studied 
ones acd it has often been found that two different foi'ms 
of effectiv^e interaction give ^Irpost equally good fit to 
level spectra. However, the vmve functions do depend 
upon the details of effective interaction. Energy 
levels arise as an interolay between effective inter- 
action and the wave functions (which have been derived 
by diagonalizing the matrix of effective Hamiltonian! ) 
There-fore, it is believed th't other properties of 
levels, e.g. their electromagnetic properties, provide 
a sensitive test to the real 'goodness’ of an effective 
interaction. Surprisingly though, they observed that 

'pseudo-experimental' data on level properties, even the 

4- i 4-3 

electromagnetic properties of nuclei Ps to 'Ps were 
excellently reproduced bv assumihg an effective inter- 
action and pure (^ 7 / 2 ^^”'^ configuration. Their results 
also simultaneously proved the non- uniqueness of effective 
interaction. It was observed that the effective inter- 
action could be obtained either as a mixture of central, 
spin-orbit and tensor components or as a superposition 
of an additional pairing force over the original 
potential, however, whose -both the S = 0 nnd S = 1 
strengths are multiplied by the Srame f.^ctor a(> 1). 
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Thus, tho two effect ivo interactions differed not only 
in effective potential parameters but nlso in the pheno- 
menology adoptc-id for the choice of effective potential. 

Wave it not fo->" the strong configuration mixing 
involved (50 per cent or more), results of Cohen et.al. 
were also one among those obts^inod by so many others 
who try to parametrize effective interaction for explain- 
ing experimental data. After all, an effective inter- 
action is designed to reproduce experiment-^l data in a 
model space. The fact that such large configuration 
mixing as occured in their calculations could also be 
conc..?aled in the effective interaction parameters was 
surprisi ng. 

In a problem like th-t of nuclear structure 
physics we can utmost establish an effective interaction 
and say that.it is the n<iture of the N-,N force - nf course, 
in the particular model space chosen. Otherwise, in 
a different model space we might a'^rive at a different 
0 P Pqq - t; int cract.iou. and so on. Therefore, it was 
felt that if the results of Cohen et.al. could be shown 
to be true in different mass regions, it would very 
much reduce computational troubles since one could 
then take just the lowest conf iguration or a small 
number of low-lying conf igurations as the model space and 
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pn.r'^met 'rize' the energy levels in this model space. 

'■’'‘■^ith these points in mind^ G-upta and Trainor'^'^ and 
i5ngel and TJnna'^'^^ did similar calculations for sd-shell 
a rid fp-shell configurations and found that results 
op Cohen et.al. did not have genera,lity over different 
sc is ot conf igurations and that nuclear properties 
depondrnt on a/ave functions could still provide a 
sensitive tr-st of true effective interaction. 

While these conclusions mav be somewhat academic 
since they are drawn from models involving only two 
subshells, they probably suggest a lesson for any study 
using a truncated conf igur-^tion space. 

In light of these observations let us discuss 

92 

ar-^in some of our results for the Zr as obtained 
in Chapter II. Let us consider only the central part 
of the SiVrS and analyse the structura of 0^ level. ' We 
haVv^ secii th'-cc the dominant conf Igur-i bion contributing 
to this level for set II and set III of single particle 
energies is just di'^feront from the one which makes up 
most of the 0^ lev'-'l with set I of single particle 

Hh* 

energies. With set I, the 0^ level receives dominant 

2 

contribution from (1d^y2^ configur-^tion whereas with 
both set II ^'nd set III, the (0gy^2^ configuratinn 
contributes most to the 0^ level.' In practice, one will 
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use the totsl interaction r-^.ther th^^n just the centre.1 
psrt . However, it should be recalled that we hed used 
set II ’,"ith the reasoning th^t single particle spectrum 

PO Q i 

of 'S'!^ can also be used for ' Zr whereas set III 

1 28 9 1 

vcas given by Cohen ' spocificallv for ' Zr nucleus. 
Th..ce ’'oro , if we consider only set I =nd sot III 



Set I 

Set III 


(MeY) 

(MeY) 

CM 

ir\ 

— 

0 

0 

2s 1/2 

1.22 

1.55 

1d3/2 

2.07 

2.70 

0g7/2 


2.70 


■we see that for both the tot-1 interaction and its 
contral p-rt the predicted structure of the 0^ level 
so crucially depends on sp-^ce dimensionality and single 
paotiel- energies being used whereas the energy of this 
level as calcul-ted with bare SBIB with set I and set III 
turns out to be 3.77369 MeV and 4.39721 MeY respectively. 
The two energy values are of the same order of magnitude. 
Suppose we can locate the 0^ level in spectrum of Zr. 
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Further, let us assume we ere able to make a best 
possible reno’^mal.ization O'P bare interaction.Unless we 
caji be defi.nite to say, there is no a priori reason to 
believe that each of these two onea.aies will not be 
clos ' onou'^h to the 8xporimv3ntal eue'^gy level. Further 
■assuming th^>t the renormaliz^’tion process does not change 
the basic structure of a wave function - as regards the 
relative dominance of main contributor configuration, 
we thus find that like this there can exist more' than 
one wave function - each of them being drastically 
different in structure from the other such that each of 
th.''*m gives satisfactory reproduction of experimental 
energies. We realize that such a, situation arises because 
relative to set I, set III contains one more configuration 
in till’ model space. 

What we have said is speculation on our part on 

the basis of assumption that the reno cma iization process 

does not change the basic structure of a wave function. 

It might also happen the other way that on renormalization, 

the ot level calculated with set III also receives 
■ 5 2 . . 

dominant contribution from (1d^^2^ configuration. 

Our experienC'.^ with microscopic calculations 
tells that every other day we have been discovering a new 
r ono 'nfla li z a t i on process and unless '"i consistent theory 
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is available, above doubts cannot be checked but they 

cert'^inly Caiso more fundamental question - how to 

properly/" choosc! a model sp'^ce? We observed the particular 
, -1- Q ? 

behaviour of the 0^ level of ' fr for set III of 

single particle energies and a mod-el space of four 

configurations. In e^ch of the model spaces we might 

find a safisfactorv effective interaction where the 

effectiveness is brought to take into account the effect 

of neglected configurations. However, not much faith 

should be put in an effective interaction which brings 

about fundamental changes in th:> structure of wave 

fun.ctions - it sounds reasonable to absorb effects of 

neglected configurations through effectiveness only so 

long as they are small. For example, for spectrum of 
92 

Zr it appears that one cannot do without including 
the Ogrj ^2 conf iguration in model space. 

Thera^-'ore, before one starts out for ,an ambitious 
prograaime of nuclear structure calculations, the dispute 
over the model sP'^ce to be worked in should be settled. 

At present there does not exist any specific prescription 
as to which configurations should be included in dia.gona- 
lization. Practice, adopted by almost everybody, till 
now has been to choose the model space as encompassing 
the major shell, the lowest partially filled orbit of 
which is occupied by the valence nucleons. In case of 
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phe nomo no logicr;! e.ffectivo interactions, the effect of 
neglected conf iguretions is forced to be -absorbed in 


por-nneters. la case of effectir.v interactions 
from realistic inter'-'ct Ions , these effects '^re 


derived 


calculated 


in per’ turbot ion theoi-y, though lately it has also been 
realic.ed tint the configurations which have hitherto 
been treated in perturbation theory should better be 
included directly in diagonaliaation to understand a 
truer picture of effective interaction . Lo ludice ot.all^’^^ 


hovre made such calculations. 


With this question in mind we undertake the present 
study. In philosophy, our study is not very much different 
from earlier Vv^orks on pseudo-nuclei. However, instead 
of pseudo-nuclei we take just one nucleus - ^^Sc. It 
has a proton and a neutron in o^t)it. This gives us 

a chance to study all the four spin-isospin components 
of a central interaction. The model snace generally 
.chosen for this nucleus is full fp - shell where 
single particle energies are : 

Single particle energy 



(MeV) 

C\J 

O 

0 

Ip 3/2 

2.0 

1Pl/2 

4.0 

°^5/2 

6.45 
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This is also a faTOurable situation that there are 
four configurations in fp-shell so that the effect of 
graduallY expanding the model space can be better 
iiiYcatigated in four spaces. The similar case can also 
be of sag - shell but in the mass region where 

valence nucleons are in one O'f" the orbits of sdg-shell, 
it is not possible to find a self-conjugate nucleus 
where both neutron and proton might be in same orbit 
in ground state so that both isospin components might be 
studied. Further ^ these nuclei ha^e neutron excess so 
that even if we piclced up a nucleus where odd neutron and 
odd proton were in different orbits, the states will not 
have good iso spin which will, hence, be projected. Above 
all, calculations are much more lengthy in higher mass 
region. Furthermore, the highest configuration contri- 
buting in fp-shell is ( Of ^^ 2 ^ with unperturbed 
energy 12.9 MeV which does not occur in any of lower 
major shell, sdg-shell and also not expected to occur 
in still higher mass region where single particle lewels 
are more closely spaced. The unperturbed energv 13MeV 
configuration is believed to be weak enough in contribution 
rss regards its effect on low-lying levels. Therefore, 
i-f the effect on low-lying levels of expanding the 

space from ''P 3 / 2 ' "*^1/2^ (*^=^7/2’ "^^ 3 / 2 ’ "'^1/2’ 

shown to be substantial, it will be 
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indicative of nonconvorgencc of effective interaction 
OS regards sp^ce dimensionality. Since, then it will 
be no guar'^'iitee that configurations lying ebove the 
Of ^/2 I'ivol reall"''' hove a small effect. 

arbitrarily choose tie parameters of a 
central two-nucleo'n potential with the only presumption 
that the level energies calculated with them do not 
turn out to be orders of magnitude off the experimental 
numbers so that qualitative results which we would be 
concluding afterwards do not turn out to be completely 
hypothetical. Then, keeping them within limits as 
allowable in this sense, we study levels arising from 

p 

the configuration as a function of interaction 

parameters and space dimensionality. 


IV. 2 METHOTi 0? CMOIJLA.TIOM: 

2 

The (0fy^2^ confi guration gives rise to 
following levels : 


^ .+ r7+ r:+ n'^ 

T = 0 : 1 , 3 , 5 , 7 


0 ^, 2 ^ 


4 , 6 


T = 1 : 


J 
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form 


We assame a central two-body potential of the 


Y^(r) = S 2T+1,2S+1 p , e ^ ^ o' 


T,S 


2T+1 , 2S+1 


(r/r^ ) 


where ' ’ P is th 5 projection operator which 

projects out two-nucleon states of isospin T and spin S 
and Aijg is strength in the (T,S) channel, define 


A = 




and 


and set 


p= (4^/^ 


P = 0.5 

which corresponds to X= 0.57735 which, in turn, 
corresponds to r^ = 1.632B fm which is acceptable 
as range of two-nucleon force. 


The two-bi.)dy matrix elements , for a gi^ren set 
of values are calculated using eq. ( 1 . 30 ) as giv( 

appendix 1 where the relative matrix elements «re cal- 
culated in terms of Talmi integrals Ip’s and the coef.i 
dents B(nl, n’l' ,p). The coefficients B(nl, n'l’,p) 


in 


are, taken from ref. 144 and Talmi integrals for 
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shape "^re calculate'''^ for the chosen value of X 


using 


= ^oX 


2 

2 _ (M)) e^~ 


Tit 


J 


p+1 




.2^ n ,;/3 , , .2- 




2P 


2p+3 


I Vit 

‘(2-^ 2hL±3_ ) I . I 


2p 


p P--1 


where , 


M= A 


and 


(^(x) = 


2X 


fF 


X .2 

/ g” dt 
0 


are error functions tabulated iu ref. 14-5 


The odd - ^ stapes are believed to contribute 
weakly in comparison to evon-^ones.. With this realiz- 
ation, peopl’C have ai.so oft-on succeeded in establishing 
phenomenological effective interactions \7hich act in 
even ^ states only (Serber force) or even in ,iust 
relativvi S-states. Including the odd ^ relative st-^tes 
just for the sake of inclusion adds at least one more 
parameter. ThoreforG, if in a problem it ca.n be shown 
that odd ^ states really contribute much less than even 
•£ ones it might as well be a good approximation to 
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altogether drop this pert of the interaction. Similerlv, 
it may also be possible that some of the conf igruations 
within the conventionally chosen model sp’^ce affect 
the l-vel energies very little. We, therefore, study 
leve;lv^ aric’ng from the configuration as a 

function of all t?ie four parameters in different 

model spaces. In such a study, the points raised ."just 
now can bo investigated. If the points can be established , 
it v/ill be reasonable to drop the less effect producing 
interaction parameters and model space configurations 
so as to make place for other more important effects. 

This has to be done in view o-^ the reality that only 
rather limited amount of experimental information is 
available for parametrization of effective interaction 
and hence best use should be made of it. 

For the nature of above study, wa should individually 
invf.’Stigato all th' eight levels of the (Ofyy' 2 )^ configur- 
ation as a function O'f’ the various parameters. However, 

to facilitate analysis , ive study the snlitting between 

2 

centroids of T = 1 and T = 0 levels of the (0fr^y2) 
c 0 nf i gur a t i o n 

AE . = - Eq 
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Z (2J +1) Ej ^ 

where E^p = _J [ 

E (2J + 1) 

J 

2 

For the (0frj^2'^ configuration, -J and T are 
related through 

J + T = Odd 

and Ej p is the energy of the particular (J,T) state. 
Thefexact calculational details '’re as follov/s : 

J* 

Suppose we are studying the effect of y.ariatlon of T=0 
strength parameters iQg* Then we fix; 

A^O = - 40 MeT and = - 50 MeV. 

Further, when we are studs^ng the e-pfect of 

component on spectra to assess its felatiye importance 

we fix Iqq = - 70 MeT and let yary between - 40 MeT 

and - 80 MoV in stops of 10 MeT. \t each of these 

parameters sets we calculate eaerge leyols arising 

P 

predominantly from tho (OfY/2' ' configuration and thonce 
E. It can be reali'^ed that v/ithin a fixed m.odel space, 
if only T = 0 interaction is being yaried, the T = 1 
centroid v^rill stay at a constant level -and vice versa, 
^fter these five runs, to study effect of Kqq component 
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fix Aq^ = - 70 MoV and let Aqq vary between - 40 MeY 
and _ 80 MeY as a.bove. After this is complete we fix: 


4oo = - 40 MeY and 


^01 


= - 50 MeY 


and Vary the parameters on exactly same ranee of 

values as for Oqq parameters, lihe this when variation 
of all the four strength parameters is studied we 
repeat these calculations in a different model space 
U-sinv the same values of streigth p'^r<a_n’ieters . Thus, we 
study the effect of configuration mixing on each of the 
lowest eight levels of = 0^, 1"^, 2^, 3'^, 4'^, 5^, 6^ 


and 7"^ which arc expected to arise predominantly from 

p 

the (Ofy^g)''” configuration. 


The general nature of the t wo two-n.ucleon 
interaction can also be discussed from behaviour of 
with respect to the various parameters. Plowever, some 
lypical behaviour of can be subanalyzed in terms of 

behaviour of individual energy levels. 


Wo have said that when variation of E in a 
given space with respect to one isospin component is 
being studied, the centroid of the other isospin levels 
stays 'stationary and is, thus, redundant. However, study 
A® a' given set in different spaces is 
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expected to reveal the relative hehaviour of the two 
isospin components in the different spaces e.g. if on 
expanding the sp^ce, the splitting A increases, 
centroids of both T = 0 and T = 1 levels will no 
doubt go down, this ’would mean that the centroid of 
T = 0 st'^'tos is lowered mo'ce th^’n that of T = 1 states 
i,e. the new ^’dded configuration affects the T = 0 
levels more than T = 1 ones. 


Within fp-shell, biggest matrices occur for 
2 and 3 states - each of dimensionality eight. We 
ha made exact diagonalization for calculations within 
fp-shell. tIowGver, for later part of the calculation where 
we include still higher conf igurations , we do not ma’'<:e 
exact diagonalization. Instead, we C'-^lculate tite effect 


of added conf igur^-tj on in an approximate way as follows 


Fcrr a. 2 

■K 2 

Hamiltoni'’'n matrix 


/ 

V \ 


I ' U 

-12 \ 


Vh2 

^ 2 ) 

ssuraing -j < H 22 j 

the lower perturbed 


+ H22) - 


2 
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Thus, th .0 effect of perturbation i.o. the of ^^-diagonel 
matrix eleraent on ^ is 



-2 ,,4 

'H 2 ^12 


Therefore, if the oft-diegoml matrix element is small 
and the sep^rntion between the tv^o diagonal matrix 
elements Quite large, app'i^oximately 


or 








- H22 


Thus, the effect of exact diagonaliz^'tion can be 
approximately treated through calculation of 

h2 

*12 » 5’or the problem of nuclear enereries, the 

Ihrh2l 

diagon'^l matrix elements contain the unperturbed excit- 
ation energy of the two-particle configuration and the 
diagonal matrix element of the two particle interaction 
in excited configuration. However, for a highly excited 
configuration, ■ the denominator may be approximated by 
just the sum e of single particle excitation energies of 
the excited configuration. 
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Such a simple derivation is not possible for 
a many-levol system. However, to first approximation, 
the e'^fect of diagonaliz'^tion can be estimated by 
calculating the second-order contribution of laach off- 
diagonal matrix element above and, for a given 

matrix, adding all of thorn. Thus, due to configuration 
mixing arising from various two-particlo conf igurations 
I ,j >, the modified energy of a configuration | i > is 
obt'^ined through 

, , , i< i iTi j > d 

/\. = < i Y i > - S — 

" ,V1 e. 

where e • is the appropriate approximate energy denominator 
«] 

for the configuration j ,j> as discussed above. 

I < i I vl i > I ^ 

In what follows we shall call ! — L— l! — L 

e . 

D 

the second-order contribution of configuration > 
and denote it by S. For calculr^tions within fp-shell, 

Ini’s for all relative P, states i.e. 0 to 6 were 
included. However, to make things speedier and easier 
in calculation, for calculation of inclusion of higher 
configurations, only £= 0,1 and 2 were considered. 

It had been checked that the effect of inclusion of 
higher partial waves would be additive on3-y so that the 
second-order contributions as reported here will increase.. 
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An. idea of relative importance of relative states 
of -£ > 2 can be had from an example. In the 12*^ 

O 

state arising from (Oi.^^y' 2 )^ conf igurarion, 12 percent 

of th.e tot:^! of the spin-ainglet aiid 20 per cent of the 

total of the spin triplet contribution to the total 

0 

matrix elorrient comes from ^ > 2 p.artial waves of the 
interaction. We have -also made similar an-ilvsis for 
6 "^ levels arising wathin the fp-shell i.e. from the 
conf igur -■'tions (Of'jy' 2 ^ *-^^ 5/2 *^^’ 7 / 2 ^* assume 

''mo ~ ~ ^'11 ” MeY. The total matrix 

clement for 6 state in the (OfYy' 2 ) configuration is 
-I .5422 MeY out of which -0.2648 MeY i.e. around 
15 per cent comes from contribution . of ^ > 2 waves. 
The perturbed 6'*' level with inclusion of contribution 
of all values is at -1.5747 MeV and the splitting 
between the tv\fo 6"^ levels is 5.6406 MeY. Ignoring 
the contribut i.on of > 2 'waves, the nerturbed 6"^ level 
is predicted at -1.32095 MeY with a splitting of 
5.5378 MeV between the two 6"*" levels. 

IY.3 RhSTJITS INh mSCIJSSION: 

First of all we see the order of magnitude of 
level energies for a randomly chosen set of parameters. 
We choose 
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■^‘10 ~ 

- 40 

MeV 

^11 

- 50 

MeV 

O 

O 

II 

- 70 

> 

CD 

■"01 

- 60 

MeV 


2 

Then, for pure (0fyy2) configuration, excitation 
energies in respective sets of levels of the two 
isospins, relative to the lowest level of that isospin, 
are tabulated in Table 4.1. We see that even v»rithout 
any reservation towards strength parameters, the T = 1 
excitation spectrum resembles the experimental spectrum 
very much. The discrepancy is unignorable for T = 0 
spectrum. Sirail'-'r situation was nlso encountered 
while working with Sussex matrix elements in Chapter II. 
However, our aim is not to do parameter fitting and 
except for the splitting of 1 and 7 levels, 
numbers for 3 aiid 5 levels are tolerable. For the 
same set of par^^meters we also calculate the excitation 
spectrum relative to the ground state 0 level. We 
present the comparison of predicted spectrum with 
experimental one in Table 4.2. Thus, except for 

“f“ ”4“ ”4“ ”f” 

inversion of 0 and 1 levels and of 5 and 2 
levels, the level ordering is correctly reproduced. 
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T’lBL'B 4.1 

Excitation energies, calcul=ite4 separately for 
T = 0 ->nd. T = 1 wsets of levels, arising from pure 
(0fY^2)^ conf igur'^tion for '\^q = - 40 MeV , -50 MeY , 

^00 = - 70 MeY, = - 60 MeY. For comparison, experi- 

mental excitation energies are also given. 



T = 

Calculated 

(MeY) 

0 

Experimental 

(MeY) 


T = 1 

Calculated 

(MeY) 

Experimental 
■ (MeY) 

1-*- 

0 

0 

0^ 

0 

0 

7"^ 

1.1346 

0.007 

2-^ 

1.8333 

1.59 

3" 

2.0548 

0,879 

4^ 

2.6441 

2.81 

Y 

2.5899 

0.899 

6-^ 

2.8154 

3.24 
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T^.BLE 4.2 

Calculatod excitatioa spectrum of level 
arising from {Of conf iguration for = ■ 
- 50 MeV, = - 70 MeV and = - 60 


j"”^ Experiinental Calculated 

(MeV) (Me 7 ) 


6-^ 

3.24 

2.8154 

4^ 

2.81 

2.6441 

2+ 

1-59 

1.8333 

5 ^ 

1.51 

1.9257 

3 " 

1.49 

1.3906 

7 ^ 

0.618 

0.4704 


0.611 

-0.6642 

0 ^ 

■ 0 

0 


.40 MeV, 
Me?. 
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It was checked, that the ordering of 2"^ and 5*^ 


levels can be correctly reproduced by taking 
instead of - 60 MeV, however, the O"^ and 1 


Aq-| = -63 MeV 
levels 


still remain inverted. 


The results for A E as a function of Aqq, Aq-, , 
A.]q and A..j ^ are plotted in Figures 4.1 - 4.4 respecti- 
vely. In each figure, results of calculation for the 


variation of that particular 4ij,g 
model spaces 


parameter in the four 


I 0^7/2 

II Of 


III Ofy/2 '’P3/2 *'^1/2 

IV 0^7/2 '’P3/2 "^^1/2 '^%/2 

are shown. 


In any of the spaces, when varied over same 
range of values, A ® found to behave most sharply with 
A.q^. k strange thing is observed to happen in Figure 4.2 
for = _ 80 MeV. In all the spaces II, III and IV, 
coiif iguration mixing increases fast with |Aq^ 1. Upto 
= - 70 MeV, the lowest levels of all J's of T = 0 



-40 ~ ' -50 -60 -70 -80 

Aoo(M«V) 

Fig. 4.1, AE plotted as q function of Aqq in different spaces 
.■ . , (denoted by I,» dfs' e:tc'. see the text) 




¥ / 



-80 


-60 -70 

A^o (MeV) 

I, of in different 


spaces 
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can be identified as 


(Of Y ^ 2 _ ) 
the l"* 


component . 
3"^ and 7"^ 


containing a predominant 
However, for .'Iq^ = - 80 MeV, 
levels still receiving dominant 


2 

contribution from the configuration, the 

5^ level receives dominant contribution from the 


(IP 3/2 ^^ 7 / 2 ^ configuration and it is the 

2 

which comes predominantly from the 


5^ level 
conf iguratibn: 


This pushes up the centroid of the T= 0 levels 
because the 5*^ level occurs with a weight 11. There- 


fore, plotting a, smooth curve between 4 q-j = - 70 MeV 
and - 80 MeV only on the basis of extrapolation of 


graph between = - 40 MeV and - 70 MeV and the 
value of Ae at Aq.| = - 80 MeV is not correct. In 
fact, it would require a point-to-point calculation of 
/\ E in the region of values between - 70 MeV and 

-80 MeV. Comparison of graphs for variation Of ■ i and 
strengths clearly demonstrates the relative 
importance of ;lg^ parameter over 4 qq. In fact , /s, E 

becOiTies more and more sharp function of as the space 
is expanded and, opposite to it, mixing seems to 
smoother the behaviour of A E with Iqq. In space IV, 
the change in value of A E in going from = -40 MeV 
to /Vqq = - 80 MeV is about half the corresponding 
value in space I which itself is much smaller than the 
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variation in in going from Aq^ = -40 MeV 

to - 80 MeV even in space I. Therefore, it can be 
concluded that if it is planned to p'^ra.metrize an 
effective interaction for the fall fp-shell, it is 
q.uite safe to ignore the \qq part altogether i.e. 
consider a centr'^1 interaction contributing to even 
states only (i.e. a Berber force). Hov/ever, the 
discussion here is restricted to comment onlv on the 
central part of effective interaction. Study of Figure 
4. 1 reveals th^^t the T = 0 interaction does not give 
saturation of ^ E as regards space dimensionality 
within fp-shell. At larger 1 Aq-]| values, Figure 4.2 also 
reveals the same thing. 

Behaviour of the interaction components 
I in different sp-'='ces is more complicated than the 

I Aqq I components. From Figure 4.4 it becomes clear 
that expanding the space from TI to III does not bring 
in anything new to splitting A®* However, it will 
be unwise to jump at conclusion that the ^'P ^/2 
conf iguratlon is just ineffective. Instead, the 
expansion from space II to III should be analyzed in 
terms of behaviour of individual J levels with respect 
to the said expansion. In Figure 4.5 we have plotted 

p 

the- levels arising predominantly from (Of ,^^2) 





f:^5h 


I 



I 

4 . 5 - 



fcv«ls ark-lnq predomfnafly from (ofy^/ configuration p 
•ctivii of space d)men?!onaiity(dcnoted by 1,11, etc ^5ce' the 




254 


configuration as calculated in the four spaces I, II, 

III and IV, including configuration mixing effects in 
the latter three for Aqq = - 70 MeV, kQ^ = - 60,MeV, 

^10 ~ “ 40 MeV and ^ - 50 MeV. A. studA/- of Vigure 4.5 

clearly establishes that the effect on levels, of inclu- 
sion oe the '^P'i/2 ^o^fig'^i'ation in model space, is 
negligible and one might as well work without configur- 

p 

ations ^'^P'l/2 "^^3/2^ and (lp-]y' 2 ) 

for T = 1 levels although it may not be said for 

^^^^5/2 "'^1/2^ configuration. Also, the effect of the 

+ + 

1 p^y 2 configuration on 1 and 3 levels is substantial. 
Thus, at least few of the configurations may be ignored. 

We make a relative study of curves for 

i 

variation of K-isI strengths (^Figures 4.3 and 4.4') at 
the two extremes of the range of values of lA^g| 
parameters considered, tabulate the corresponding 

values of A E iu Table 4.3. 4 relative study of the 

results, for difference between the values of A E at the 
tv/o extremities of the curves in different spaces, 
for the two spin components of the T = 1 interaction 
shows that variation of ^ (odd Jl ) causes sharper, 
change in AE than does the 1 .,q component (even ^). 

This is for (0fr^^2^ configuration (space l). However, 
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'l ABini; 4.3 

In the fi-rst two colu’ons the parai^ebers are 

indi-cnteh- I, II etc. stand for space I, space II etc. 

In each o the spaces, the valu.e of A fi corresponding 
to the particular ' set is shown. ■■'M.thin the s^rae 
space, the change in value of 'is on^ goes -rroni 

-40 heV to - 80 MsV for the particular spin coniponent 
of the force, is also indicated. 


■‘10 


-70 


iVll 

(Me 7) 

I 

(Me7) 

Splitting 

II 

(MeT ) 

III 
(Me 7) 

IV 

(MeV) 

o 

1 

o 

• 

CD 

o.eoj 

1 

0.56 

I 

I 

0.28 7 


>0.95 


^ 0.98 

1 

1 

I 0.92 

j 0.72 

-80 

-0.47 I 

CD 

NA 

• 

O 

1 


-0.36_, 


-0.443 


-40 


O. 30 I 

-80 

-70 

- 0 . 39 ] 


0-69 


0,50 


-0.35 


0.85 


0.55 


-0.37- 


' 0.92 


0.60 


-0.60 


1-.20 
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this difference in behaviour is reduced ■^s the sp^^ce 
expends, to the extent th-=it in. sp‘-^ce If, .a^q causes 
sb'^.rper v^^riatln’-' than does the co-irorent mean.inp' 

tliereby thot , effectivel''?' , in space li even states 
aTs dominant o\rer odd ones. This is also the general 
belief that contribut.ion of even Jc states is stronger 
than that O''^ odd ones. Therefore, it v/-= hold it, vv^e can 
say that, for invastigation of T = 1 levels at least, 
one should adopt the full fp-shell as model space. 

Further, even in the sp'^ce IV, although the even states 
dominate over those of odd , the contributions of the 
two types 0 ri states a !’e competitive so that unlike the 
c-’se for T = 0 interaction, it will not be justified 
to UPC; a Serber force for T = 1 spectroscopy. 

For a given set ot interaction parameters, energy 
levels plotted a function of space dimensionalitv 
(Figuie 4.5) clearly show that as the sp^ce is expanded 
from I to IV, at lea.c’t the O'*" and l"*” levels continuously 
keep on moving down. Rather, fo::- these levels the effect 
of including the Of^/o configurati on in model space is 
much more than that of including the 

It moans that if an effective interaction is parametrized 
to fit the excitation spectra relative to ground state, 
we will end up with different parameters in different 
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spaces. This is, however, an established fact that 
phenomenological effective interactions are space 
dependent but the importance of the configuration 

puts a. doubt over physical validity of parameters derived 
in space 11 or III. This is because we feel the effects 
of neglected configurations should be meaningfully 
absorbed in parameters only if these effects can be 
proved to be reallv sraall. Otherwise, as in the present 
problem, say for 1*^ level, we are seeina: that the Of ^y '2 
configuration is making a sharp effect on energy for a 
random set of parameters. An effective interaction can 
be derived in space III ^which might satisfactorily 
reproduce the l"*" level but it can ver]'' well happen 
that if the 1"^ level position is calculated in space IV 
with same 'Effective interaction parameters, the agreement 
would be l.ost. It onlv means that wo ra=>v get the correct 
energy but not the correct ivave function. 

We have discussed the importance of ^^5/2 
We have seen th-^t for 1"^ level, this configuration 
affects much more than any of other low-lying conf igurations . 
This further raises a doubt that similarly the effect of 
some still higher configuration might be too strong to be 
treated in effective interaction pa-ameters. .As such, 
it is appearing that whichever configurations contribute 
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strongly should bo inoludod directly in diagonalization 
so as to be explicitly included as a component in tb i 
.’.'■avo f'l:: nc ti ■'n. Jt, therefore, b cc'" -s nacossarv to 
c''vci-c th’O conrerganc..' of contributions of higb';r lying 
CO if i"tir ■■ tions os the space is expanded. 

'..s mC'n''ioned oo'"’l.ier, we hove not made "'xact 
d i.-'go'.io 1 i ■-•atlon to c^nsidi-T the, effect of configurations 
Ivin-'’: btJyond fn-al'iel 1 . Hovvovi-r, tlien, we should first 
c*''‘ck as till 'now close is the appro Kim at ion of replacing 
t^U‘ ctjrt'cet diayonolised eigenvalue by the one calculated 
-..'nrOKlmatol’'' through colciil-tion of second-order contri- 
h’,'.ti,i)ri S. In Table 4.4 we list the second-order 
CO 'it ri butionn of dif f'orent conf igurations within fp-shell, 

! ■ t?i.' sc-con'-ordar cslcul'itions have been done for 


'"Of) 

= - 70 

May 

'"01 

r-- _ /JO 

K/IeV 

*10 

= _ 40 

Me/ 

■11 

= - 50 

Mo\r 


Tht- approximate level energies obtained in this 
way are compared against those obtained by diagonalization 
in different spaces in Table 4.5. We sea th-'t approximation 



259 


TmL^. 4.4 

’''i-r c"'Cul-^tion O'P soponrl-oriier contri- 


h'ltiofiJ"! O'f' c Tfif li 

i n ri , i n '•>; : ■> 1 ii ;l ii i:'- 

ijri-s P.''y)Vf' l0V/.,?3t (Ofyyp) 

p-f ? \ ■ ■:■: 11 n nd eruit r ’’ b , 1 1 i iit^ to 

cuaf ieur- 

O”^ =ind l"** 

,L • • / 1 . 'T h 'o r r o ■ it 1 si 

'-10 = - 40 IWO’/ n I,i 

■.'■'01 • 

= - 50 

= -70 1 

Q., = -40 WsY, 

•n; '’ixcitefi oonPi,'^- 

nr-’tiori 

1 > 

IT’iporturberl 
oxc Lt'-'t.Lop 
e rie 
(WcJV) 

< 0 f'f'lv 1 

^ (MeV) 

' Second-or- 
der contri- 
bution (s) 

(Me 7) 

( 1 P 3 / 2 O 

4 

-1.1396 

0.3247 

( iPi/jP' 

8 

-0.5348 

0.0357 


12.9 

-1 . 1576 

0. 10335 

(lp- 5 / 2 )^ 

4 

-0.8280 

0. 1714 

’'^^ 1/2 **^3/2 

6 

-0. 34^^ 

0.01977 

ClPi/2)^ 

8 

0. 2198 

O.OO 604 

Of^/2 ^^j/2 

6.45 

-0.9198 

0.15117 

Ofg /2 IP 3/2 

8.45 

0.43B2 

0.02?-72_ 

(Of5/2)"' 

12.9 

- 0.4124 

0.013184 
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T\B L0 4.5 

Compnrlso'i of incorpor"^ ' i jn of ef^oct of 

^ . H- -f 

co)'i ^ LtfUf'- 1. l.m mixin'? on. n '^nl 1 orising 

O 

pr ’rtor,]in-'nt ly i’rorri to,;. ( 0 fr 7 y 2 )'^ configi^ration, os 

c" 1 bv OK-'ct di"'rron"’l. i'^-^tion -^ni bv c'^lcalotlon of 

sc'coti '-cyrdor contrlbut .Ions. 


Sp'^Hh- 

’"Envirgy IK 

( I 2 ^ ~ 

corifigurotion 

(MeV) 

Enor.'yf by 
exact diogo- 
n"’].ia<^tion 

(McV) 

\ppr 0 ximotel 7 

calculated 

energy 

(MeV) 

Tl 

-4.3576 

-4.6014 

-4.6523 

ETT 


-4. 6054 

-4.7180 

rv 


- 4.7915 

-4.8218 

i.r 

-4 . 1 246 

-4.7561 

-4.2960 

in 


- 4.3501 

-4.3218 

IV 


-4-5004 

-4.4889 
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works much better for 1*^ level where the predicted 

opcroximctc euersiy fells above that obtai’ced by diogonali- 

'■^etidri la both sp-^ce III ancl IV. It means that if second- 

order cfi^dct 0 conf iyuT'^t ions boycnd. fn- shell can be 

shtAvn to be subs+'antlel, it wi].! imean that the e'^fect of 

those conf Lffuretions will be still more on dlagomlisabion. 

however, ^or lev:;! the secoad-order contribution 

ovwesti^atec the effect of mixing. Yet it does not 

mean that third-ordc" contribution is important. For 

2 

oxemrU.), the second-order contribution of 
oonfigur-vLlon is 0.^247 MeV and the next correction is 
oriLv 0.026 Mo7 (to be included with an opposite sign to 
that of the second-order contribut ion) which is really 
na j-lifTlblo and does not much account for the discrepancy 
b'‘tw('!en th.’ love! acergias in sp^ce II as obtained by 
d iavnna 1 i'sat ion and by approximrote ca I cul-^tion. Instead, 
t>hn di. scrnpancy might do attributed to underestimation of 
energv denominator which mi^ht be more than just the 

O 

o.xoitatloii ahorgv of (iP^y'2^ com igUi’ation. The 
discrepancy is observed in all spaces. However, that can 
be seen to be because o-^' the discrepancy in space II. 

Sinc<'-' the energy by diagonalization in space III is 

2 

-4.6054 MeY and tho second-order contribution of (Of5/2^ 
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arid 0.13117 Me? respectively, Sir'^ilsrl'?' , the coitri- 
butioiiwS 0" (lp-iy'2)^ snci ('^^5/2 ”'^3/2^ coof igurstions 


1 be Gonfipcrer] . Suppose we p'-'ranietri 'iie effective 
Int G "'cctlo ri in srece' TEI, Ihoc vve would be including 
»il c inf i-’-uret ions upto (113^/2) ot 8 Mo? without 
re's] ini th'^t coufiyur-^ tion 

if- cent hut j.n,"' about four tinier tVi ' 8 iWeV configuration 
'I I'i if b-; ' n-T >e,xcladed from riiodel sp^tco iust because 
the 6)f,-^2 oonfiffurc tion is not betrig included in th? 
mod.’! spoc- . All the rnorf, the (Of^y'2 ^"^’^7/2^ configur- 
ation -it 6.13 Me? contributes about 22 times tho 8 ivie? 
coiif iau.-’t ion -'nd we are not including it for the siadlar 


r.iasoti. On tho other hand, within spare III, we are 

2 

Including tho configuration in diagonaliaation 

bocauf.o it has been a convention to include all cotifi- 

0 

durations incl'.,ding upto tho higheft configuration 

wi.tlout realising that even within, space III, tho 

2 

(1Pl/^,)' corf igura clan conxi'-ilah ■’ucl'. less than tho 
leading contributor (lP'^,/2^ ' configuration and o.nce 
-ar '' .all out to determine an e ^'■•‘■ectivG interaction 
only, it might as well be vorv goad aoproxiinatio 1 to 
absorb so snail effects in interact Ion p'-Tarnutors and 
drop this conf igurntlon at all from diagonalizatlon. 

This puts a serious doubt on the wav the model space 
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is arbitrarily chosea oust by tha consideration of 
oinyle particle energy of an orbit anb rest al-1 are 
igri^raf] ,1 iat on the ple^^ that rbov iirrolye high single 

9 

particle excitations. The contribntlon of 

CO iClgurntion to O'*' level is 0.10385 OeT while that 

of (ip-j/p) conf iniiration is 0.0357 WeV. There is no 

2 

reor.on to believe that contribution of, sav, 

2 

conf i .•Tur'''tlo 11 will not be comparable to that of 

conf Iguration .at least. In Table 4.6 we present second- 

orfier contributions to levels other tlian o'" and 1"*". 

The numbers varv over a wide range. Therefore, in 
Figure 4.6 we plot the logarithms of second-order contri- 
butionc to di''f'erent J levels versus the single particle 
excitation enervy of the excited conf iguration. It is 

clear that the second-order contributions for all J 
levels vary erratical3.v and single particle energy of 
a conf iguration can certnlnlv not be the sole criterion 
for Inclusion or exclusion of an excited conf iguration 
in model sp"ce. Sven co-vionsense tolls that instead of 
single particle energy of an orbit, the unperturbed 
energy of an excited two-particlo conf iguration 
would be a more honest criterion for selection of 
model space. This can be seen even from the basic 
philosophy as to why do we have to include configurat- 
ion mixing- 
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T\BL i: 4.6 

'Results fir calculation of second-or'^er contri- 

2 

buttons or configurations above lowest {0frj^2> coiifi- 
nuration, Ivirip; within fp-shell and contributin'^ to 

+ 4-4- 4. 

2 , 4- , 3 and 5 lev 'Is recooctivel"’. 'These results 
arc for = - 70 Pde7, Aq^ = - 40 VieV, = -40 MeV 

and = - 50 MeV. 


E xc i t e d 0 0 nf i 2 ; u- 
ration 

i 1 1 -1 2 ^ 

Hnper turbo cl 

excitation 

energy 

(Me 7) 

<0fy|7| 4 ,l 2 >j 

2 “ 

(Me 7) 

Second-order 

contribution 

(s) 

(Me 7) 

2 

3 

4 

5 

Ofyy'2) 

2 

-0.5504 

0.151470 

(IP3/2)''' 

4 - 

-P .4104 

0.042106 

(lPl/2 

6 

0.2255 

0.008474 

(nf 3/2 ory/ 2 ) 

6.45 

- 0.4251 

0.028017 

(Of 3/2 IP 3 / 2 ) 

0 . 45 

0.3010 

0.010722 

— k 

10.45 

-0.3164 

0.009701 

(Of^/g)^ 

12.90 

- 0.0512 

0.000203 


Oontd 



Contd. Table 4.6 
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2 

3 

4 

5 


0 

£™ 

-0.4011 

0.120580 

ro 

O 

ro 

4 

0. 1762 

0. "^07761 

(Of,-/^ Ofy/g) 

6 . 13 

■ -0.0261 

0.000106 

(0^3/2 ^P3/2) 

P.4d 

0.3009 

0.010714 

(OPg/P^ 

12.90 

-0.059B 

0.000277 

(IP3/2 0^7/2) 

2 

-0.67S1 

0.229900 

(lP^/2)'^ 

4 

-0.40Rn 

0.041300 

(lPl/2 

4 

-0.3432 

0.029445 

(Ofg/g Of^/P 

6.45 

-0.0723 

0.000810 

(OC^/2 IP3/2) 

R.45 

0. 1444 

0.002468 

(0r3/2 1 Pt/2^ 

10.45 

0.1296 

0.001607 

^ 'l ^ 

12.90 

-0.1515 

0.001779 


2 

-0.7547 

0.284800 

(OP3/2 

6.45 

0.1044 

001690 

(Of 5/2)^- 

12.9 

0.0517 

0.000207 



- fO 


V 


-2.0h 


-3.0i~ 


! 


l 


t 




Logarfthrn cf hccnrvi-crclcf i'cnt/ i'jufions (55) tcwnr^is different 


J-ieV'.dfe fjfisinn if err; 

of Ide o-npcr1orbo;d encitatlnn 


crin-fiooration f.dr'tted as c func 
enenpy cf tne two-particlti ror 
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In view of these doubts it becomes very 

necessary to see if saturation in respect of model 

space dimensionalit-'p is ever achievable. T?or this 

purpose \¥e neeb to calculate the second-order contributio: 

of con'^igurations beyond fp-shell. However, correct 

single pa-''’ticle e’lergies in this region are not 

estoblisheb. There‘'^ove, we adont the level scherno as 

prooof'ed by I'layer and Jensen^"^^ and arbitrarily assign 

single pa tide energies. Those are shown in Table 

d7, kt■^ such, the numbers which we would arrive at 

should not be taken too literallv. Yet t^'ls calculation 

is expected to give an apnroximete order-ot-mognitude 

ostim^-^te of e-ffect of high - lying configurations on 

+ 

low-lylnv levels. The results for 0 and 1 levels 
are tabulated in Tables 4.8 .and 4.9 respectively. The 
loga.ritlims of the numbe'^^s are plotted versus configur- 
ation number in Picu-'-'cr. 4.7 and 4.8 respectively. It 
is vi’pv do'' thut rh'- contrib'^t’d' is keep on oscili ating 
instead of monotonical L'’.’ deoreasinv as one hopes. 

Yg did not go b-'vond orbit but the trend of 

gf-aphs raak'^s ifc dear tha*" convergence in level 
energies is not achieved upto SO Mo7 conf iguration 
•■^t least. Although no strict relationship could ben 
established but, generally, for an e.xcited |.i.].l 2 > 
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T.A3L5 4.7 

SiTitQrle particle levels' seq.ueiace (ta^cen trom 
M-'iye''’ and Je'iisen ' ) used fer calculntion of second- order 
CO itrlbutions from higher conf igur'-^tions. dxcept for first 
four conf igurations , single particle energi ns haAre been 
arbitr'^rilv assigned. 


S . "To . 


Orbit 


“Sin'^’e" ParticTe" 
(MeV) 


Enex gy 


1 Of^/2 0 

2 IP3/2 2 

7 V 1/2 4 

4 '^'^5/2 6.45 


Og 9/2 


6 


15 

7 

^'^5/2 

17 

p 

^3 /2 

19 

9 

2^1/2 

21 

10 

Ohi 1/2 

23 


Contd - 
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Co'itd. Table 4.7 


'lo . 


Oi-b It 


Single Particle Energy 
(Me 7) 


11 

Ohg/2 

30 

12 

^■^7/2 

32 

13 

'-^5/2 

34 

14 

2 p ^/2 

36 

15 

2 Pi /2 

33 

16 

Oil 5/2 

40 
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T^.BLB 4.8 

Results for CTlculetion of seco id-orrier contribut- 
ioTi(s) of all conf iffur-^ttons lyiag upto 

contributing to 0 l>"'V>j!l arisine fro’''-' (OfY/2) confi- 

guration. Calculations have been done ^^or 4 .^q = -4-1 MeV, 
= - 50 Me7. 


Con^i- 

guration 

Index 

Two-part- 
iclc confi- 
guration 

ipp > 

Unpertur- 
bed enorgv 2 

of the two- <0fy Iv 1-1-1 

parti, c].e 

conf iguration , 

(Mel) (Me 7) 

Second-order 
c )ntribution 
2) + (.") 

(Me 7) 

1 

2 

3 

4 

5 

1 

2 

( 1 P*^ y/2 ^ 

4 

-1 .1396 

0.3247 

2 

( 1 PV 2 )" 

8 

-0,5348 

0.0357 

3 


12.9 

-1.1576 

0.10385 

4. 

y'p) ^ 

18 

3.52341 

0.68967 

5 

(0^7/2)^- 

30 

1.02577 

0.035072 

f. 

^"^'^’7/2 *^^7/2^ 

32 

-1.15856 

0.06349 

7 

(ld^/2)‘" 

34 

1.076165 

0.034062 

R 

(1^3/2)'^ 

38 

0.423915 

0.0047289 

9 

(2P3/2 ''P3/2) 

38 

-O.23O86O 

0.014026 

10 

(If 5/2 0'f5/2^ 

40.45 

-0.182364 

0.00082218 


Contd. 
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Contd, Table, 4.8 


1 

2 

3 

1 1 

(2s^/2)^ 

42 

12 

(2Pl/2 ■'Pl/2) 

42 

1? 


46 

14 

(0hy/2)2 

60 

15 


64 

16 

(1 1-5/2 

68 

17 

(21)3/28 

72 

18 

(eni/g)^ 

76 

19 

(01,5/2)2 

80 


4 

5 

-0.289824 

0.002 

-0.045978 

0^000050334 

-3.31757 

0,23925 

0. 11 258 

0,000211 23 

-1.69937 

0.045128 

0.06972 

0,000071481 

-0.14019 

0,00027295 

-0,636875 

0.0062078 

1.594075 

0.031765 



273 


TA.BLE 4.9 


Results !^or ce.i,ci:iT qtiou q-p second-order 

contribution (S) of all conf i^ura tions lyine upto 

+ / \2 
orbit contributin.'^ to 1 level arising t'rom (OfY/ 2 ) 

C'.'infiguration. Calculations have b ■•''n done tor Kqq = -70 MeV . 

= - 40 Me 7 . 


Oonfigu 

ration 

Inde\ 

Two- particle 
confinurat io'i 

1 ;) 1 i 2 > 

linpertui’be'd 
energy of 
ttie two-par- 
ticle config- 
uration 
(Me7) 

<0f^|T| 1^: 

(MeV ) 

Second-order 
. V. , contribution 

2 1 ^3) 

(MeV) 

■~T 

2 

■y 

4 

5 

1 


4 

-0.8280 

0.1714 

2 

(lp ^/2 ''P 3 / 2 ) 

6 

-0.3444 

0.01077 

3 

(Of 5/2 Of 7 / 2 ) 

6.45 

-0.9198 

0.13117 

4 

(1Pl/27 

8 

0.2198 

O.OO 6 O 4 

5 

( 0 fr ,/2 1P3/2) 

8.45 

0.4382 

n . 02272 

6 

(Ofj/j)® 

12.9 

-0.4124 

0.0131B4 

7 

P 

(OS 3/2 ^ 

18 

3.28617 

0.59994 

8 

(Og ^/2 ^^ 9 / 2 ^ 

24 

0.89231 

0.033176 

9 

( 0 ^ 7 / 2 )^ 

30 

0.35501 

0.0042011 


Gontd. 
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Contd. Table 4*9 


1 

2 

3 

4 

5 

10 

(011^^2 ^--7/2^ 

30 

-0.66186 

0.014601 

11 

(IO5/2 ^§7/2) 

32 

0.19106 

0.0011407 

12 

(107/2 O frj /2 ) 

32 

-1. 50360 

0. 070661 

13 

(Idj/g)*'’ 

34 

0. 88 9 49 5 

0.023271 

14 

(lf^/2 0^7/0) 

34 

-0. 112685 

0.00037345 

15 

(ld-^/2 1^5/2) 

36 

0.300596 

0.0025099 

16 

(l'=‘ 5/2 ■’P3/2) 

36 

-0.093983 

0.00024536 

17 

(id^/g)^ 

38 

0.035665 

0.000033474 

1 B 

(2P3/2 '>P3/2) 

38 

-0.214010 

0,0012052 

I'j 

^“^^7/2 ^^^5/2^ 

38,45 

0.1126B5 

0.00033022 

20 

(23 1/2 '1^13/2) 

40 

0. 180978 

0.00081893 

21 

(2P3/2 IP1/2) 

40 

0. 0091583 

0.0000020969 

22 

(2v-|/2 ’^^3/2^ 

40 

-0.0091583 

0.0000020969 

23 

(105/2 ■'^•^5/2) 

40.45 

-0.27'7804 

0.001908 

'r 

(23^/2)^ 

42 

-0.189732 

0.00085714 

25 

(2Pl/2 ''Pl/2^ 

42 

0.00299314 

0.0000002153 

26 

(2P3/2 0%/2^ 

42.45 

-0.26451 

0.0016482 


Contd 
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1 

> 

cl 

1-y 

:/ 

4 

5 

27 


46 

- 3.700176 

0. 29762 


On ^ ^ ) 

53 

-0.470480 

0.0041766 

2 ') 

(Oh,, /,,)■- 

60 

-1.543502 

0.039700 

•50 

('‘ 7/2 '’'’V 2 ' 

62 

0.029784 

0.000014309 

31 

("'7/2’“ 

61 

-1.911334 

0.057079 

32 

(' d /2 "’ 7 / 2 ( 

66 

-0.114136 

0.00019739 

33 

( 1 ' 5 / 2 )" 

68 

-0.710628 

0.0074264 

3 1 

(2li.,/g Ifj/g) 

70 

-0.264500 

0.00099945 

3 '; 

( 2 p ,/22 

72 

0. 297558 

0.0012297 

36 

(2r,/2 2p,/g) 

74 

-0.5 66 220 

0,0043326 

37 

{2Pl/g)'' 

76 

0.51134B 

0. '■'■034404 

3 ^ 

(Oj.i,/g)2 

80 

1.399082 

0.02447 






fi'j 4,7. LcgvJrithrr record -orn-^r cofitributions (5) towotds O^b'vti 

2 

■'urrsinc] fr'srr (c^ ■'•i.'rjftqumtior plotted as a funclicr 
cf intkx number c-f the excited two-partj'clc cofvtiguroti'- tt 
(Cf. table ■ 4.8)^ 



legorithm of second-order contribuiions (S) towards levc 

, ' 

arising from (of 7 / 2 )“'' configuration plotted as c function of 
•f'cto number of the, excited two -particle configuration 
{Cf,., table , ,''4.9 ) ' ■ ' 
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conf i ^'ir=itio)i, the co itribution is small for small 
’h "** *2 40 atsY 42 MeV conf i^ura tions for 

“f* 

1 Ilh/, 1. jufTh cant rihut ion;; of coth tho 42 Me? 

CO afj 'nj r- tm ou'-' are .le.alio' ..blo yet such a l'^r,o;e differeace-- 

■ 1 r-ir.j-, ,r of 4000 !:■. no h imd--;rstnuc1a'jl8 on purely single 

p'vr'tip'’ ' ori'.'rgies ' considv'rat''.ons . The contribution of 

th ■' ;■'?} ''T conf iguration to 1 leyel is 

2 4 4 7 Mmr am’ in view of practical limitations it can 

hi.> ne*''!. ctod in comparison to 0.1714 Maf contribiitio:i of 

ttie 4 'to? i;orif ivura t;ion (IP 3 /;;) wo should view it 

i.'i (;!p'''n*iriu,o,ii bo t’^.e eontribut ion O.OO 6 O 4 Me? of the 

2 

8 ’V'V t' 0 !i'''iyu-!'ation (ip^y'p) w'>icb pe have alway,a been 


i nt’lU'! ’ UP ! The coiitri 'out ion O.5OO94 Me? o*^ 18 Mo? 
t;n nf i 'Curat ion level is much more tlian 

any 0“ th'' fn-s'h,ell configurations, ’‘'a, theref oi-e , 
ar.-rue that tho (lP'[/'2^*^ configuration con bo neglected 


P 

i,;i favour 0''' the (01--|3/2) configuration and similarly 
the ('"*79/2^'' ^^^^11/?^' configuratio ns should also 

t:)C e,^'| 5 l ic jT.lv Included in diagonalination. We present 
soni'-; o*’ such contrasting numbero! in' Table 4-10 for 


0'*’ level . 


We have done these cnlculations for a randomly 
chom.n not of parameters and for one nucleus and one 
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TASL9 4.10 

order contributions of some of the excited 
I'v'u i.'Ti) ti one contributing to f lercl arising from 
(f* f 'j, CO nf Igurntio n piched up so as to sHow fiat the 
uruji’ct u rb ::d ■niergv of tlie two-porbicl? c mf Ihuration cannot 
!a; tlu' ool- criterimi for incl-i.<sion of a tv^o-partiole 
t>(Vi I'i -ui,* ■ f; o:i in moriel sp'ice. 


U‘ uKCltod Sto.O 
'o-i'ort' Lcl e con- 

Unperturbed energy 

0 " the tv'o-T)o,rticie 

Second-order 
contribution (S 

• tu r-'it i un 

conf iaurotlon 
(MeV) 

(MeY) 


42 

0,002 


42 

0.000050334 


64 

0.045128 

1 / 

68 

0.0000714B1 

■, ^ 
l/P 

76 

0.0062078 


80 

0.031765 

°r'i/2 

o o 

0.035072 

0.031765 

0i\t^/2 





280 


Wight suspect if it is uot going to be the case as 
faccfl by Oohen et . al merelv sn accident as proved 
later i'.f G-upt-^ and Trainer =’nJ 07 jiia'-..! and Unna. 

In the first plac.e, even if such a thing happened on 
actual calculations in other shells, our results still 
stand vuiere thev a^’e - that the nonconvergenco may or 
may not come. Chancea se!‘;m to be slim that such a 

thj.ag would ha'caen. Just to check, we calculated the 

/ \ P 4. -j. 

effect of (0fyy2/" conf iguration on 0 and 1 levels 

2 2 
of (Od^y'2^ conf iguration and that of (Ogqy'2) confi- 

gura.tion on 0 and 1 levels of (lp^/2) configuration. 

Same parameters as those used here were taken. The 

Of 7/2 orbit was assumed 12 Me 7 above the /2 

and thi^ 0 ^c 5/2 assumed to be 6 MeY above the 

1p^/2 orbit. The results are presented in Table 4 -. 11 . 

42 

The results onlv substantiate our findings for the Sc 

nucleus. Ohoico 0 ■' anv other radial shape or the 

potrmtial o'ill certainly not make any drastic change in 

rf'si.ilts. The strengths were also arbitrary and there 

is no reason why the general trend of ros'ilts should not 

hold for otben oar.ameter value.a, Sven the use of 

r^aalistic interact 'ons is not expected to make the 

1 47 

things go the other waj^. "^ork of Joshi has shown 
that realistic potentials can be satisfactorily simulated 
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TOL'g 4.11 

? 

Second-oriler contributions of (0fy^2) 

(Ogq/g)'^ configurot loni,.f to (Od^^y.^)^ cai (Ip^y2^"^ 

conf iff ur'^t ions for O''’ and l"^ levels. Calcul'-itions 
Ir’. ffo been periorffu.-d. A-| q = -40 Mef, = -50 MeV, 

/Iqo = - '70 an-T = -40 IvIeV. 


Matrix olefn.mt 


(Mel) 


'TTnpc rTiirb 3T"'" * 
e tier op'' of the 
two-particle 
configuration 


Sec ond -order 
Contribution 
(S) 

(Mel) 


2 

<0dq 


^■■^7/2 ^Q+ 

3.60574 


24 


0.54169 


<0ch^/2l''l ^*7/2 >1 

- -S. 279268 


+ 


24 


0.44,8097 


>0"- 

1 . 37289 


12 


0.15708 


:ip^/2lll Ogq/2 >^ + 


12 


0.08125 


= 0.93742 
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by phenomenological potentials and we can reali'^je 
that in a stadv like the present onn what matters is 
the parameters ’ values and not the wav the^^ have been 
derived - they may h'^- phenomenological e'^f'ectime inter- 
-^ction derived to fit spectra or phenomenological potentials 
eauiyalc'iit of realistic interactions. Conversely, the 
trend is also expected to hold i^ the realistic inter- 
action itself is used. 


IT. 4 COIICLITSIOTS : 

Our study shows that a Berber type effective 
potential can bo taken for T = 0 interaction but not 
f'jr T=1 states. Further, it can be possible to rouahly 
ectimato beforehand th-^t some of the conf igura+'ions , 
whicls wi:i might otherwise ho going to include, really 
contribute very littlfD and they can as M<iell be dropped 
from calculation. Lil':.:; this by making a preliminary 
invostivation, we can get a physicallv more meaningful 
o f f 'i c t i '/ ' i nt e ra c t i o n by inc 1 ud tng an equal num b er 
of more important effects in place of those dropped 
even without adding to compr'tation'^l difficultiee. 
Piirthor, it is not legitimate to arbltraril-”- choose a 
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model spoce .iust on the basis of single particle 
energy of: highest orbit in the model space. Instead, 
our calculations sho'.v that even a liighlv excited two- 
part? cle configuration can contribute more stronvly than 
a comparatively much low excited one. However, many of 
the configurations in going upto the highest conf i.guration 
also contribute very weakly so that they are better 
dropped. 
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AP PENDIX A 


W.^ve function oP a Giaityle nucleon o-p mass ra, moving 
in a harmonic oscillator pot'-ntial, is given bv 

R -j (r'l 

'^nO. ra --"T 

Th‘^ novraal i^.ed sDl'itton, corresoonling to energy 
eigenvalue 


-n I 

is given by 


(2n + ^ + -I CO 


(A. 2) 
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nS, 


(r) 




f^i(~2l V 1 ) !~! 


2 (2-^.+ 2n + 1 )! ! 

J ^ n! 

(^.3) 


(2-yr^) 


wh'/)ve , 




(x) 


n 

= r 

k=0 


(-f s'" 


I 

U) 


(2iL+1 ) ! ’ 

(2 £+ 2 k +1 ) ! ! 


(11.4) 


and 


mco 


(A. 5) 


The wa'/e function for t'vo nucleons relnUve 
motion is also given by above expressions provided V is 



. . . i' 

replaced by V = where /A. is the reduced mass 

o'’ tie pnrticles. For two nucleons, /-K = ~ . Hence, 


hi h) 


X -n+2 


IVX ^ (2 i+2n +1 

f% )_(2t +1 ) M ] 2 n! 


1 ? 

r” e v^.(Vr^) 


(A..6 ) 


where V is same ao defined, above and r = I r^ - ^ 2 ], 


4 properly antio'rT^'oetrlced and norna^ iced two-particle 


wave function in 1.i - coupling is glAren bv 


11 ^ , lip £ 2^2 = > 


T 


N hn^ £ .,D^» Ho ^2'’2* > 

:]■ +,ip+1+J+T . n 

(-) ' ^ Inp-tp.ip.H^ £^,1^ :JM/1’M^ > 


(A. 7) 


v/here the conyention is maintained that angular momentum 
of" particle 1 is counled io annul t momentum 0 ^ particle 2 
(in this order) and H i-". given by 


(2+2 ^ 6£ 0 6. . ) 


(A-B) 



i- 3 


The expression for properly antisvi.iiinetri.zeci 
end noryn'i.irjed metrix element of t#o-bod'/" interaction 
V betv/een two such stetes is ffiven bv 

^ 1 '^1 ’ -^232 '^'1 1 ’ '^ 2^2 ^2 

= ?'Ti' [<ir ■■ii 

.1' + ,1A + J+T , , 

+ (-) ^ <.1^,l2:JTlV| 1^' :JT> 

(4.9) 

where otter quantum numbers have bee i drooped. 


The expression for nntrix element can also be 
miven in a form where matrix elements appearing on' right 
side of eq- (4.9) are calculated in L S coupling 
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2Tf'r 


i: 

L3 
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2 
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X 



^2 ^ 2 ' ^^1 1 I "^1 ’ ^ 2 ^ 2 * 


^2 + 1+S+T 


<n^ , Ug -^2 :LMj^| Vj n2 ^ 2 ’ 




( 4 . 10 ) 



A- 4 


However, ( 1 . 10 ) Hoida strictly for C'^nti'al ■ forces 
wliicH conserArc L -^nc! it is ■assnme''' th spin, and isospin 
deifeiid'-'iioe of 7 has bvcen c.-tLculatcd sep irately. 

In th} method of Slater integrals, the direct and 
exchange te-^ms ii eq. (4.10) are separately colcolated using 


k;=o 




(A. 11 ) 


wliere 
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i ( 2 -^,, + 1 )( 2 -^ 2 +'' ^2 ^\ + l)(2 ^ 2 +^) 
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1 ^ f- L 
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y 1 2 
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1^2 
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(4.12) 


and 

r= r / R 

0 b 




I 

(4.13) 
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2 k' + 1 
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R^, (r2) ^^1 ‘^'^2 d(Cosc<J.|2^ 

2 2 


(A. 14) 



where , 


CosCO^ 2 



Definin'? , 


r 




md 




0 0 n a 


r . R 
rR 


Hence , 


Cos 



"*2 


(A. 15) 


CA..16) 




The Tjegendre function. CCosiiJ.j2) 7 e expresssed 

es e polynomiel is CosiC^^- ^ volumt-* element 

dr.j 1*2 dr2 d(OoBCO^p) in (r..j, representation 

o 2 

becomes dr R dR d (Cos a) in (r,R,a) representation. 
Hence, rinolly, the inherrend in en. (A. 14 ) becomes a 
f'inct ion oP r,R and a. Thu.a, 
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r f / 

■ 0 0-1 




R„ n (ro)R, 


■ n , 5 .r 


n ^ 


^1 


T / 


\/(^)P^(CosC0^2) dr dR d (Cos a) (A. 18) 



CB Lculnt -cl usinfY tho iTiethod of Sister ioteffrals. 


<n.^ 3^2 ^2 '■ ■ I "1 ■^1’ ^2 ’^ 2 ' ^"'■L ^ ~ 

( 'V. 24) 

where, f2 is c-'lcul'^ted from. (.1,12) s’ld the Sinter 
0 

inte;Trsl 7'" is cnlnulited from (\.13) 




2 2 

r!; R 


2 "nj “"2 


dr. 


(/I. 25) 

ThoAr .-r^re tabulohed in the append-i x C. In eq. ( 4 , 23 ) 
% is the strenrt'o p'='rsmeter. 


Matrix elements of the spin-dependent 6-funotion 


,v'o roe 


V " '^0 \J ■■ + a ( ^2^ 6(r.^ " ^2^ (4.26) 


BT’o '■^ivon hv 


<n. JT | n'!' 1 , n^-d'oli 

^ Q I I 1 ^ C. C I ' ' d 

% O 


;J1 > 


'\S 


2 It 







Expression Contd- on next page 
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( 4 . 27 ) 


where, the notn+ion [(_xl = x +1 and Yq^ is stren^-yfch 

in the particular spi i-state S. For two ideutlcal nucleons 
T = 1 and hence, via a 6-function force, the antisvirroetry 
condition requires that they can intvsract in spin singlet 
state onlv and tVv' expression final 1v rednces to 
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"I .1 n't ^ t do 1 1 
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^L^.nu4L^ri.o. 
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( 4 . 28 ) 


The radial integrals are tabulated in appendix B, 



In the method of Moshinshy brackete , 





B 


ni 


n 


. <■» 
i 3 




n'i' n' I' 
^ 1 1 , "2 2 

B 


n ■ 


^ ^ L 


5 L 

Is 


■P /V '» 

( <(nJ2,S)J'|vi (n^ I', S)J' 
J JM 


For- central forces it reduces to 
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<n^\Y\n'Jl' > =I 




Z B(nt, n'^',p) I (4.01) 


(£.+£)_<p_< ■^(-^+^)+n + n 


(4.32) 


The inte.-ir^l Is snme ps 

P Up 


Vi;),(,) IopIi! 


r dr 


(4.33) 


Tn ffpner'Hl, 


cx.) / N 


v(r) r dr (4.34) 


where the radial wave funotlons are for relatl.?'e motion. 


Phnr, (4.33) reduces to 


,pP+2 


p + 2 °° 
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Ylt ( 2p +1 ) ! ! 0 


\/'(r)clr (A. 35) 


?hlch is conai'’tc -it with (4.10). 


The raiMal i.'ii.e^ral I 


n , nt 


is also denoted by I. 


liatrix elements O'” ti\e tvvo-narticlo spin-orbit 


V(r) g. 


Interacti on 



k. 1 1 


pre given by 
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v/here P ■= 2n +£,+ 2:T + L 
with the condition 


1 ?. 




X X , A i ' 




Motrix elements of the tensor -force 


( P v(r) X (e, ([))-:x,. 


( 4 ,? 8 ,) 



4. 12 


where is the Racah tensor whose components are 
the spherical harmonics y2^p(0,(t^') 
two R'.,cah t-nsm- whose m - 0 coraponent i 


ani X2 is the rank 
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(4.39) 


are cs'cnlated from 
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1 

with the condition 


^ = A,A± n X± 2 
.£^= 5 , 5 + 2 
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In the formalism of second quo nti^at ion j confinin'?; 
t;: t' interactions onlv, th-‘ Ha\iiiltonian is defined 

in, 
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a [3 Y 6 
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(A. 41 ) 
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1 1; .vi'i'C' ri? I' (■'^rorico e'lorsv 


’ -ap 

r. !•> ' 1 t'l ni qn - f'lo Ylneti.o energv, 


if-i tlie one-body 
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(A,. 42) 


f»*! ' \»4 unantlPVTnrTietr]>.ed ip^trix elements 

'/ l.M't v.'een the two-perticle states j ctp > anri 1 Y 

iiof. ds 0 innlHcr a closed shell nucleus whose yround 
:Y )*'.!: h-‘ 'ieentu'd bv | 0 > r?ucb f'lat all states below a sharpl-'" 
j, i’,- 4 p'iii m 'inv.'Ptnf' ane com{>'ctelv filled and those 
■p. -V',-' comuletel'/ empty. Tficn, rel-ilive to the ground state, 
the (!•') ^n' th"’ hole (h) states are defined through 

n I fl > - 0 

nj t 

1 0 > = 0 

where a and a"*" are respectively the particle annihilation 
anil creation operatorrs » t:he subscript collect- 

ively denotes the quantum number o.e state on which they act. 
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The creation of a hole state anfl a particle state 
is then defined by respect i rely , 


■’n". >=(-)" !o> 


I :lp”'p > - 3 jpHp I 


The particle and hole states defined like this behave 

sirailarl^r under rotationsy and hence can be coupled 

throuyb C.G. coefficients, where- -as the states ™ ln> 

''h\ 

^ j 

and a. ^ 0> do not. The somewhat different structure 

•'p p 

of the hole state creation operator from the corresponding 
particle state creation operator arj.aes from th"' fact 
that relative to a ^^round state of total angular momentum 
'r,ero , a hole state I ^ vrlll be created when a 
particle ■f'rom the state 11]^“’^]^ > is destroyed. Thus, 


= I (y - "'h^ ' ^ 

where the symbol signifies vacancy in the particular 

nu.antum e-tate. In shell model description, an excited 
state in a closed shell nucleus occurs when a particle 
from th'-' Fermi sea is lifted to one o^' the empty orbits. 
Therefore, a total angular momentum coupled hole -particle 
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ev.cited state will be obtsini^d as 
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( . A? ) 


''’'e make a term - by - te!"ra calculation, of matrix 
olemento ot P between two sucb h~p states. 
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If iso3pln is slso included 




y ( 2."'+1 ) ( 2 T ' +1 ) 


•I f 


.1h- -ip. J' 


C 1 
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1 ; 
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h'-p 


IS 


(A. 47) 


r *1 e detniled derivation where all the terms are 
retained on the ri^ht of eqs.(A.45), (a . 46 ) and (\.47), 
tl'C'-’o coffu-? terms which represent the core energy and the 
i nte'f’act l 0 "i the hole and the particle v/ith the core- 
The litter Interactions can be ignored it one tahes energies 
of particle levels from e^xneriment since experimental 

numbers absorb all effects to all. orders. In a calculation 
^her’C one needs to calculato the total binding energies 
o'* tfv le/els, one also renuires the core energy. However, 
in cases where one is required oniv to calculate ths 
excitation energies of levels relative to the ground state , 
tha co'^e energy term simply drops out. 
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OTEWDIX B 


Tnhl>'v’. O'” the inte.rr.'^ls 
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\11 oth'5r I'T-tegra^-U needed in eny of the 
ei-eva meMtinned shelle and not written here 
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APP3NDIX C 
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